summary.jointPenal {frailtypack}R Documentation

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


## Not run: 

  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)
## End(Not run)


[Package frailtypack version 2.2-12 Index]