sigmashrink function

Empirical Bayes shrinkage estimate of sigma^2