# \donttest{
# Bivariate fit on PBC data -----------------------------------------
data(PBC)
# Subset data and remove NAs
PBC <- subset(PBC, select = c('id', 'survtime', 'status', 'drug', 'time',
'albumin', 'platelets'))
PBC <- na.omit(PBC)
# Specify bivariate fit
long.formulas <- list(
albumin ~ time*drug + (1 + time|id),
platelets ~ time * drug + (1 + time|id)
)
surv.formula <- Surv(survtime, status) ~ drug
fit <- joint(long.formulas, surv.formula, PBC, family = list('gaussian', 'poisson'))
fitted(fit)
# }
Run the code above in your browser using DataLab