+ Site Statistics
+ Search Articles
+ PDF Full Text Service
How our service works
Request PDF Full Text
+ Follow Us
Follow on Facebook
Follow on Twitter
Follow on LinkedIn
+ Subscribe to Site Feeds
Most Shared
PDF Full Text
+ Translate
+ Recently Requested

Derivative-free restricted maximum likelihood estimation in animal models with a sparse matrix solver



Derivative-free restricted maximum likelihood estimation in animal models with a sparse matrix solver



Journal of Dairy Science 74(12): 4337-4343



Estimation of (co)variance components by derivative-free REML requires repeated evaluation of the log-likelihood function of the data. Gaussian elimination of the augmented mixed model coefficient matrix is often used to evaluate the likelihood function, but it can be costly for animal models with large coefficient matrices. This study investigated the use of a direct sparse matrix solver to obtain the log-likelihood function. The sparse matrix package SPARSPAK was used to reorder the mixed model equations once and then repeatedly to solve the equations by Cholesky factorization to generate the terms required to calculate the likelihood. The animal model used for comparison contained 19 fixed levels, 470 maternal permanent environmental effects, and 1586 direct and 1586 maternal genetic effects, resulting in a coefficient matrix of order 3661 with .3% nonzero elements after including numerator relationships. Compared with estimation via Gaussian elimination of the unordered system, utilization of SPARSPAK required 605 and 240 times less central processing unit time on mainframes and personal computers, respectively. The SPARSPAK package also required less memory and provided solutions for all effects in the model.

Please choose payment method:






(PDF emailed within 0-6 h: $19.90)

Accession: 002066494

Download citation: RISBibTeXText

PMID: 1787202

DOI: 10.3168/jds.s0022-0302(91)78629-3


Related references

Restricted Maximum Likelihood Estimation of Variance Components for Univariate Animal Models Using Sparse Matrix Techniques and Average Information. Journal of Dairy Science 78(2): 449-456, 1995

Restricted maximum likelihood estimation of variance components in animal model using sparse matrix inversion and a supercomputer. Journal of Dairy Science 73(1): 163-172, 1990

Restricted maximum likelihood estimation of variance components in animal model using sparse matrix inversion and a supercomputer. Journal of Dairy Science 72(Suppl. 1): 30, 1989

A derivative free approach for estimating variance components in animal models by restricted maximum likelihood. Journal of Animal Science 64(5): 1362-1370, 1987

Univariate and multivariate animal models for evaluation of Sahiwal sires using derivative-free restricted maximum likelihood method. Indian Journal of Animal Sciences 82(5): 501-504, 2012

Restricted maximum likelihood to estimate variance components for animal models with several random effects using a derivative free algorithm. Genetics Selection Evolution 21(3): 317-340, 1989

Restricted maximum likelihood estimation of covariances in sparse linear models. Genetics Selection Evolution (Paris) 30(1): 3-26, 1998

Sparse matrix inversion for restricted maximum likelihood estimation of variance components by expectation-maximization. Journal of dairy science 76(5): 1479-1483, 1993

Restricted maximum likelihood estimation for animal models using derivatives of the likelihood. Genetics Selection Evolution 28(1): 23-49, 1996

Computation of all eigenvalues of matrices used in restricted maximum likelihood estimation of variance components using sparse matrix techniques. Genetics Selection Evolution 28(1): 51-65, 1996

Restricted maximum likelihood estimation of genetic parameters for number of pigs born alive estimated from field data using sparse matrix inversion. Journal of Animal Science 70(Suppl. 1): 147, 1992

Sequential transformation for multiple traits for estimation of (co)variance components with a derivative-free algorithm for restricted maximum likelihood. Journal of Animal Science 71(4): 836-844, 1993

Estimation of components of genetic variance and heritability for flowering time and yield in gerbera using Derivative-Free Restricted Maximum Likelihood (DFRML). Tag. Theoretical and Applied Genetics. Theoretische und Angewandte Genetik 86(2-3): 234-236, 1993

Estimation of genetic parameters by restricted maximum likelihood under multi-trait animal models in selected layer lines. Japanese Poultry Science 33(3): 185-192, 1996

Evaluation of derivative-free marginal maximum likelihood estimation of variance components in threshold mixed models. Proceedings, 5th World Congress, University of Guelph, Guelph, Ontario, Canada, 7-12 August, 1994 Volume 22 Computing strategies and software: 15-18, 1994