search for: explp

...ita = 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) smessobs <- subset(smess2,seccens==0) # Produce plot # plot(smesscens$explp2...

...kp]         tframe <- mframe[rep(1,length(var.values)),]         tframe[,var.num] <- var.values         xmat <- model.matrix(cfit,data=tframe)[,-1]         tlps <- as.vector(xmat%*%cfit$coef)         names(tlps) <- var.values         obs.off <- tlps[as.character(obs.var)]          explp.off <- exp(outer(lps - obs.off ,tlps,"+")) bfit <- survfit.coxph(cfit, se.fit = F)         fits <- outer(bfit$surv,explp.off,function(x,y) x^y)         avg.fits <-             apply(sweep(fits,2,weights,"*"),c(1,3),sum)/sum(weights)         dimnames(avg.fits) <- l...

...tframe <- mframe[rep(1,length(var.values)),] tframe[,var.num] <- var.values xmat <- model.matrix(cfit,data=tframe)[,-1] tlps <- as.vector(xmat%*%cfit$coef) names(tlps) <- var.values obs.off <- tlps[as.character(obs.var)] explp.off <- exp(outer(lps - obs.off ,tlps,"+")) bfit <- survfit.coxph(cfit, se.fit = F) fits <- outer(bfit$surv,explp.off,function(x,y) x^y) avg.fits <- apply(sweep(fits,2,weights,"*"),c(1,3),sum)/sum(weights) dimnames(avg.fits) <...