Robust estimation in multivariate heteroscedastic regression models with autoregressive covariance structures using EM algorithm

GÜNEY Y., ARSLAN O., GÖKALP YAVUZ F.

Journal of Multivariate Analysis, cilt.191, 2022 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 191
  • Basım Tarihi: 2022
  • Doi Numarası: 10.1016/j.jmva.2022.105026
  • Dergi Adı: Journal of Multivariate Analysis
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus, Academic Search Premier, ABI/INFORM, CAB Abstracts, MathSciNet, Psycinfo, zbMATH
  • Anahtar Kelimeler: Heteroscedastic regression models, Joint mean and scale covariance model, Laplace distribution, Modified decomposition, Prediction, Robustness, GENERALIZED ESTIMATING EQUATIONS, VARIABLE SELECTION, MATRIX, DISTRIBUTIONS, MIXTURES
  • Orta Doğu Teknik Üniversitesi Adresli: Evet

Özet

© 2022 Elsevier Inc.In the analysis of repeated or clustered measurements, it is crucial to determine the dynamics that affect the mean, variance, and correlations of the data, which will be possible using appropriate models. One of these models is the joint mean–covariance model, which is a multivariate heteroscedastic regression model with autoregressive covariance structures. In these models, parameter estimation is usually carried on under normality assumption, but the resulting estimators will be very sensitive to the outliers or non-normality of data. In this study, we propose a robust alternative method and an EM algorithm for estimating the parameters of joint mean–covariance models. Robustification is achieved using a multivariate heavy-tailed distribution with the same number of parameters as the multivariate normal distribution. To simplify the estimation procedure, a modified Cholesky decomposition is adopted to factorize the dependence structure in terms of unconstrained autoregressive and scale innovation parameters. Also, the technique for the prediction of future responses is given. An intensive simulation study and a real data example are provided to demonstrate the performance of the proposed method.