Lineplot of densities

plotDensities.R

plotDensities <- function(s, col=NULL, xlim=NULL, ylim=NULL, na.rm=TRUE, lty=rep(1, length(s)), ...) {
    if (!is.list(s)) s <- list(s)
    if (is.null(col)) col <- rainbow(length(s))
    sD <- lapply(s, density, na.rm=na.rm)
    if (is.null(xlim)) xlim <- range(sapply(sD, function(x) range(x$x)))
    if (is.null(ylim))ylim <- range(sapply(sD, function(x) range(x$y)))
    plot(0, type="n", xlim=xlim, ylim=ylim, ...)
    for (i in 1:length(sD)) lines(sD[[i]], col=col[i], lty=lty[i], ...)
}