Skip to main content

A general approach to mixed effects modeling of residual variances in generalized linear mixed models


We propose a general Bayesian approach to heteroskedastic error modeling for generalized linear mixed models (GLMM) in which linked functions of conditional means and residual variances are specified as separate linear combinations of fixed and random effects. We focus on the linear mixed model (LMM) analysis of birth weight (BW) and the cumulative probit mixed model (CPMM) analysis of calving ease (CE). The deviance information criterion (DIC) was demonstrated to be useful in correctly choosing between homoskedastic and heteroskedastic error GLMM for both traits when data was generated according to a mixed model specification for both location parameters and residual variances. Heteroskedastic error LMM and CPMM were fitted, respectively, to BW and CE data on 8847 Italian Piemontese first parity dams in which residual variances were modeled as functions of fixed calf sex and random herd effects. The posterior mean residual variance for male calves was over 40% greater than that for female calves for both traits. Also, the posterior means of the standard deviation of the herd-specific variance ratios (relative to a unitary baseline) were estimated to be 0.60 ± 0.09 for BW and 0.74 ± 0.14 for CE. For both traits, the heteroskedastic error LMM and CPMM were chosen over their homoskedastic error counterparts based on DIC values.

(To access the full article, please see PDF)

Author information



Corresponding author

Correspondence to Robert J Tempelman.

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Kizilkaya, K., Tempelman, R.J. A general approach to mixed effects modeling of residual variances in generalized linear mixed models. Genet Sel Evol 37, 31 (2005).

Download citation


  • Bayesian analysis
  • genetic evaluation
  • heterogeneous variances
  • threshold model