An alternative specification of generalized linear mixed models

Consider stratified data in which Y-i1, ... Y-inj, denote real-valued response variables corresponding to the observations from stratum i, i = 1, ... , m and suppose that Y-ij follows an exponential family distribution with canonical parameter of the form theta(ij) = x(ij)beta + gamma(i). In analyzing data of this type, the stratum-specific parameters are often modeled as random effects; a commonly-used approach is to assume that gamma(1), .... gamma(m) are independent, identically distributed random variables.The purpose of this paper is to consider an alternative approach to defining the random effects, in which the stratum means of the response variable are assumed to be independent and identically distributed, with a distribution not depending on beta . It will be shown that inferences about beta based on this formulation of the generalized linear mixed model have many desirable properties. For instance, inferences regarding beta are less sensitive to the choice of random effects distribution, are less subject to bias from omitted stratum-level covariates and are less affected by separate between- and within-cluster covariate effects. (C) 2009 Elsevier B.V. All rights reserved.