¹The equations shown here for fast augmentation of the saturated posterior mean and variances were derived by Professor John Haslett.