Statistical tests for NormalGamma and NormalInverseGamma

This pull request adds `mean` and `var` for `NormalGamma` and `NormalInverseGamma`. In addition, it uses those moments to verify that `rand` is working correctly for those distributions.
You might want to keep the indentation consistent here.