Update EM algorithm after checking

I tried the EM algorithm with the true values of the CCPs and the true pi_s and got the same pi_s (0.249265 instead of 0.25), so it looks like it works. The issue with getting pi_s close to 1 or zero is something that happens with a bad CCP estimate. We may have to watch that.