rstrans {relsurv} | R Documentation |
The function transforms each person's time to his/her probability of dying at that time according to the ratetable. It then fits the Cox proportional hazards model with the transformed times as a response. A possible side effect are the transformed times.
rstrans(formula, data, ratetable, int,na.action,init,control,...)
formula |
a formula object, with the response on the left of a ~ operator, and
the terms on the right. The terms consist of predictor variables separated by
the + operator, along with a ratetable term. The ratetable term
matches each subject to his/her expected cohort.If the variables are organized and named
in the same way as in the population tables, the ratetable term can be omitted.
The response must be a survival object as
returned by the Surv function. The time must be in days.
|
data |
a data.frame in which to interpret the variables named
in the formula .
|
ratetable |
a table of event rates, such as survexp.us .
|
int |
the number of follow-up years used for calculating survival(the rest is censored). If missing, it is set the the maximum observed follow-up time. |
na.action |
a missing-data filter function, applied to the model.frame,
after any subset argument has been used. Default is
options()$na.action . |
init |
vector of initial values of the iteration. Default initial value is zero for all variables. |
control |
a list of parameters for controlling the fitting process.
See the documentation for coxph.control for details.
|
... |
other arguments will be passed to coxph.control . |
an object of class coxph
. See coxph.object
and coxph.detail
for details.
y |
an object of class Surv containing the transformed times.
|
Stare J., Henderson R., Pohar M. "An individual measure for relative survival." Journal of the Royal Statistical Society: Series C , 54 115-126, 2005.
Pohar M., Stare J. "Relative survival analysis in R." Computer Methods and Programs in Biomedicine, 81: 272-278, 2006.
rsmul
, invtime
, rsadd
, survexp
.
data(slopop) data(rdata) rstrans(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365, sex=sex,year=year),ratetable=slopop,data=rdata)