Nothing
err <- function(n, maxit) {
if (n == 0)
msg <- ""
if (n < 0) {
msg <-
paste(
"convergence for ",-n, "th lambda value not reached after maxit=",
maxit, " iterations; solutions for larger lambdas returned",
sep = ""
)
n <- -1
msg <- paste("From KERE fortran code -", msg)
}
list(n = n, msg = msg)
}
cvcompute <- function(mat, foldid, nlams) {
###Computes the weighted mean and SD within folds, and
# hence
# the se of the mean
nfolds <- max(foldid)
outmat <- matrix(NA, nfolds, ncol(mat))
good <- matrix(0, nfolds, ncol(mat))
mat[is.infinite(mat)] <- NA
for (i in seq(nfolds)) {
mati <- mat[foldid == i,]
outmat[i,] <- apply(mati, 2, mean, na.rm = TRUE)
good[i, seq(nlams[i])] <- 1
}
N <- apply(good, 2, sum)
list(cvraw = outmat, N = N)
}
ercls <- function(r, omega) {
abs(omega - (r < 0)) * r ^ 2
}
error.bars <- function(x, upper, lower, width = 0.02,
...) {
xlim <- range(x)
barw <- diff(xlim) * width
segments(x, upper, x, lower, ...)
segments(x - barw, upper, x + barw, upper, ...)
segments(x - barw, lower, x + barw, lower, ...)
range(upper, lower)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.