powered by
Maximizes log(likelihood) wrt phi.
McCullagh_conditional_symmetry_maximize_phi(n)
phi matrix
matrix of observed counts