Displaying 1 result from an estimated 1 matches for "smess".
Did you mean:
mess
2012 Feb 07
2
Actual vs. predicted survival times
...hat, 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,type="lp&...