powered by
gradient of negative log-likelihood function for poisson endpoints
grad_poisson_likelihood(theta, n_exp, n_ref, n_pla, sum_exp, sum_ref, sum_pla)