BUG: fix an interaction term in the calculation of expected numbers of transitions along edges, and add chatty tests, also the analytical expected log likelihood now works when leaf states are held constant, as tested by comparison to Monte Carlo Rao-Teh estimates