rstrans {relsurv}R Documentation

Fit Cox Proportional Hazards Model in Transformed Time

Description

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.

Usage

rstrans(formula, data, ratetable, int,na.action,init,control,...)

Arguments

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.

Value

an object of class coxph. See coxph.object and coxph.detail for details.

y an object of class Surv containing the transformed times.

References

Stare J., Henderson R., Pohar M. "An individual measure for Relative Survival." Journal of the Royal Statistical Society: Series C , 54 115-126, 2005.

See Also

rsmul, invtime, rsadd, survexp.

Examples

data(slopop)
data(rdata)
rstrans(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365,
        sex=sex,year=year),ratetable=slopop,data=rdata)

[Package relsurv version 1.2 Index]