rsadd {relsurv} | R Documentation |
The function fits an additive model to the data. The methods implemented are the maximum likelihood method,
a glm model with a binomial
error and a glm model with a poisson
error.
rsadd(formula, data=parent.frame(), ratetable = survexp.us, int, na.action, method, init,centered,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, organized as a ratetable object, such as survexp.us .
|
int |
either a single value denoting the number of follow-up years or a vector
specifying the intervals (in years) in which the hazard is constant (the times that are
bigger than max(int) are censored. If missing, the intervals are set to be one year
long and include 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 . |
method |
glm.bin or glm.poi for a glm model and max.lik for the maximum likelihood model (default). |
init |
vector of initial values of the iteration. Default initial value is zero for all variables. |
centered |
if TRUE , all the variables are centered before fitting and the baseline excess hazard
is calculated accordingly. Default is FALSE . |
control |
a list of parameters for controlling the fitting process.
See the documentation for glm.control for details.
|
... |
other arguments will be passed to glm.control . |
An object of class rsadd
. In the case of method="glm.bin"
and method="glm.poi"
the class also
inherits from glm
which inherits from the class lm
.
Objects of this class have methods for the functions print
and summary
.
An object of class rsadd
is a list containing at least the following components:
data |
the original data, along with the variables defined in the rate table |
warnme |
the number of groups with the observed number of deaths is smaller than the expected. |
ratetable |
the ratetable used. |
int |
the maximum time (in years) used. All the events at and after this value are censored. |
method |
the fitting method that was used. |
linear.predictors |
the vector of linear predictors, one per subject. |
groups |
life tables for groups specified by the covariates. |
Dickman P.W., Sloggett A., Hills M., Hakulinen T. "Regression models for relative survival." Statistics in Medicine, 2004.
data(slopop) data(rdata) rsadd(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365,sex=sex, year=year),ratetable=slopop,data=rdata,int=5)