## James-Stein shrinkage
shrink <- function(x, n) {
p <- length(x)
lam <- (p-3)/(n*p-p)
xbar <- mean(x)
xbar + (x-xbar)/(1+lam)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.