Results are quite different from the original in R. I can't seem to locate the issue. Parameter g should be much higher, around 0.28 and its HPD interval should be completely above zero to prove the point that we see a gender advantage when in reality there is none.
They are different because the code in R has a seed (set.seed(1913)) and that changes the number for g and d. If you rerun that model (changing G, S, D, A), the numbes will be almost the same.
You wrote
They are different because the code in R has a seed (
set.seed(1913)
) and that changes the number forg
andd
. If you rerun that model (changingG, S, D, A
), the numbes will be almost the same.