[R] R error using Survr function with gcmrec

Caroline King ccostell at ucla.edu
Tue Jun 17 16:56:26 CEST 2008


Would someone be able to help with this question?  I'm using the  
Gcmrec, Survrec, and Design packages to do a power analysis on  
simulated data.  I'm receiving an error after using the Survr function  
that all data must have a censoring time even after using the gcmrec  
function:  newdata<-addCenTime(olddata).  My program is below.  I'd  
greatly appreciate any help!

id<-c(seq(1,288,by=1),seq(1,79,by=1),seq(1,11,by=1))
x<-c(rep(0,5),rep(1,6),rep(0,45),rep(1,23),rep(0,124),rep(1,85),
  +rep(0,4),rep(1,1),rep(0,1),rep(1,5),rep(0,31),rep(1,14),rep(0,5),
  +rep(1,18),rep(0,8),rep(1,3))
myrates<-((1-x)*0.0639 + (x)*0.0320)
y<-c(rexp(378,rate=myrates))
cen<-c(rexp(378,0.0385))
time<-pmin(y,cen)
event<-as.numeric(y<=cen)
x2<-(x-1)*(-1)
bvdata<-data.frame(id,event,time,x2)
bvdata2<-addCenTime(bvdata)
fit<-cph(Survr(id,event,time)~x2,data=bvdata2)



More information about the R-help mailing list