## Demonstrate that t-distribution quantiles get closer to the ## standard normal quantile as the degrees of freedom increases. M = array(0, c(100,2)) q = 0.975 ## The quantile to analyze. ## Cycle through various degrees of freedom. for (d in 1:100) { ## The qth quantile of the t-distribution with d degrees of freedom. M[d,1] = qt(q, d) ## The qth quantile of the standard normal distribution ## (doesn't depend on d). M[d,2] = qnorm(q) }