Binomial Approx to Normal

Purpose
Binomial approximates to Normal if np > 5 and nq > 5

> par(mfrow = c(2, 1))
> n <- 100
> p <- 0.5
> q <- 1 - p
> plot(dbinom(1:n, n, p), type = "l", lwd = 2, col = "sienna",
+     ylab = "", main = "Bin(n,p)")
> normal.mean <- n * p
> normal.std <- sqrt(n * p * q)
> plot(dnorm(1:100, normal.mean, normal.std), main = "Normal(np, npq)",
+     type = "l", ylab = "", lwd = 2, col = "sienna")