search for: time30

How do you get the design matrix R used when calculating ANOVA? -- View this message in context: http://www.nabble.com/how-to-get-design-matrix--tp23464638p23464638.html Sent from the R help mailing list archive at Nabble.com.

...I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong? My code is reproduced below, my figure is attached to this EMail message. John > #Create simple survival object > GVHDdata<-list(Time=GVHD$Time,Time30=(GVHD$Time)/30, + Age=GVHD$Age,Drug=GVHD$Drug,Died=GVHD$Died, + AgeGrp=cut(GVHD$Age,breaks=c(0,15,25,45))) > > summary(GVHD$Drug) MTX MXT+CSP 32 32 > > > > fit0<-coxph(Surv(Time30,Died)~Drug,data=GVHDdata) > summary(fit0) Call: coxph(formula = Surv(Time30...

...should give two lines (one for each treatment) gives only one line and a warning message. I would appreciate help getting two lines, and an explanation of the warning message. My problem may the that I have very few events in one of my strata, but I don't know. Thanks, John fit0<-coxph(Surv(Time30,Died)~strata(Rx)+Age,data=GVHDdata) plot(survfit(fit0),fun="cloglog") WARNING: Warning in xy.coords (x, y, xlabel, ylabel, log) : 2 x values <=0 omitted from logarithmic plot > print(survfit(fit0),fun="cloglog") Call: survfit.coxph(object = fit0) n events med...