Nothing
denoverplot1 <- function(..., ci = NULL, col=NULL, lty=1, xlim=NULL, ylim=NULL, xlab = "", ylab = "Density", main = NULL, style=c("gray", "plain"), gpar=NULL){
dat <- list(...)
style <- match.arg(style)
if (length(dat)==1 && is.list(dat[[1]])) dat <- dat[[1]]
n <- length(dat)
if (n==1 && is.list(dat[[1]])) {
dat <- dat[[1]]
n <- length(dat)
}
if (is.null(col)){
col <- mcmcplotsPalette(n)
}
denout <- lapply(dat, density, bw="SJ")
xx <- sapply(denout, function(den) den$x)
yy <- sapply(denout, function(den) den$y)
if (style=="plain")
do.call("matplot", c(list(x=xx, y=yy, col=col, lty=lty, xlim=xlim, ylim=ylim, xlab=xlab, ylab=ylab, main=main, type="l"), gpar))
if (style=="gray"){
do.call("matplot", c(list(x=xx, y=yy, type="n", bty="n", xaxt="n", yaxt="n", xlim=xlim, ylim=ylim, xlab=xlab, ylab=ylab, main=main)))
.graypr()
do.call("matlines", c(list(x=xx, y=yy, col=col, lty=lty), gpar))
}
if (!is.null(ci)){
lb <- sapply(dat, quantile, (1 - ci)/2)
ub <- sapply(dat, quantile, ci + (1 - ci)/2)
do.call("abline", c(list(v=lb, col=col, lty=lty), gpar))
do.call("abline", c(list(v=ub, col=col, lty=lty), gpar))
}
}
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.