An alternative specification of generalized linear mixed models

Consider stratified data in which Yi1,…,YiniYi1,…,Yini denote real-valued response variables corresponding to the observations from stratum ii, i=1,…,mi=1,…,m and suppose that YijYij follows an exponential family distribution with canonical parameter of the form θij=xijβ+γiθij=xijβ+γ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 γ1,…,γmγ1,…,γ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 ββ. It will be shown that inferences about ββ based on this formulation of the generalized linear mixed model have many desirable properties. For instance, inferences regarding ββ 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.