[R] competing risk analysis
Ravi Varadhan
rvaradhan at jhmi.edu
Fri Apr 30 16:26:42 CEST 2010
The object `xx' is a list with 3 components time, est, and var. You need
the `est' component of this list. So, you simply do the following:
bbb <- xx$est
There is no function called `cuminc.est' in the "cmprsk" package.
Hope this helps,
Ravi.
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
Behalf Of paaventhan jeyaganth
Sent: Thursday, April 29, 2010 3:42 PM
Cc: r
Subject: [R] competing risk analysis
Deart All,
I am doing competing risk analysis, I am tring store the estimate value in
bbb
but it gives error, it's anybody knows how to store in estimated value.
Thanks
Paaveen
xx <- cuminc(wkdt$fu.year, wkdt$status)
bbb <- cuminc.est(xx,F)
Error in cuminc.est: Argument number 2 not matched: cuminc.est(xx, F)
cuminc.est(xx)
[1] "Competing-OtrMrt: 38/171 patients"
estimated lower(95%) upper(95%)
0 0.00 0.00 0.00
1 5.65 2.05 9.25
2 6.39 2.53 10.25
3 7.23 3.07 11.39
4 12.04 6.32 17.75
5 18.53 11.23 25.83
6 22.39 14.21 30.58
7 26.87 17.69 36.05
8 33.08 22.27 43.89
9 41.92 29.47 54.38
10 41.92 29.47 54.38
11 41.92 29.47 54.38
12 45.61 31.90 59.31
13 49.58 34.70 64.45
[1] "Competing-CssMrt: 3/171 patients"
estimated lower(95%) upper(95%)
0 0.00 0.00 0.00
1 0.00 0.00 0.00
2 1.48 -0.56 3.53
3 1.48 -0.56 3.53
4 1.48 -0.56 3.53
5 1.48 -0.56 3.53
6 2.80 -0.48 6.07
7 2.80 -0.48 6.07
8 2.80 -0.48 6.07
9 2.80 -0.48 6.07
10 2.80 -0.48 6.07
11 2.80 -0.48 6.07
12 2.80 -0.48 6.07
13 2.80 -0.48 6.07
_________________________________________________________________
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list