plot.cv.gcdclust <-
function(x, sign.lambda = 1, ...) {
cvobj <- x
xlab <- "log(Lambda)"
if (sign.lambda < 0)
xlab <- paste("-", xlab, sep = "")
plot.args <- list(x = sign.lambda * log(cvobj$lambda), y = cvobj$cvm,
ylim = range(cvobj$cvupper, cvobj$cvlo), xlab = xlab,
ylab = cvobj$name, type = "n")
new.args <- list(...)
if (length(new.args))
plot.args[names(new.args)] <- new.args
do.call("plot", plot.args)
error.bars(sign.lambda * log(cvobj$lambda), cvobj$cvupper,
cvobj$cvlo, width = 0.02, col = "darkgrey")
points(sign.lambda * log(cvobj$lambda), cvobj$cvm, pch = 20,
col = "red")
axis(side = 3, at = sign.lambda * log(cvobj$lambda), labels = paste(cvobj$nz),
tick = FALSE, line = 0)
abline(v = sign.lambda * log(cvobj$lambda.min), lty = 3)
abline(v = sign.lambda * log(cvobj$lambda.1se), lty = 3)
invisible()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.