plot.frailtyPenal {frailtypack} | R Documentation |
Plots estimated baseline survival and hazard functions from an object of class `frailtyPenal'. Confidence bands are allowed.
plot.frailtyPenal(x, type.plot = "hazard", conf.bands=TRUE, ...)
x |
Object of class frailtyPenal (output from calling frailtyPenal
function). |
type.plot |
a character string specifying the type of curve. Possible value are "hazard", or "survival". The default is "hazard". Only the first words are required, e.g "haz", "su" |
conf.bands |
logical value. Determines whether confidence bands will be plotted. The default is to do so. |
... |
Other graphical parameters |
Print a plot of class frailtyPenal
data(readmission) # Let's compare shared frailty model with Cox proportional hazards model mod.sha<-frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id), n.knots=10,kappa1=10000,data=readmission) plot(mod.sha,type="surv",conf=FALSE) mod.cox<-frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id), n.knots=10,kappa1=10000,data=readmission,Frailty=FALSE) lines(mod.cox,type="surv",conf=FALSE,col=2) # Stratified model mod<-frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id)+strata(sex), n.knots=10,kappa1=10000,kappa2=10000,data=readmission) plot(mod) # no confidence bands plot(mod,conf.bands=FALSE)