A function to calculate posterior quantities of bivariate normals. See pages 79-86.
normal.posterior.summary(reps)
reps