updated bayesian, still need to look up E[X] of Gamma()