Test the Proportional Hazards Assumption for Relative Survival Regression Models
Test the proportional hazards assumption for relative survival models (rsadd
, rsmul
or rstrans
)
by forming a Brownian Bridge.
rs.br(fit,sc,rho=0,test="max",global=TRUE)
fit |
the result of fitting a relative survival model, using the
|
sc |
partial residuals calculated by the |
rho |
a number controlling the weigths of residuals. The weights are the number of individuals at
risk at each event time to the power |
test |
a character string specifying the test to be performed on Brownian bridge.
Possible values are |
global |
should a global Brownian bridge test be performed, in addition to the per-variable tests |
an object of class rs.br
. This function would usually be followed by both a
print and a plot of the result. The plot gives a Brownian bridge for each of the variables.
The horizontal lines are the 95
of the Brownian bridge
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.
data(slopop) data(rdata) fit <- rsadd(Surv(time,cens)~sex,rmap=list(age=age*365.241), ratetable=slopop,data=rdata,int=5) rsbr <- rs.br(fit) rsbr plot(rsbr)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.