×

An alternative derivation of Harville’s restricted log likelihood function for variance component estimation. (English) Zbl 1412.62188

Summary: Estimation of variance components in linear mixed models is important in clinical trial and longitudinal data analysis. It is also important in animal and plant breeding for accurately partitioning total phenotypic variance into genetic and environmental variances. Restricted maximum likelihood (REML) method is often preferred over the maximum likelihood (ML) method for variance component estimation because REML takes into account the lost degree of freedom resulting from estimating the fixed effects. The original restricted likelihood function involves a linear transformation of the original response variable (a collection of error contrasts). Harville’s final form of the restricted likelihood function does not involve the transformation and thus is much easier to manipulate than the original restricted likelihood function. There are several different ways to show that the two forms of the restricted likelihood are equivalent. In this study, I present a much simpler way to derive Harville’s restricted likelihood function. I first treat the fixed effects as random effects and call such a mixed model a pseudo random model (PDRM). I then construct a likelihood function for the PDRM. Finally, I let the variance of the pseudo random effects be infinity and show that the limit of the likelihood function of the PDRM is the restricted likelihood function.

MSC:

62P10 Applications of statistics to biology and medical sciences; meta analysis
62H11 Directional data; spatial statistics
62J05 Linear regression; mixed models

Software:

SAS/STAT; SAS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Box, G. E. P., & Tiao, G. C. (1973). Bayesian inference in statistical analysis. New York: Wiley & Sons, Inc. · Zbl 0271.62044
[2] Hager, W. (1989). Updating the inverse of a matrix. SIAM Review, 31, 221-239. · Zbl 0671.65018
[3] Hartley, H. O., & Rao, J. N. K. (1967). Maximum‐likelihood estimation for the mixed analysis of variance model. Biometrika, 54, 93-108. · Zbl 0178.22001
[4] Harville, D. A. (1974). Bayesian inference for variance components using only error contrasts. Biometrika, 61, 383-385. · Zbl 0281.62072
[5] Harville, D. A. (1977). Maximum likelihood approaches to variance component estimation and to related problems. Journal of American Statistical Association, 72, 320-338. · Zbl 0373.62040
[6] Henderson, C. R. (1950). Estimation of genetic parameters (Abstract). The Annals of Mathematical Statistics, 21, 309-310.
[7] Henderson, C. R. (1953). Estimation of variance and covariance components. Biometrics, 9, 226-252.
[8] Henderson, C. R. (1975). Best linear unbiased estimation and prediction under a selection model. Biometrics, 31, 423-447. · Zbl 0335.62048
[9] Kang, H. M., Zaitlen, N. A., Wade, C. M., Kirby, A., Heckerman, D., Daly, M. J., & Eskin, E. (2008). Efficient control of population structure in model organism association mapping. Genetics, 178, 1709-1723.
[10] Patterson, H. D., & Thompson, R. (1971). Recovery of inter‐block information when block sizes are unequal. Biometrika, 58, 545-554. · Zbl 0228.62046
[11] Institute Inc, S. A. S. (2009). SAS/STAT: Users’ Guide, Version 9.3. Cary, NC: SAS Institute Inc.
[12] VanRaden, P. M. (2008). Efficient methods to compute genomic predictions. Journal of Dairy Science, 91, 4414-4423.
[13] Yu, J., Pressoir, G., Briggs, W. H., Vroh Bi, I., Yamasaki, M., Doebley, J. F., … Buckler, E. S. (2006). A unified mixed‐model method for association mapping that accounts for multiple levels of relatedness. Nature Genetics, 38, 203-208.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.