summary.jointPenal {frailtypack} | 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%. |
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. |
Prints HR and its confidence intervals for each covariate. Confidence level is allowed (level argument)
frailtyPenal for Joint frailty models
,
plot.jointPenal
,
print.jointPenal
,
## 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)