lattice.normal <-
function(x, mu, SE, ...) {
    panel.histogram(x, ...)
    panel.mathdensity(dmath = dnorm, col = "black", lwd=2,
    args = list(mean = mu, sd = SE))
}