Open jar398 opened 4 years ago
In normal.py, the comment says rho ~ Gamma(nu/2, s/2) but the code says rho = rng.gamma(nu/2., scale=2./s)
rho ~ Gamma(nu/2, s/2)
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.
Gamma
rng.gamma
In normal.py, the comment says
rho ~ Gamma(nu/2, s/2)
but the code saysrho = 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
andrng.gamma
do different things.