search for: smesscens

...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 Time",xlab="Predicted Survival Time",main="Survival Times",xlim=c(0,3500),ylim=c(0,3500))...