Estimates the distribution of true effects (or any latent
parameter) from estimates that each have a known, possibly different,
standard error. The main function, em_smooth(), implements the smoothed
EM (EMS) algorithm of Silverman, Jones, Wilson and Nychka (1990)