search for: positive1

Hi there, I have a rather large data set and perform the following cox model: test1 <- list(tstart,tstop,death1,chemo1,radio1,horm1) out1<-coxph( Surv(tstart,tstop, death1) ~ chemo1+chemo1:log(tstop+1)+horm1+horm1:log(tstop+1)+age1+grade1+grade1:log(tstop+1)+positive1+positive1:log(tstop+1)+size1+size1:log(tstop+1), test1) out1 Up to here everything works fine (with each covariate having a length of 289205) Now I want to see a specific profile of the above model so I ask for: x11() profilbig2=(survfit(out1,newdata=data.frame(chemo1=rep(0,length(chemo1)),horm1...