Employing a Monte Carlo algorithm in Newton-type methods for restricted maximum likelihood estimation of genetic parameters

A - Papers appearing in refereed journals

Matilainen, K., Mantysaari, E. A., Lidauer, M. H., Stranden, I. and Thompson, R. 2013. Employing a Monte Carlo algorithm in Newton-type methods for restricted maximum likelihood estimation of genetic parameters. PLOS ONE. 8 (12), p. e80821.

AuthorsMatilainen, K., Mantysaari, E. A., Lidauer, M. H., Stranden, I. and Thompson, R.
Abstract

Estimation of variance components by Monte Carlo (MC) expectation maximization (EM) restricted maximum likelihood (REML) is computationally efficient for large data sets and complex linear mixed effects models. However, efficiency may be lost due to the need for a large number of iterations of the EM algorithm. To decrease the computing time we explored the use of faster converging Newton-type algorithms within MC REML implementations. The implemented algorithms were: MC Newton-Raphson (NR), where the information matrix was generated via sampling; MC average information(AI), where the information was computed as an average of observed and expected information; and MC Broyden's method, where the zero of the gradient was searched using a quasi-Newton-type algorithm. Performance of these algorithms was evaluated using simulated data. The final estimates were in good agreement with corresponding analytical ones. MC NR REML and MC AI REML enhanced convergence compared to MC EM REML and gave standard errors for the estimates as a by-product. MC NR REML required a larger number of MC samples, while each MC AI REML iteration demanded extra solving of mixed model equations by the number of parameters to be estimated. MC Broyden's method required the largest number of MC samples with our small data and did not give standard errors for the parameters directly. We studied the performance of three different convergence criteria for the MC AI REML algorithm. Our results indicate the importance of defining a suitable convergence criterion and critical value in order to obtain an efficient Newton-type method utilizing a MC algorithm. Overall, use of a MC algorithm with Newton-type methods proved feasible and the results encourage testing of these methods with different kinds of large-scale problem settings.

Year of Publication2013
JournalPLOS ONE
Journal citation8 (12), p. e80821
Digital Object Identifier (DOI)doi:10.1371/journal.pone.0080821
Open accessPublished as ‘gold’ (paid) open access
Publisher's version
Output statusPublished
Publication dates
Online10 Dec 2013
Publication process dates
Accepted09 Dec 2013
Copyright licenseCC BY
ISSN1932-6203
PublisherPublic Library of Science (PLOS)

Permalink - https://repository.rothamsted.ac.uk/item/8qyv7/employing-a-monte-carlo-algorithm-in-newton-type-methods-for-restricted-maximum-likelihood-estimation-of-genetic-parameters

17 total views
4 total downloads
0 views this month
0 downloads this month
Download files as zip