pub fn mle(counts: &[u16], p: usize, q: usize, relerr: f64) -> f64