Past exam of the mathematical and Computational Biology course of the University of Cambridge / 2025 / Paper 2 / 3 / c / Solution

Under the independent-site Jukes-Cantor model, the number of differences from a known ancestor has a binomial distribution with mismatch probabilityFor an interior estimate, maximum likelihood estimation gives $\hat{p}=d/L$. Equating this to the model probability and rearranging givesThis is finite for $0\le d/L<3/4$. At $d/L=3/4$, the estimate tends to infinity: the endpoint distribution has reached its limiting saturation and no finite rate can reproduce it exactly. For $d/L>3/4$, the displayed logarithm is not defined as a real number; within this model the constrained likelihood is maximized only as $\mu \to \infty$. For $d=0$ the boundary estimate is zero. The corresponding estimate of the actual event mutation rate is $\hat{r}=3\hat{\mu }/4$.