polyaHessian function

Obtain Hessian of the marginal Dirichlet-multinomial likelihood