[R] R error using Survr function with gcmrec
Caroline King
ccostell at ucla.edu
Tue Jun 17 18:26:59 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