Fitting mixture of binomials

That would be as in DGP#1 in my code here.

When mix_frac is 0.5, the model behaves as seen here.