Risanje zamaknjenih normalnih porazdelitev

Funkcija za risanje ene normalne porazdelitve

plotNormal <- function(x0 = 0, y0 = 0, dx = 0, x = seq(-3, 3, 0.1), lwd = 3, 
    ...) {
    y <- dnorm(x)
    lines(x0 + dx + x, y0 + y, lwd = lwd)
    segments(x0 + dx + min(x), y0, x0 + dx + max(x), y0)
    segments(x0, y0, x0, y0 + (max(y) * 1.5), lty = 2)
    segments(x0 + dx, y0, x0 + dx, y0 + (max(y) * 1.5))
    if (dx != 0) 
        arrows(x0, y0 + (max(y) * 1.5), x0 + dx, y0 + (max(y) * 1.5), lwd = 4, 
            col = "darkgreen", code = 3, length = 0.15)
}

Narišimo končno sliko

plot(c(-4, 10), c(0, 4), type = "n", axes = FALSE, xlab = "", ylab = "")
plotNormal(0, 0)
plotNormal(2, 1, dx = -1)
plotNormal(4, 2, dx = +1.5)
segments(-1, -0.5, 5, 2.5)