Fitting mixture of binomials

Lets say I want to write a mixture to model what I am simulating (which sounds like maybe the issue is I wrote the model wrong in pymc). Could you help me with that?