×

Robustness for general design mixed models using the \(t\)-distribution. (English) Zbl 07257703

Summary: The \(t\)-distribution allows the incorporation of outlier robustness into statistical models while retaining the elegance of likelihood-based inference. In this paper, we develop and implement a linear mixed model for the general design of the linear mixed model using the univariate \(t\)-distribution. This general design allows a considerably richer class of models to be fit than is possible with existing methods. Included in this class are semi-parametric regression and smoothing and spatial models.

MSC:

62-XX Statistics

Software:

SemiPar; WinBUGS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Aitkin M (1999) A general maximum likelihood analysis of variance components in generalized linear models . Biometrics, 55, 117-28 . · Zbl 1059.62564 · doi:10.1111/j.0006-341X.1999.00117.x
[2] Booth JG and Hobert JP (1999) Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm . Journal of the Royal Statistics Society, Series B, 61, 265-85 . · Zbl 0917.62058 · doi:10.1111/1467-9868.00176
[3] Breusch TS , Robertson JC and Welsh AH (1997) The emperor’s new clothes: a critique of the multivariate t regression model . Statistica Neerlandica, 51, 269-86 . · Zbl 0929.62062 · doi:10.1111/1467-9574.00055
[4] Caffo BS , Jank W and Jones GL (2005) Ascent-based Monte Carlo expectation maximization . Journal ofthe Royal Statistical Society, Series B: Statistical Methodology, 67, 235-51 . · Zbl 1075.65011 · doi:10.1111/j.1467-9868.2005.00499.x
[5] Crainiceanu C , Ruppert D and Wand MP (2005) Bayesian analysis for penalized spline regression using WinBUGS . Journal of Statistical Software, 14. · doi:10.18637/jss.v014.i14
[6] Dempster AP , Laird NM and Rubin DB (1977) Maximum likelihood from incomplete data via the EM algorithm . Journal ofthe Royal Statistics Society, Series B, 39, 1-22 . · Zbl 0364.62022
[7] Ghidey W , Lesaffre E and Eilers P (2004) Smooth random effects distribution in a linear mixed model . Biometrics, 60, 945-53 . · Zbl 1274.62238 · doi:10.1111/j.0006-341X.2004.00250.x
[8] Hagobian TA and Braun B (2006) Interactions between energy surplus and short-term exercise on glucose and insulin responses in healthy people with induced, mild insulin insensitivity . Metabolism, 55, 402-408 . · doi:10.1016/j.metabol.2005.09.017
[9] Henderson CR (1975) Best linear unbiased estimation and prediction under a selection model . Biometrics, 31, 423-47 . · Zbl 0335.62048 · doi:10.2307/2529430
[10] Jara A and Quintana F (2006) Unpublished MS. Linear mixed models with skew-elliptical distributions: a Bayesian approach.
[11] Kammann EE and Wand MP (2003) Geoadditive models . Journal ofthe Royal Statistical Society,SeriesC, 52, 1-18 . · Zbl 1111.62346 · doi:10.1111/1467-9876.00385
[12] Kleinman K and Ibrahim JG (1998) A semi-parametric Bayesian approach to the random effects model . Biometrics, 54, 921-38 . · Zbl 1058.62513 · doi:10.2307/2533846
[13] Kotz S and Nadarajah S (2004) Multivariate t-distributions and their applications. Cambridge: Cambridge University Press . · Zbl 1100.62059 · doi:10.1017/CBO9780511550683
[14] Laird NL and Ware JH (1982) Random-effects models for longitudinal data . Biometrics, 38, 963-74 . · Zbl 0512.62107 · doi:10.2307/2529876
[15] Lange KL , Little RJA and Taylor JMG (1989) Robust statistical modeling using the t-distribution . Journal of the American Statistical Association, 84, 881-96 .
[16] Levine RA and Casella G (2001) Implementations of the Monte Carlo EM algorithm . Journal ofComputational and Graphical Statistics, 10, 422-39 . · doi:10.1198/106186001317115045
[17] Lin X and Zhang D (1999) Inference in generalized additive mixed models by using smoothing splines . Journal ofthe Royal Statistical Society, Series B, 61, 381-400 . · Zbl 0915.62062 · doi:10.1111/1467-9868.00183
[18] Liu C , Rubin DB and Wu YN (1998) Parameter expansion to accelerate EM: The PX-EX algorithm . Biometrika, 85, 755-70 . · Zbl 0921.62071 · doi:10.1093/biomet/85.4.755
[19] Louis TA (1982) Finding the observed information matrix when using the EM algorithm . Journal of the Royal Statistical Society, Series B, 44, 226-33 . · Zbl 0488.62018
[20] McCulloch CE and Searle SR (2001) Generalized, linear, and mixed models. London: Chapman & Hall . · Zbl 0964.62061
[21] Meng XL and Rubin DB (1993) Maximum likelihood estimation via the ECM algorithm: a general framework . Biometrika, 80, 267-78 . · Zbl 0778.62022 · doi:10.1093/biomet/80.2.267
[22] O’Connell MA and Wolfinger RD (1997) Spatial regression models, response surfaces, and process optimization . Journal of Computational and Graphical Statistics, 6, 224-41 .
[23] Pinheiro JC , Liu C and Wu YN (2001) Efficient algorithms for robust estimation in linear mixed-effects models using the multivariate t-distribution . Journal ofComputational and Graphical Statistics, 10, 249-76 . · doi:10.1198/10618600152628059
[24] Robinson GK (1991) That BLUP is a good thing: the estimation of random effects . Statistical Science, 6, 15-51 . · Zbl 0955.62500 · doi:10.1214/ss/1177011926
[25] Rosa GJM , Gianola D and Padovani CR (2004) Bayesian longitudinal data analysis with mixed models and thick-tailed distributions using MCMC . Journal ofApplied Statistics, 31, 855-73 . · Zbl 1121.62479
[26] Rosa GJM , Padovani CR and Gianola D (2003) Robust linear mixed models with normal/independent distributions and Bayesian MCMC implementation . Biometrical Journal, 45, 573-90 . · Zbl 1441.62474 · doi:10.1002/bimj.200390034
[27] Ruppert D , Wand MP and Carroll RJ (2003) Semiparametric regression modelling. Cambridge: Cambridge University Press . · Zbl 1038.62042 · doi:10.1017/CBO9780511755453
[28] Smith M and Kohn R (1996) Nonparametric regression using Bayesian variable selection . Journal of Econometrics, 75, 317-44 . · Zbl 0864.62025 · doi:10.1016/0304-4076(95)01763-1
[29] Speed T (1991) Comment on paper by Robinson . Statistical Science, 6, 42-44 .
[30] Stranden I and Gianola D (1998) Mixed effects linear models with t-distributions forquantitative genetic analysis: a Bayesian approach . Genetics, Selection, Evolution, 31, 25-42 . · doi:10.1186/1297-9686-31-1-25
[31] Tao H , Palta M , Yandel B and Newton MA (1999) An estimation method for the semiparametric mixed effects model . Biometrics, 55, 102-110 . · Zbl 1059.62572 · doi:10.1111/j.0006-341X.1999.00102.x
[32] van Dyk DA (2000) Nesting EM algorithms for computational efficiency . Statistica Sinica, 10, 203-25 . · Zbl 1009.62018
[33] Wahba G (1978) Improper priors, spline smoothing and the problem of guarding against model errors in regression . Journal ofthe Royal Statistical Society, Series B, 40, 364-72 . · Zbl 0407.62048
[34] Wei GCG and Tanner MA (1990) A Monte Carlo implementation of the EM algorithm and the poor man s data augmentation algorithm . Journal ofthe American Statistical Association, 85, 699-704 . · doi:10.1080/01621459.1990.10474930
[35] Welsh AH and Richardson AM (1997) Approaches to the robust estimation of mixed models. In Maddala GS and Rao CR . (eds) Handbook ofStatistics, Volume 15, New York: Elsevier Science B.V . · Zbl 0908.62071 · doi:10.1016/S0169-7161(97)15015-5
[36] Wu CFJ (1983) On the convergence properties of the EM algorithm . The Annals ofStatistics, 11, 95-103 . · Zbl 0517.62035 · doi:10.1214/aos/1176346060
[37] Zhang D and Davidian M (2001) Linear mixed models with flexible distributions of random effects for longitudinal data . Biometrics, 57, 795-802 . · Zbl 1209.62087 · doi:10.1111/j.0006-341X.2001.00795.x
[38] Zhao Y , Staudenmayer J , Coull BA and Wand MP (2006) General design Bayesian generalized linear mixed models . Statistical Science, 21, 35-51 . · Zbl 1129.62063 · doi:10.1214/088342306000000015
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.