powered by
This and other similar internal functions calculate the Hessian for trivariate binary models.