Simple reparameterization to improve convergence in linear mixed models | Enostavna reparametrizacija za izboljsanje konvergence linearnih mesanih modelov
2010
Gorjanc, G., Ljubljana Univ. (Slovenia). Biotechnical Fac., Zootechnical Dept. | Flisar, T., Ljubljana Univ. (Slovenia). Biotechnical Fac., Zootechnical Dept. | Martínez-Ávila, J.K., Departamento de Mejora Genética, Instituto Nacional de Investigación Agraria, Madrid (Spain) | García-Cortés, L.A., Departamento de Mejora Genética, Instituto Nacional de Investigación Agraria, Madrid (Spain)
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.
显示更多 [+] 显示较少 [-]