Displaying 1 result from an estimated 1 matches for "statrep".
Did you mean:
statret
2006 Aug 17
1
putting the mark for censored time on 1-KM curve or competing risk curve
...urvfit". 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",...