s/2 vs. 2/s

[jar398]

In normal.py, the comment says rho ~ Gamma(nu/2, s/2) but the code says rho = rng.gamma(nu/2., scale=2./s)

This is confusing. Either the comment is wrong, or else the presentation would be improved by explaining that Gamma and rng.gamma do different things.