polyaGradient function

Obtain gradient of the marginal Dirichlet-multinomial likelihood