View source: R/summary.joint.R
summary.joint | R Documentation |
joint
object.Generate summary of a fitted multivariate joint model.
## S3 method for class 'joint'
summary(object, ...)
object |
a joint model fit by the |
... |
additional arguments (none used). |
Object of class summary.joint
.
James Murray j.murray7@ncl.ac.uk
joint
and joint.object
# Simple univariate on log(serum bilirubin) ----------------------------
data(PBC)
long.formulas <- list(
log(serBilir) ~ drug * (time + I(time^2)) + (1 + time + I(time^2)|id)
)
surv.formula <- Surv(survtime, status) ~ sex + drug
fit <- joint(long.formulas = long.formulas,
surv.formula = surv.formula,
data = PBC, family = list("gaussian"))
summary(fit)
# Bivariate with a dispersion model ------------------------------------
PBC <- na.omit(PBC[,c('id', 'survtime', 'status', 'sex',
'drug', 'platelets', 'albumin', 'time')])
long.formula <- list(
platelets ~ time * drug + (1 + time|id),
albumin ~ time * drug + (1 + time|id)
)
surv.formula <- Surv(survtime, status) ~ sex + drug
fit <- joint(long.formula, surv.formula, PBC,
family = list("negbin", "gaussian"),
disp.formula = list(~time, ~1))
summary(fit)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.