[R] ploting survival curves

[raymond chiruka]

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


______________________________________________
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.