compute the hyperparameters of an Normal-Inverse-Gamma distribution using an empirical Bayes approach