Very different sampling with ~beta_binomial than with += beta_binomial_lpmf

Yes this is know bug. See Rstan works but cmdstan does not for a beta binomial model

This was fixed for the 2.22 version of Stan. As far as I understand it was introduced somewhere between 2.19 and 2.21.

4 Likes