R/se-grpsurv.R

# Bootstrap method for calculating CVSE for Cox models
se.grpsurv <-function (y, eta, B = 100) {
  cve <- matrix(NA, B, ncol(eta))
  for (b in 1:B) {
    ind <- sample(1:nrow(eta), replace=TRUE)
    cve[b,] <- loss.grpsurv(y[ind,], eta[ind,])
  }
  apply(cve, 2, sd)
}
pbreheny/grpreg documentation built on Nov. 29, 2023, 1:16 a.m.