Finished up the EM algorithm today. Key was getting the right function to maximize. Turns out wikipedia has a very nice write up of this very example, but in our notation:
Eplog(pf1)+(1−Ep)log((1−p)f2)
Where Ep comes from the expectation step.
Jamie has joined us and has set up a github repository for the discussion group. Find the successful abstract algorithm there.