search for: messcen

...are predict.coxph and predict.cph options available to me. My first attempt was as follows: # fit Cox model # coxfita = coxph(Surv(tsecond,seccens)~stroke(smess)+rels(smess)+asleep(smess)+eeg1(smess)+eeg2(smess)+ct1(smess)+ct2(smess)+treat(smess),data=smess) # Find censored and observed groups # messcens <- subset(smess,seccens==1) messobs <- subset(smess,seccens==0) # Obtain predicted survival times # explp <- exp(predict(coxfita,type="lp")) explp2 <- mean(ssmess$tsecond,na.rm=TRUE)*explp smess2 <- data.frame(ssmess,explp2) # Find censored and observed groups # smesscen...