| summary.jointPenal | R Documentation |
This function returns hazard rations (HR) and its confidence intervals.
## S3 method for class 'jointPenal'
summary(object, level = 0.95, len = 6, d = 2,
lab="HR", ...)
object |
output from a call to frailtyPenal for joint models |
level |
significance level of confidence interval. Default is 95%. |
len |
the total field width. Default is 6. |
d |
the desired number of digits after the decimal point. Default of 6 digits is used. |
lab |
label of printed results. |
... |
other unused arguments. |
Prints HR and its confidence intervals for each covariate. Confidence level is allowed (level argument).
frailtyPenal
## Not run:
data(readmission)
#-- gap-time
modJoint.gap <- frailtyPenal(
Surv(time, event) ~ cluster(id) + sex + dukes +
charlson + terminal(death),
formula.terminalEvent = ~ sex + dukes + charlson,
data = readmission, n.knots = 14, kappa = c(9.55e+9, 1.41e+12)
)
#-- calendar time
modJoint.calendar <- frailtyPenal(
Surv(t.start, t.stop, event) ~ cluster(id) +
sex + dukes + charlson + terminal(death),
formula.terminalEvent = ~ sex + dukes + charlson,
data = readmission, n.knots = 10, kappa = c(9.55e+9, 1.41e+12), recurrentAG = TRUE
)
#-- It takes around 1 minute to converge
summary(modJoint.gap)
summary(modJoint.calendar)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.