Displaying 1 result from an estimated 1 matches for "chaz2".
Did you mean:
chaz
2013 Feb 12
0
error message from predict.coxph
...<- (chaz * newx[indx2, ]) - xbar[indx2, ] # *** SHOULD BE JUST xbar
# se[indx2] <- sqrt(varh + rowSums((dt %*% object$var) * dt)) * newrisk[indx2]
# }
#}
#else {
# j2 <- approx(afit$time, 1:afit.n, newy[indx2, 2], method = "constant", f = 0, yleft = 0, yright = afit.n)$y # chaz2 <- approx(-afit$time, afit$cumhaz, -newy[indx2, 2], method = "constant", rule = 2, f = 0)$y # chaz2 <- c(0, afit$cumhaz)[j2 + 1] # pred[indx2] <- (chaz2 - chaz) * newrisk[indx2] # if (se.fit) {
# varh2 <- c(0, cumsum(afit$varhaz))[j1 + 1]
# xbar2 <- rbind(0, afit$x...