hyperEM function

Estimate hyperparameters using an EM algorithm