Displaying 1 result from an estimated 1 matches for "messobs".
Did you mean:
smessobs
2012 Feb 07
2
Actual vs. predicted survival times
...s 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 #
smesscens <- subset(smess2,seccens==1)
smesso...