residuals.rsadd: Calculate Residuals for a "rsadd" Fit

Description

Calculates partial residuals for an additive relative survival model.

Usage

# S3 method for rsadd
residuals(object, type = "schoenfeld", ...)

Value

A list of the following values is returned:

res: a matrix containing the residuals for each variable.
varr: the variance for each residual
varr1: the sum of varr.
kvarr: the derivative of each residual, to be used in rs.zph function.
kvarr1: the sum of kvarr.

Arguments

object: an object inheriting from class rsadd, representing a fitted additive relative survival model. Typically this is the output from the rsadd function.
type: character string indicating the type of residual desired. Currently only Schoenfeld residuals are implemented.
...: other arguments.

References

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.

Goodness of fit: Stare J.,Pohar Perme M., Henderson R. (2005) "Goodness of fit of relative survival models." Statistics in Medicine, 24: 3911--3925.

Examples

Run this code


data(slopop)
data(rdata)
fit <- rsadd(Surv(time,cens)~sex,rmap=list(age=age*365.241),
       ratetable=slopop,data=rdata,int=5)
sresid <- residuals.rsadd(fit)