Hi All,
I have been working on a Competing Risks analysis using the mstate package
and when I run my analysis I was surprised to see the cumulative incident
function ending substantially before my data set did.
You can reproduce the behavior (kind of) using the data that was provided
with mstate:
library(mstate)
data(aidssi)
ci <- Cuminc(time=aidssi$time, status=aidssi$status)
and you can see that ci goes to 10.448
However, when you use the cmprsk package, the CIF goes pretty much to the
end of the data
library(cmprsk)
CI.overall <- cuminc(ftime =aidssi$time, fstatus=aidssi$status)
and you can see that CI.overall ends at 12
Using the data provided in the package, a difference of < 2 years is not
terrible but with my dataset mstate ends at 1.02 while cmprsk ends at 4 (my
max for event time is 4.9) and that is a substantial difference.
I believe I must be missing something here. Is there a way to get mstate to
report the entire CIF ?
--
- Don
Donald Catanzaro PhD
dgcatanz...@gmail.com
16144 Sigmond Lane
Lowell, AR 72745
479-721-2533
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
