i'm trying to plot survival curves using the following code
n=20
n1=n/2
n2=n/4
a11=4;a12=4 ;a21=4 ;a22=4
t1<-array(1,c(n1))
t2<-array(2,c(n1))
treatgrp=matrix(c(t1,t2))
st11<-array(1,c(n2))
st12<-array(2,c(n2))
st21<-array(1,c(n2))
st22<-array(2,c(n2))
strata=matrix(c(st11,st12,st21,st22))
time11=array(rweibull(n2,a11,1))
time12=array(rweibull(n2,a11,1))
time21=array(rweibull(n2,a21,1))
time22=array(rweibull(n2,a22,1))
time=matrix(c(time11, time12, time21, time22))
censorTime=runif(n,0,3)
m=cbind(treatgrp,strata,time,censorTime)
colnames(m)=c(“treat”,”strata”,”time”,”censorTime”)
act.surv.time<-pmin(m[,"time"],m[,"censorTime"])
m<-cbind(m,act.surv.time)
m<-cbind(m,0)
m[m[,3]<m[,4],6]<-1
colnames(m)[6]<-"censoring"
m
b=data.frame(m)
c1= survdiff(Surv(act.surv.time)~treatgrp ,data=b)
c2= survdiff(Surv(act.surv.time,censoring)~treatgrp ,data=b)
c1
c2
plot(c2,lty=2:3,fun=”cumhaz”,xlab=”time”,ylab=”cumulative hazard”)
but keep getting the following error
Error in plot.window(xlim, ylim, log, asp, ...) :
need finite 'xlim' values
In addition: Warning messages:
1: no non-missing arguments to min; returning Inf
2: no non-missing arguments to max; returning -Inf
3: no non-missing arguments to min; returning Inf
4: no non-missing arguments to max; returning -Inf
5: "fun" is not a graphical parameter in: plot.window(xlim, ylim,
log, asp, ..
what is the problem i have looked at the examples given.
__________________________________________________
[[alternative HTML version deleted]]