library(pammtools)
data("patient")
ped <- patient %>%
dplyr::slice(1:100) %>%
as_ped(Surv(Survdays, PatientDied)~ ApacheIIScore + CombinedicuID, id="CombinedID")
pam <- mgcv::gam(ped_status ~ s(tend) + ApacheIIScore + s(CombinedicuID, bs="re"),
data=ped, family=poisson(), offset=offset)
gg_re(pam)
plot(pam, select = 2)
Run the code above in your browser using DataLab