Displaying 1 result from an estimated 1 matches for "smessobs".
Did you mean:
messobs
2012 Feb 07
2
Actual vs. predicted survival times
...essobs <- 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 Time",xlab="Predicted Survival Time",main="Survival Times",xlim=c(0,3500),ylim=c(0,3500))
points(smessobs$explp2,messobs$tsecond,pc...