search for: seccens

...at all. Given that, I was concerned about proceeding but I would like to have access to a plot to make a decision on its usefulness. I appreciate that there 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,...

...line hazard function and not the baseline survival estimate. Is there a way to obtain the baseline survival estimate or do I have to use the formula which does something like S(t) = exp[- the integral from 0 to t of h(u) du]? Thank you for your assistance, Laura fit1 <- coxph(Surv(tsecond/365,seccens)~stroke(smess1)+othnd(smess1)+relat(smess1)+asleep(smess1)+abeeg1(smess1)+treat(smess1),data=smess1) basehaz(fit1) where stroke is a function which creates a binary variable from the dataset smess1 etc. [[alternative HTML version deleted]]