Multinomial GLMM Maximum likelihood

This is an internal function called by mnGLMM.

Usage

mnGLMM.ml(par, par.fixed, Y, X, fitted.values = FALSE, REML = FALSE)