Dear all,

It's fairly straightforward to plot cumulative histograms using the hist()
function. You do something like:

h <- hist(rnorm(100), plot=FALSE)
h$counts<- cumsum(h$counts)
plot(h)

However, I have failed to find any example where this is done using the
lattice histogram() function. I realize I need to slightly alter the panel
function panel.histogram. Specifially I would like to add the following line
in red, just like I did above:

function (x, breaks, equal.widths = TRUE, type = "density", nint =
round(log2(length(x)) +
    1), alpha = plot.polygon$alpha, col = plot.polygon$col, border =
plot.polygon$border,
    lty = plot.polygon$lty, lwd = plot.polygon$lwd, ...)
{
    plot.polygon <- trellis.par.get("plot.polygon")
    xscale <- current.panel.limits()$xlim
    panel.lines(x = xscale[1] + diff(xscale) * c(0.05, 0.95),
        y = c(0, 0), col = border, lty = lty, lwd = lwd, alpha = alpha)
    if (length(x) > 0) {
        if (is.null(breaks)) {
            breaks <- if (is.factor(x))
                seq_len(1 + nlevels(x)) - 0.5
            else if (equal.widths)
                do.breaks(range(x, finite = TRUE), nint)
            else quantile(x, 0:nint/nint, na.rm = TRUE)
        }
        h <- hist.constructor(x, breaks = breaks, ...)

        h$counts<- cumsum(h$counts)

        y <- if (type == "count")
            h$counts
        else if (type == "percent")
            100 * h$counts/length(x)
        else h$intensities
        breaks <- h$breaks
        nb <- length(breaks)
        if (length(y) != nb - 1)
            warning("problem with 'hist' computations")
        if (nb > 1) {
            panel.rect(x = breaks[-nb], y = 0, height = y, width =
diff(breaks),
                col = col, alpha = alpha, border = border, lty = lty,
                lwd = lwd, just = c("left", "bottom"))
        }
    }
}
<environment: namespace:lattice>


My problem is I'm too unexperienced in handling these panel functions to
achieve this. Simply copying the panel function, appending the line, and
giving it another name obviously doesn't work (it won't find the function
hist.constructor). I would very much appreciate help on how this could be
done, or some other way to draw cumulative lattice histograms.

Thanks in advance,
Ola Caster

	[[alternative HTML version deleted]]

