Various changes:  added hrss calculation, changed discrete summation so that calculated hrss matches injected hrss, fixed normalizations of prior and determinant terms, calculate baye's factor