# Bootstrap method for calculating CVSE for Cox models
se.ncvsurv <-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.ncvsurv(y[ind,], eta[ind,])
}
apply(cve, 2, sd)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.