Hierarchical model varying indexing of distributions

Very small nitpick: Those are statistically the same but need not be numerically the same (allclose). It just happens that location/scale doesn’t influence numpy RNG, but they could in other algorithms/ loc-scale families (e.g., when the number of random bits consumed depends on the parameter values).

1 Like