Marginal maximum likelihood estimation is commonly used to estimate logistic-normal models. In this approach, the contribution of random effects to the likelihood is represented as an intractable integral over their distribution. Thus, numerical methods such as Gauss