Basic.y function

Helper function for computing the posterior mean, posterior variance