Efficient ReML inference in variance component mixed models using a Min-Max algorithm

PLoS Comput Biol. 2022 Jan 24;18(1):e1009659. doi: 10.1371/journal.pcbi.1009659. eCollection 2022 Jan.

Abstract

Since their introduction in the 50's, variance component mixed models have been widely used in many application fields. In this context, ReML estimation is by far the most popular procedure to infer the variance components of the model. Although many implementations of the ReML procedure are readily available, there is still need for computational improvements due to the ever-increasing size of the datasets to be handled, and to the complexity of the models to be adjusted. In this paper, we present a Min-Max (MM) algorithm for ReML inference and combine it with several speed-up procedures. The ReML MM algorithm we present is compared to 5 state-of-the-art publicly available algorithms used in statistical genetics. The computational performance of the different algorithms are evaluated on several datasets representing different plant breeding experimental designs. The MM algorithm ranks among the top 2 methods in almost all settings and is more versatile than many of its competitors. The MM algorithm is a promising alternative to the classical AI-ReML algorithm in the context of variance component mixed models. It is available in the MM4LMM R-package.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms*
  • Computational Biology / methods*
  • Models, Genetic*
  • Models, Statistical*

Grants and funding

This work was funded by the French National Research Agency (ANR) as part of the “Investissements d’Avenir” program, ANR-10-BTBR-0001, Amaizing (https://anr.fr) to FL, AC, and TMH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.