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]]