+ 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

Application of supernodal sparse factorization and inversion to the estimation of (co)variance components by residual maximum likelihood



Application of supernodal sparse factorization and inversion to the estimation of (co)variance components by residual maximum likelihood



Journal of Animal Breeding and Genetics 131(3): 227-236



We demonstrated that supernodal techniques were more efficient than traditional methods for factorization and inversion of a coefficient matrix of mixed model equations (MME), which are often required in residual maximum likelihood (REML). Supernodal left-looking and inverse multifrontal algorithms were employed for sparse factorization and inversion, respectively. The approximate minimum degree or multilevel nested dissection was used for ordering. A new computer package, Yet Another MME Solver (YAMS), was developed and compared with FSPAK with respect to computing time and size of temporary memory for 13 test matrices. The matrices were produced by fitting animal models to dairy data and by using simulations from sire, sire-maternal grand sire, maternal and dominance models for phenotypic data and animal model for genomic data. The order of matrices ranged from 32,840 to 1,048,872. The YAMS software factorized and inverted the matrices up to 13 and 10 times faster than FSPAK, respectively, when an appropriate ordering strategy was applied. The YAMS package required at most 282 MB and 512 MB of temporary memory for factorization and inversion, respectively. Processing time per iteration in average information REML was reduced, using YAMS. The YAMS package is freely available on request by contacting the corresponding author.

Please choose payment method:






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

Accession: 051613381

Download citation: RISBibTeXText

PMID: 24906028

DOI: 10.1111/jbg.12058


Related references

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 of variance components in animal model using sparse matrix inversion and a supercomputer. Journal of Dairy Science 72(Suppl. 1): 30, 1989

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

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 Variance Components for Univariate Animal Models Using Sparse Matrix Techniques and Average Information. Journal of Dairy Science 78(2): 449-456, 1995

Application of singular value decomposition to restricted maximum likelihood estimation of variance components. Journal of Dairy Science 70(12): 2680-2684, 1987

Restricted maximum likelihood estimation of variance components for multiple traits with missing observations and an application to beef cattle. Dissertation Abstracts International B, Sciences and Engineering 49(10): 4071, 1989

Restricted maximum likelihood method to variance components estimation. Ciencia e Pratica 18(4): 425-437, 1994

Maximum likelihood estimation of variance components using human pedigrees. Genetics 88(4 Part 2): S68-S69, 1978

Estimation of variance components with unbalanced data in a maternal effects model ii. maximum likelihood estimators of variance components for direct and maternal effects of 8 week weights of large white pigs. Journal of Animal Breeding and Genetics 105(6): 462-467, 1988

Property of trace in restricted maximum likelihood estimation of variance components. Journal of Dairy Science 73(12): 3583-3585, 1990

A Monte-Carlo algorithm for maximum likelihood estimation of variance components. Genetics Selection Evolution 28(4): 329-343, 1996

Multiple lactation variance components estimation using restricted maximum likelihood. Dissertation Abstracts International B, Sciences and Engineering 49(6): 2007, 1988

A universal formula of maximum likelihood estimation of variance-covariance components. Journal of Geodesy 70(4): 233-240, 1996

Maximum likelihood estimation of variance components in dairy cattle breeding research. Journal of Dairy Science 59(12): 2146-2150, 1976