Simple reparameterization to improve convergence in linear mixed models

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...

Full description

Saved in:
Bibliographic Details
Main Authors: Gregor GORJANC, Tina FLISAR, Jose Carlos MARTÍNEZ-ÁVILA, Luis Alberto GARCÍA-CORTÉS
Format: Article
Language:English
Published: University of Ljubljana Press (Založba Univerze v Ljubljani) 2010-12-01
Series:Acta Agriculturae Slovenica
Subjects:
Online Access:https://journals.uni-lj.si/aas/article/view/14699
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:1854-1941