Learn R Programming

frailtypack (version 2.7.1)

summary.jointPenal: summary of parameter estimates of a joint frailty model

Description

This function returns hazard rations (HR) and its confidence intervals

Usage

## S3 method for class 'jointPenal':
summary(object, level = 0.95, len = 6, d = 2, lab="hr", ...)

Arguments

object
output from a call to frailtyPenal for joint models
level
significance level of confidence interval. Default is 95%.
d
the desired number of digits after the decimal point. Default of 6 digits is used.
len
the total field width. Default is 6.
lab
label of printed results.
...
other unusued arguments.

Value

  • Prints HR and its confidence intervals for each covariate. Confidence level is allowed (level argument)

See Also

frailtyPenal

Examples

Run this code
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)

Run the code above in your browser using DataLab