Simple reparameterization to improve convergence in linear mixed models
DOI:
https://doi.org/10.14720/aas.2010.96.2.14699Keywords:
statistics, mixed model, Bayesian analysis, McMC, reparameterization, convergenceAbstract
Slow convergence and mixing are one of the main problems of Markov chain Monte Carlo (McMC) algorithms applied to mixed models in animal breeding. Poor convergence is to a large extent caused by high posterior correlation between variance components and solutions for the levels of associated effects. A simple reparameterization of the conventional model for variance component estimation is presented which improves McMC sampling and provides the same posterior distributions as the conventional model. Reparameterization is based on the rescaling of hierarchical (random) effects in a model, which alleviates posterior correlation. The developed model is compared against the conventional model using several simulated data sets. Results show that presented reparameterization has better behaviour of associated sampling methods and is several times more efficient for the low values of heritability.
Downloads
Published
Issue
Section
License
Copyright (c) 2010 University of Ljubljana, Biotechnical Faculty

This work is licensed under a Creative Commons Attribution 4.0 International License.