Inference based on mixed-effect models, including generalized linear mixed models with spatial correlations, multivariate responses, and models with non-Gaussian random effects (e.g., Beta). Various approximations of likelihood or restricted likelihood are implemented, in particular Laplace approximation and h-likelihood (Lee and Nelder 2001