Fitting mixture of binomials

I simplified your code to fix the fraction parameter at 0.5 instead of letting it vary across observations. When you let it vary, what variation is causing what (or caused by what) becomes ambiguous. But the point about incompatibility between your DGP and the model remains.