frailtypack (version 2.2-13)

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 for Joint frailty models, plot.jointPenal, print.jointPenal,

Examples

data(readmission)

## Gap-time ##
modJoint_gap<-frailtyPenal(Surv(time,event)~cluster(id)+sex+as.factor(dukes)
		+as.factor(charlson)+terminal(death),
		formula.terminalEvent=~sex+as.factor(dukes)+as.factor(charlson),
		data=readmission,n.knots=14,kappa1=9550000000,
		kappa2=1410000000000,Frailty=TRUE,joint=TRUE,recurrentAG=FALSE)

## Calendar time ##
modJoint_calendar<-frailtyPenal(Surv(t.start,t.stop,event)~cluster(id)+sex
		+as.factor(dukes)+as.factor(charlson)+terminal(death),
		formula.terminalEvent=~sex+as.factor(dukes)+as.factor(charlson),
		data=readmission,n.knots=10,kappa1=9550000000,
		kappa2=1410000000000,Frailty=TRUE,joint=TRUE,recurrentAG=TRUE)

						 
						 
# It takes around 1 minute to converge  

  summary(modJoint_gap)
  summary(modJoint_calendar)