\begin{Schunk}
\begin{Sinput}
par(mfrow=c(2,2))
N <- 2000
x <- seq(0.0000001, 2, length.out = 4000)
y1 <- rgamma(N, shape = 2, scale = 1/2)
h1 <- hist(log(y1), main = "shape = 2", breaks = 50, prob = TRUE,
# xlim = rockchalk::magRange(log(y1), c(1.25,1.25)),
xlim = c(-4, 2),
xlab = expression(log(delta[i])))
y2 <- rgamma(N, shape = 2.5, scale = 1/2.5)
hist(log(y2), main = "shape = 2.5", breaks = 50, prob = TRUE,
# xlim = rockchalk::magRange(log(y1), c(1.25,1.25)),
xlim = c(-4, 2),
xlab = expression(log(delta[i])))
y3 <- rgamma(N, shape = 5, scale = 1/5)
hist(log(y3), main="shape = 5", breaks = 50, prob = TRUE,
xlim = c(-4, 2),
xlab = expression(log(delta[i])))
y4 <- rgamma(x, shape = 10, scale = 1/10)
hist(log(y4), main="shape = 10", breaks = 50, prob = TRUE,
xlim = c(-4, 2),
xlab = expression(log(delta[i])))
par(mfrow=c(1,1))
\end{Sinput}
\end{Schunk}