[R] How to obtain a 95% envelope for the estimated cumulative risk function from bootstrap samples
wanr@ucalgary.ca
wanr at ucalgary.ca
Fri Sep 10 21:38:13 CEST 2004
Hi all,
I am trying to replicate the results from a paper. The problems are in the
setting of survival analysis.
The hypothetical data is displayed as follows.
ID time status
1 81 0
2 42 1
3 37 1
4 54 0
5 35 0
6 38 1
7 29 0
8 40 0
Question 1: How to obtain a 95% envelope for the estimated cumulative risk
function from bootstrap samples? I guess the output of this step consists of
the envelope and the estimated cumulative risk function.
Question2: bootstrap process will be repeated for n times; then the average
cumulative risk function was estimated as the median of the n empirical
cumulative risk functions. My question is how to plot this so called Average
Risk Function vs. time easily instead of using lines() and points() to draw
the above step function?
Thanks in advance.
Rui
More information about the R-help
mailing list