Graphical Inspection of Proportional Hazards Assumption in Relative Survival Models
Displays a graph of the scaled partial residuals, along with a smooth curve.
## S3 method for class 'rs.zph'
plot(x, resid=TRUE, df = 4, nsmo = 40,var,cex=1,add=FALSE,col=1,
lty=1,xlab,ylab,xscale=1,...)x |
result of the |
resid |
a logical value, if |
df |
the degrees of freedom for the fitted natural spline, |
nsmo |
number of points used to plot the fitted spline. |
var |
the set of variables for which plots are desired. By default, plots are produced in turn for each variable of a model. Selection of a single variable allows other features to be added to the plot, e.g., a horizontal line at zero or a main title. |
cex |
a numerical value giving the amount by which plotting text and symbols should be scaled relative to the default. |
add |
logical, if |
col |
a specification for the default plotting color. |
lty |
the line type. |
xlab |
x axis label. |
ylab |
y axis label. |
xscale |
units for x axis, default is 1, i.e. days. |
... |
Additional arguments passed to the |
Goodness of fit: Stare J.,Pohar Perme M., Henderson R. (2005) "Goodness of fit of relative survival models." Statistics in Medicine, 24: 3911-3925.
Package: Pohar M., Stare J. (2006) "Relative survival analysis in R." Computer Methods and Programs in Biomedicine, 81: 272-278.
Relative survival: Pohar, M., Stare, J. (2007) "Making relative survival analysis relatively easy." Computers in biology and medicine, 37: 1741-1749, 2007.
data(slopop)
data(rdata)
fit <- rsadd(Surv(time,cens)~sex+as.factor(agegr),rmap=list(age=age*365.241),
ratetable=slopop,data=rdata,int=5)
rszph <- rs.zph(fit)
plot(rszph)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.