Nothing
#' Comparison Kernel Density Estimate Plot
#'
#' Produces side-by-side kernel density estimate plots.
#'
#'
#' @param x a \code{fit.models} object.
#' @param fun a function to extract the desired quantity from \code{x}.
#' @param \dots additional arguments are passed to
#' \code{\link[lattice]{xyplot}}.
#' @return the \code{trellis} object is invisibly returned.
#' @keywords hplot
#' @importFrom lattice densityplot panel.densityplot panel.abline strip.default
#' @export
sideBySideKernelDensityPlot <- function(x, fun, ...)
{
n.models <- length(x)
mod.names <- names(x)
y <- lapply(x, fun)
n.y <- sapply(y, length)
panel.special <- function(x, y, ...)
{
panel.densityplot(x, ...)
panel.abline(v = 0, lty = 2)
invisible()
}
mod <- factor(rep(mod.names, n.y), levels = mod.names)
tdf <- data.frame(y = unlist(y), mod = mod)
p <- densityplot(~ y | mod,
data = tdf,
n = 256,
bw = "SJ",
plot.points = FALSE,
ref = TRUE,
panel = panel.special,
strip = function(...) strip.default(..., style = 1),
layout = c(n.models, 1, 1),
...)
print(p)
invisible(p)
}
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.