powered by
Helper function for computing the posterior mean, posterior variance
Basic.y(y, tau, k, data.var)