ENH: added mean method for NormalInverseGamma

I suggest that you directly compute the mean using the formula, instead of constructing an intermediate InverseGamma instance.