Graficar histograma de aleatorios y densidad.

set.seed(999)
exponenciales <- rexp(10000)
hist(exponenciales, col = "yellowgreen", freq = FALSE)
curve(dexp(x, 1), 0, 6, n = 101, add = TRUE, lwd = 3, col = "red")

plot of chunk unnamed-chunk-1

Generar números aleatorios exponenciales a partir de aleatorios uniformes.

nsim = 10000  # número de aleatorios uniformes
u = runif(nsim)
y = -log(u)
par(mfrow = c(1, 2))
hist(y, col = "tomato", freq = FALSE, xlim = c(0, 4), breaks = 19, main = "Aleatorios exponenciales", 
    sub = "A partir de aleatorios uniformes", ylab = "Densidad")
hist(rexp(nsim), col = "blue", freq = FALSE, xlim = c(0, 4), main = "Aleatorios generados por R", 
    ylab = "Densidad")

plot of chunk unnamed-chunk-2

date()

## [1] "Mon Oct  7 20:15:04 2013"