rs.surv {relsurv}R Documentation

Compute a Relative Survival Curve

Description

Computes an estimate of the relative survival curve using the Kaplan-Meier method for the observed and the Hakulinen or Ederer II method for the expected survival.

Usage

rs.surv(formula, data,ratetable=survexp.us,fin.date,method="hakulinen",...) 

Arguments

formula a formula object, with the response as a Surv object on the left of a ~ operator, and, if desired, terms separated by the + operator on the right. If the variables are not organized and named in the same way as in the population tables, a ratetable term must be added to match each subject to his/her expected cohort. The time must be in days. For a single survival curve the ~ 1 part of the formula is not required.
data a data.frame in which to interpret the variables named in the formula.
ratetable a table of event rates, organized as a ratetable object, such as survexp.us.
fin.date The date of the study ending, used for calculating the potential follow-up times in the Hakulinen method. If missing, it is calculated as max(year+time). Only to be specified when using the method=hakulinen.
method The method to be chosen for the estimation of the population survival. hakulinen for the Hakulinen method, conditional for the Ederer II method.
... other arguments will be passed to the survfit function that calculates the observed survival.

Details

The potential censoring times needed for the calculation of the expected survival by the Hakulinen method are calculated automatically. The times of censoring are left as they are, the times of events are replaced with fin.date - year.

Value

a survfit object; see the help on survfit.object for details. The survfit methods are used for print, plot, lines, and points.

See Also

survfit, survexp

Examples

data(slopop)
data(rdata)
rs.surv(Surv(time,cens)~sex+ratetable(age=age*365,sex=sex,
      year=year),ratetable=slopop,data=rdata)

[Package relsurv version 1.3 Index]