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