powered by
Compute sums too use in maximizing log(likelihood)
McCullagh_compute_s_plus(n)
list of s_i_plus and s_plus_i
matrix of observed counts