Genet. Sel. Evol.
Volume 40, Number 4, July-August 2008
|Page(s)||395 - 413|
|Published online||17 June 2008|
Bayes factor between Student t and Gaussian mixed models within an animal breeding contextJoaquim Casellas1, Noelia Ibáñez-Escriche1, Luis Alberto García-Cortés2 and Luis Varona1
1 Genètica i Millora Animal, IRTA-Lleida, 25198 Lleida, Spain
2 Departamento de Mejora Genética Animal, SGIT-INIA, Carretera de la Coruña, km. 7, 28040 Madrid, Spain
Received 2 April 2007; accepted 19 December 2007; published online 17 June 2008
Abstract - The implementation of Student t mixed models in animal breeding has been suggested as a useful statistical tool to effectively mute the impact of preferential treatment or other sources of outliers in field data. Nevertheless, these additional sources of variation are undeclared and we do not know whether a Student t mixed model is required or if a standard, and less parameterized, Gaussian mixed model would be sufficient to serve the intended purpose. Within this context, our aim was to develop the Bayes factor between two nested models that only differed in a bounded variable in order to easily compare a Student t and a Gaussian mixed model. It is important to highlight that the Student t density converges to a Gaussian process when degrees of freedom tend to infinity. The twomodels can then be viewed as nested models that differ in terms of degrees of freedom. The Bayes factor can be easily calculated from the output of a Markov chain Monte Carlo sampling of the complex model (Student t mixed model). The performance of this Bayes factor was tested under simulation and on a real dataset, using the deviation information criterion (DIC) as the standard reference criterion. The two statistical tools showed similar trends along the parameter space, although the Bayes factor appeared to be the more conservative. There was considerable evidence favoring the Student t mixed model for data sets simulated under Student t processes with limited degrees of freedom, and moderate advantages associated with using the Gaussian mixed model when working with datasets simulated with 50 or more degrees of freedom. For the analysis of real data (weight of Pietrain pigs at six months), both the Bayes factor and DIC slightly favored the Student t mixed model, with there being a reduced incidence of outlier individuals in this population.
Key words: Bayes factor / Gaussian distribution / mixed model / Student t distribution / preferential treatment
Corresponding author: Joaquim.Casellas@irta.es
© INRA, EDP Sciences 2008