Description Usage Arguments Details Value See Also Examples
Estimates the crude probability of death due to disease and due to population reasons
1 2 
formula 
a formula object, with the response as a NOTE: The time must be in days, and the same
is required for the ratetable variables (the variables used in the population tables), for example age and year (year must be
given in the 
data 
a data.frame in which to interpret the variables named in the

ratetable 
a table of event rates, organized as a 
na.action 
a missingdata filter function, applied to the model.frame,
after any subset argument has been used. Default is

tau 
the maximum followup time of interest, all times larger than 
conf.int 
the level for a twosided confidence interval on the survival curve(s). Default is 0.95. 
precision 
the level of precision used in the numerical integration of variance. Default is 1, which means that daily intervals are taken, the value may be decreased to get a higher precision or increased to achieve a faster calculation. The calculation intervals always include at least all times of event and censoring as border points. 
add.times 
specific times at which the value of estimator and its variance should be evaluated. Default is all the event and censoring times. 
NOTE: All times used in the formula argument must be specified in days. This is true for the followup time as well as for
any variables needed ratetable
object, like age
and year
.
Note that numerical integration is required to calculate the variance estimator. The integration precision is set with argument precision
, which defaults to daily intervals, a default that should give enough precision for any practical purpose.
The area under the curve is calculated on the interval [0,tau
].
Function summary
may be used to get the output at specific points in time.
An object of class cmp.rel
.
Objects of this class have methods for the functions print
and plot
. The summary
function can be used for printing output at required time points.
An object of class cmp.rel
is composed of several lists, each pertaining the cumulative hazard function for one risk and one strata. Each of the lists contains the following objects:
time 
the timepoints at which the curves are estimated 
est 
the estimate 
var 
the variance of the estimate 
lower 
the lower limit of the confidence interval 
upper 
the upper limit of the confidence interval 
area 
the area under the curve calculated on the interval [0, 
index 
indicator of event and censoring times among all the times in the output. The times added via paramater 
add.times 
the times added via parameter 
1 2 3 4 5 6 7 8 9  data(slopop)
data(rdata)
#calculate the crude probability of death
#note that the variable year is in the date format  the number of days since 01.01.1960 and that
#age must be multiplied by 365 in order to be expressed in days.
fit < cmp.rel(Surv(time,cens)~sex+ratetable(age=age*365,sex=sex,
year=year),ratetable=slopop,data=rdata,tau=3652.41)
fit
plot(fit,col=c(1,1,2,2),xscale=365.241,xlab="Time (years)")

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.