[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