[R] putting the mark for censored time on 1-KM curve or competing risk curve

From: Linda Lei <llei_at_bccrc.ca>
Date: Thu 17 Aug 2006 - 10:33:50 EST


Hi All,  

I'm trying to figure out the cumulative incidence curve in R in some limited time. I found in package "cmprsk", the command "plot.cuminc" can get this curve. But I noticed that there is no mark for the censored time there, comparing with the KM curve by "plot.survfit". Here are my codes (attached is the data):  


dat<-read.table("F://wendy/BMT data analysis/final data.txt",header=TRUE,sep="\t")  

library(cmprsk)  

library(survival)  

attach(dat)  

par(mfrow=c(2,1))  

curve<-cuminc(SURVEFP/365.25,STATREP,Ritux)  

plot(curve,xlab="years",curvlab=c("No rituximab", "Rituximab"),lty=1:2,col=c("red","blue"),mark.time=TRUE)  

a<-survfit(Surv(SURVEFP/365.25,STATREP)~Ritux,type="kaplan-meier")  

plot(a,main="OS for whole
group",conf.int=F,xlab="years",ylab="probability",lty=1:2,col=c("red","b lue"))


 

Could you help me with how to put the mark for censored time on cumulative incidence curve, or maybe how to get the 1-KM curve with the mark for censored time? This will help me a lot!  

Thank you very much!



R-help@stat.math.ethz.ch 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.
Received on Thu Aug 17 10:39:54 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 17 Aug 2006 - 12:22:21 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.