Mixed-effects models with random cluster sizes

Mixed effects models are often used in situations where responses are clustered. In this paper, we show that in the case of generalized linear mixed models where the cluster sizes are assumed to be independent random variables, whose joint distribution is unknown but does not depend on the parameters involved in the mixed model, the standard maximum likelihood estimation procedure assuming nonrandom cluster sizes can be applied without modification. Some extensions of the result are discussed.