Model comparison of generalized linear mixed models

Generalized linear mixed models (GLMMs) have been widely appreciated in biological and medical research. Maximum likelihood estimation has received a great deal of attention. Comparatively, not much has been done on model comparison or hypotheses testing. In this article, we propose a path sampling procedure to compute the observed-data log-likelihood function, so that the Bayesian information criterion (BIC) can be applied to model comparison or hypothesis testing. Advantages of the proposed path sampling procedure are discussed. Two medical data sets are analysed for providing illustrative examples of the proposed methodology.