Correcting for nonrandom ascertainment in generalized linear mixed models (GLMMs), fitted using Gibbs sampling