Dear All, here is the example for cumulative incidence analysis with cmprsk package:
set.seed(2) ss <- rexp(100) gg <- factor(sample(1:3,100,replace=TRUE),1:3,c('a','b','c')) cc <- sample(0:2,100,replace=TRUE) strt <- sample(1:2,100,replace=TRUE) print(xx <- cuminc(ss,cc,gg,strt)) plot(xx,lty=1,color=1:6) When I perform this example, I have 6 curves. a 1 b 1 c 1 a 2 b 2 c 2 I would like to plot only 3 first curves of risk cc=1 not all of 6 curves. a 1 b 1 c 1 How I can do this with R ? Many thanks Jan [[alternative HTML version deleted]] ______________________________________________ R-help@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.