## Run radon model
## source("http://web.as.uky.edu/statistics/users/pbreheny/701/S13/notes/3-19.R")

## Finite population standard deviations
s <- apply(a, 1, function(x) sqrt(mean((x-mean(x))^2)))
psm(s)
psm(sigma)