Displaying 1 result from an estimated 1 matches for "ssmess".
Did you mean:
smess
2012 Feb 07
2
Actual vs. predicted survival times
...)+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)
smessobs <- subset(smess2,seccens==0)
# Produce plot #
plot(smesscens$explp2,messcens$tsecond,pch=4,col="blue",ylab="Actual Survival T...