[1]    Aguilar I., Legarra A., Cardoso F., Masuda Y., Lourenco D., Misztal I. Frequentist p-values for large-scale-single step genome-wide association, with an application to birth weight in American Angus cattle. Genet. Sel. Evol. 51 (2019) 28. doi: 10.1186/s12711-019-0469-3.

[2]    Amestoy P.R., Davis T.A., Duff I.S. An approximate minimum degree ordering algorithm. SIAM J. Matr. Anal. Appl. 17 (1996) 886–905.

[3]    Anderson E., Bai Z., Bischof C., Blackford S., Demmel J., Dongarra J., Du Croz J., Greenbaum A., Hammarling S., McKenney A., Sorensen D. LAPACK Users’ Guide. Society for Industrial and Applied Mathematics, Philadelphia, PA, Third edn. (1999). ISBN 978-0-89871-447-0.

[4]    Bijma P. Multilevel selection 4: Modeling the relationship of indirect genetic effects and group size. Genetics 186 (2010) 1029–1031. doi: 10.1534/genetics.110.120485.

[5]    Bondari K., Willham R.L., Freeman A.E. Estimates of direct and maternal genetic correlations for pupa weight and family size of Tribolium. J. Anim. Sci. 47 (1978) 358–365.

[6]    Boyd S., Vandenberghe L. Convex Optimization. Cambridge University Press (2004).

[7]    Christensen O.F. Compatibility of pedigree-based and marker-based relationship matrices for single-step genetic evaluation. Genet. Sel. Evol. 44 (2012) 37. doi: 10.1186/1297-9686-44-37.

[8]    Dennis J.E., Schnabel R.B. Numerical Methods for Unconstrained Optimization and Nonlinear Equations. SIAM Classics in Applied Mathematics. Society for Industrial and Applied Mathematics, Philadelphia (1996).

[9]    Dongarra J.J., Croz J.D., Hammarling S., Hanson R.J. An extended set of FORTRAN Basic Linear Algebra Subprograms. ACM Trans. Math. Softw. 14 (1988) 1–17. doi: 10.1145/42288.42291.

[10]    Erisman A.M., Tinney W.F. On computing certain elements of the inverse of a sparse matrix. Commun. ACM 18 (1975) 177–179. ISSN 0001-0782. doi: 10.1145/360680.360704.

[11]    Eskow E., Schnabel R.B. Algorithm 695: Software for a new modified Cholesky factorization. ACM Trans. Math. Software 17 (1991) 306–312.

[12]    Fernando R.L., Grossman M. Genetic evaluation with autosomal and X-chromosomal inheritance. Theor. Appl. Genet. 80 (1990) 75–80. doi: 10.1007/BF00224018.

[13]    Fernando R.L., Cheng H., Golden B.L., Garrick D.J. Computational strategies for alternative single-step Bayesian regression models with large numbers of genotyped and non-genotyped animals. Genet. Sel. Evol. 48 (2016) 96. doi: 10.1186/s12711-016-0273-2.

[14]    Ferris M., Lucidi S., Roma M. Nonmonotone curvilinear line search methods for unconstrained optimization. Comput. Optim. Applic. 6 (1996) 117–136.

[15]    Forsgren A., Gill P.E., Murray W. Computing modified Newton directions using a partial Cholesky factorization. SIAM J. Sci. Statist. Comp. 16 (1995) 139–150.

[16]    Garcia-Baccino C.A., Legarra A., Christensen O.F., Misztal I., Pocrnic I., Vitezica Z.G., Cantet R.J.C. Metafounders are related to F_ST fixation indices and reduce bias in single-step genomic evaluations. Genet. Sel. Evol. 49 (2017) 34. doi: 10.1186/s12711-017-0309-2.

[17]    Grippo L., Lampariello F., Lucidi S. A nonmonotone line search technique for Newton’s method. SIAM J. Numer. Anal. 23 (1986) 707–716. ISSN 0036-1429. doi: 10.1137/0723046.

[18]    Gualdrón Duarte J.L., Cantet R.J.C., Bates R.O., Ernst C.W., Raney N.E., Steibel J.P. Rapid screening for phenotype-genotype associations by linear transformations of genomic evaluations. BMC Bioinformatics 15 (2014) 246. doi: 10.1186/1471-2105-15-246.

[19]    Gustavson F., Waśniewski J., Dongarra J., Langou J. Rectangular full packed format for Cholesky’s algorithm: factorization, solution, and inversion. ACM Trans. Math. Softw. 37 (2010) 18. doi: 10.1145/1731022.1731028.

[20]    Harville D.A. Use of the Gibbs sampler to invert large, possibly sparse, positive definite matrices. Lin. Alg. Appl. 289 (1999) 203–224.

[21]    Ihaka R., Gentleman R. R: A language for data analysis and graphics. J. Comp. Graph. Stat. 5 (1996) 299–314.

[22]    Karypis G., Kumar V. MeTis A software package for partioning unstructured graphs, partitioning meshes, and computing fill-in reducing ordering of sparse matrices Version 4.0. Department of Computer Science, University of Minnesota, Minneapolis, MN 55455 (1998). 44 pp.

[23]    Koivula M., Negussie E., Mäntysaari E.A. Genetic parameters for test-day somatic cell count at different lactation stages of Finnish dairy cattle. Livest. Prod. Sci. 90 (2004) 145–157. doi: 10.1016/j.livprodsci.2004.03.004.

[24]    Legarra A., Aguilar I., Colleau J. Short communication: Methods to compute genomic inbreeding for ungenotyped individuals. J. Dairy Sci. 103 (2020) 3363–3367. doi: 10.3168/jds.2019-17750.

[25]    Legarra A., Christensen O.F., Vitezica Z.G., Aguilar I., Misztal I. Ancestral relationships using metafounders: finite ancestral populations and across population relationships. Genetics 200 (2015) 455–468. doi: 10.1534/genetics.115.177014.

[26]    Liu J.W.H. Modification of the minimum degree algorithm by multiple elimination. ACM Trans. Math. Soft. 11 (1985) 141–153.

[27]    Mäntysaari E.A. Derivation of multiple trait reduced random regression (RR) model for the first lactation test day records of milk, protein and fat. In: Proceedings of the 50th Annual Meeting of the European Association of Animal Production. Europ. Ass. Anim. Prod. (1999).

[28]    Meuwissen T.H.E., Luo Z. Computing inbreeding coefficients in large populations. Genet. Sel. Evol. 24 (1992) 305–313.

[29]    Meyer K. DfReml — a set of programs to estimate variance components under an individual animal model. In: Proceedings Animal Model Workshop, Vol. 71 Supplement 2 of J. Dairy Sci. Edmonton, Canada, June 25–26, 1988 (1988), pp. 33–34. doi: 10.1016/S0022-0302(88)79977-4.

[30]    Meyer K. DfReml version 3.0. CD-ROM of the Sixth World Congress on Genetics Applied to Livestock Production (1998).

[31]    Meyer K. Ordering strategies to reduce computational requirements in variance component estimation. Proc. Ass. Advan. Anim. Breed. Genet. 16 (2005) 282–285.

[32]    Meyer K. Random regression analyses using B–splines to model growth of Australian Angus cattle. Genet. Sel. Evol. 37 (2005) 473–500. doi: 10.1051/gse:2005012.

[33]    Meyer K. Multivariate analyses of carcass traits for Angus cattle fitting reduced rank and factor–analytic models. J. Anim. Breed. Genet. 124 (2007) 50–64. doi: 10.1111/j.1439-0388.2007.00637.x.

[34]    Meyer K. A penalized likelihood approach to pooling estimates of covariance components from analyses by parts. J. Anim. Breed. Genet. 130 (2013) 270–285. doi: 10.1111/jbg.12004.

[35]    Meyer K., Houle D. Sampling based approximation of confidence intervals for functions of genetic covariance matrices. Proc. Ass. Advan. Anim. Breed. Genet. 20 (2013) 523–526.

[36]    Meyer K., Kirkpatrick M. Restricted maximum likelihood estimation of genetic principal components and smoothed covariance matrices. Genet. Sel. Evol. 37 (2005) 1–30. doi: 10.1051/gse:2004034.

[37]    Meyer K., Kirkpatrick M. Better estimates of genetic covariance matrices by ‘bending’ using penalized maximum likelihood. Genetics 185 (2010) 1097–1110. doi: 10.1534/genetics.109.113381.

[38]    Meyer K., Kirkpatrick M., Gianola D. Penalized maximum likelihood estimates of genetic covariance matrices with shrinkage towards phenotypic dispersion. Proc. Ass. Advan. Anim. Breed. Genet. 19 (2011) 87–90.

[39]    Meyer K., Smith S.P. Restricted maximum likelihood estimation for animal models using derivatives of the likelihood. Genet. Sel. Evol. 28 (1996) 23–49. doi: 10.1051/gse:19960102.

[40]    Meyer K., Tier B. "SNP Snappy": A strategy for fast genome wide association studies fitting a full mixed model. Genetics 190 (2012) 275–277. doi: 10.1534/genetics.111.134841.

[41]    Nelder J.A., Mead R. A simplex method for function minimization. Computer J. 7 (1965) 308–313.

[42]    Nocedahl J., Wright S.J. Numerical Optimization. Springer Series in Operations Research. Springer Verlag, New York, Berlin Heidelberg (1999). ISBN 0-38798793-2.

[43]    Powell M.J.D. An efficient method for finding the minimum of a function of several variables without calculating derivatives. Computer J. 7 (1965) 155–162.

[44]    Quaas R.L. Computing the diagonal elements of a large numerator relationship matrix. Biometrics 32 (1976) 949–953. doi: 10.2307/2529279.

[45]    Schnabel R.B., Estrow E. A new modified Cholesky factorization. SIAM J. Sci. Statist. Comp. 11 (1990) 1136–1158.

[46]    Schnabel R.B., Estrow E. A revised modified Cholesky factorization algorithm. SIAM J. Opt. 9 (1999) 1135–1149.

[47]    Tier B. Computing inbreeding coefficients quickly. Genet. Sel. Evol. 22 (1990) 419–425.

[48]    Westell R.A., Quaas R.L., Van Vleck L.D. Genetic groups in an animal model. J. Dairy Sci. 71 (1988) 1310–1318. doi: 10.3168/jds.S0022-0302(88)79688-5.

[49]    Wilson A.J., Reale D., Clements M.N., Morrissey M.B., Postma E., Walling C.A., Kruuk L.E.B., Nussey D.H. An ecologist’s guide to the animal model. Journal of Animal Ecology 79 (2010) 13–26. doi: 10.1111/j.1365-2656.2009.01639.x.