search for: pfit2

...pfit1 <- coxph(Surv(time, status > 0) ~ trt + log(bili) + log(protime) + age + platelet, data=pbc) #this works as expected survexp(~ trt, ratetable = pfit1, data = pbc) #now I fit a similar model, simply adding the 'sex' term #note sex is now the only factor in the model pfit2 <- coxph(Surv(time, status > 0) ~ trt + log(bili) + log(protime) + age + platelet + sex, data = pbc) #however, I now get an error with survexp survexp(~ trt, ratetable = pfit2, data = pbc) Error in ratetable(trt = trt, bili = bili, protime = protime, age = age, : object 'sexf&...

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]