Employing a Monte Carlo algorithm in expectation maximization restricted maximum likelihood estimation of the linear mixed model

Journal of Animal Breeding and Genetics

Volumn 129, Issue 6, 2012, Pages 457-468

  Matilainen, K ^a   Mantysaari E A ^a   Lidauer, M H ^a   Stranden I ^a   Thompson, R ^b  

^a FINNISH FOREST RESEARCH INSTITUTE   (Finland)

^b ROTHAMSTED RESEARCH   (United Kingdom)

Author keywords

Convergence criterion; Monte Carlo sample size; Variance components

Indexed keywords

ALGORITHM; ANIMAL; ARTICLE; BREEDING; CATTLE; FEASIBILITY STUDY; MONTE CARLO METHOD; STATISTICAL MODEL;

ALGORITHMS; ANIMALS; BREEDING; CATTLE; FEASIBILITY STUDIES; LIKELIHOOD FUNCTIONS; LINEAR MODELS; MONTE CARLO METHOD;

EID: 84869166290     PISSN: 09312668     EISSN: 14390388     Source Type: Journal    
DOI: 10.1111/j.1439-0388.2012.01000.x     Document Type: Article

References (27)

1. Booth J.G., Hobert J.P. (1999) Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm. J. R. Statist. Soc. B, 61, 265-285.
2. Celeux G., Diebolt J. (1985) The SEM algorithm: a probabilistic teacher algorithm derived from the EM algorithm for the mixture problem. Comp. Statist. Q., 2, 73-82.
3. Delyon B., Lavielle M., Moulines E. (1999) Convergence of a stochastic approximation version of the EM algorithm. Ann. Statist., 27, 94-128.
4. Dempster A.P., Laird N.M., Rubin D.B. (1977) Maximum likelihood from incomplete data via the EM algorithm. J. R. Statist. Soc. B, 39, 1-38.
5. Foulley J.-L., van Dyk D.A. (2000) The PX-EM algorithm for fast stable fitting of Henderson's mixed model. Genet. Sel. Evol., 32, 143-163.
6. García-Cortés L.A., Sorensen D. (2001) Alternative implementations of Monte Carlo EM algorithms for likelihood inference. Genet. Sel. Evol., 33, 443-452.
7. García-Cortés L.A., Moreno C., Varona L., Altarriba J. (1992) Variance component estimation by resampling. J. Anim. Breed. Genet., 109, 358-363.
8. García-Cortés L.A., Moreno C., Varona L., Altarriba J. (1995) Estimation of prediction-error variances by resampling. J. Anim. Breed. Genet., 112, 176-182.
9. Gilmour A., Thompson R., Cullis B. (1995) Average information REML: an efficient algorithm for variance parameter estimation in linear mixed models. Biometrics, 51, 1440-1450.
10. Guo S.W., Thompson E.A. (1991) Monte Carlo estimation of variance component models for large complex pedigrees. IMA J. Math. Appl. Med. Biol., 8, 171-189.
11. Harville D.A. (2004) Making REML computationally feasible for large data sets: use of Gibbs sampler. J. Stat. Comp. Sim., 74, 135-153.
12. Henderson C.R. (1984) Estimation of variances and covariances under multiple trait model. J. Dairy Sci., 67, 1581-1589.
13. Hickey J.M., Veerkamp R.F., Calus M.P.L., Mulder H.A., Thompson R. (2009) Estimation of prediction error variances via Monte Carlo sampling methods using different formulations of the prediction error variance. Genet. Sel. Evol., 41, 23.
14. Jaffrézic F., Meza C., Lavielle M., Foulley J.-L. (2006) Genetic analysis of growth curves using the SAEM algorithm. Genet. Sel. Evol., 38, 583-600.
15. Klassen D.J., Smith S.P. (1990) Animal model estimation using simulated REML. In: Proceedings 4th World Congress on Genetics Applied to Livestock Production XIII: 472-475. Edinburgh, 23-27 July 1990.
16. Levine R.A., Casella G. (2001) Implementations of the Monte Carlo EM algorithm. J. Comp. Graph. Stat., 10, 422-439.
17. Levine R.A., Fan J. (2004) An automated (Markov chain) Monte Carlo EM algorithm. J. Stat. Comp. Sim., 74, 349-360.
18. Lidauer M.H., Matilainen K., Mäntysaari E.A., Strandén I. (2011) Technical reference guide for MiX99, Release IV/2011. MTT Agrifood Research Finland. http://www.mtt.fi/BGE/Software/MiX99.
19. Madsen P., Jensen J. (2000) A user's guide to DMU. A Package for Analysing Multivariate Mixed Models, Ver. 6, release 4. Danish Institute of Agricultural Sciences, Denmark.
20. Mäntysaari E., Van Vleck L.D. (1989) Restricted maximum likelihood estimates of variance components from multitrait sire models with large number of fixed effects. J. Anim. Breed. Genet., 106, 409-422.
21. McCulloch C.E. (1997) Maximum likelihood algorithms for generalized linear mixed models. J. Am. Stat. Ass., 92, 162-170.
22. Meza C., Jaffrézic F., Foulley J.-L. (2009) Estimation in the probit normal model for binary outcomes using the SAEM algorithm. Comp. Stat. Data Anal., 53, 1350-1360.
23. Misztal I. (2008) Reliable computing in estimation of variance components. J. Anim. Breed. Genet., 125, 363-370.
24. Ripatti S., Larsen K., Palmgren J. (2002) Maximum likelihood inference for multivariate frailty models using an automated Monte Carlo EM algorithm. Lifetime Data Anal., 8, 349-360.
25. Strandén I., Lidauer M. (1999) Solving large mixed linear models using preconditioned conjugate gradient iteration. J. Dairy Sci., 82, 2779-2787.
26. Thompson R. (1994) Integrating best linear unbiased prediction and maximum likelihood estimation. In: Proceedings 5th World Congress on Genetics Applied to Livestock Production XVIII: 337-340. Guelph, 7-12 August 1994.
27. Wei G., Tanner M. (1990) A Monte Carlo implementation of the EM algorithm and the poor man's data augmentation algorithms. J. Am. Stat. Ass., 85, 699-704.