fayHerriotEB function

Obtain empirical Bayesian estimates for group j