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.