Maximum Likelihood Estimation in Shape and Skew Scale Mixtures of Normal Linear Regression Models: The Impact of Reparameterization

Document Type : Original Article

Authors
1 Department of Statistics, Marv.C., Islamic Azad University, Marvdasht, Iran
2 Department of Statistics, Ra.C., Islamic Azad University, Rasht, Iran
3 Department of Statistics, Das.C., Islamic Azad University, Borazjan, Iran
Abstract
This study introduces a novel class of Linear Regression Models (LRM), the Shape and Skew Scale Mixtures of Normal (SHSSMN) distributions, to effectively handle skewed and heavy-tailed error terms. By extending the Expectation-Maximization (EM) algorithm, we explore two distinct scenarios based on without or with reparametrization. Simulation results demonstrate that the parameter estimates under both scenarios are consistent and nearly identical, with the second scenario offering a computational advantage by reducing the number of iterations required for convergence, especially for complex distributions. The models are applied to real-world data from the Australian Institute of Sport (AIS) to evaluate their practical utility. The results show that the SHSSMN models, particularly the Modified Skew-t Cauchy (STC) distribution, provide superior fit and effectively address the skewness and heavy-tailed nature of the residuals, making them a valuable tool for regression analysis in diverse fields.
Keywords
Subjects

Alizadeh Ghajari Z, Zare K, Shokri S. Estimation in shape mixtures of skew-normal linear regression models via ECM coupled with Gibbs sampling. Monte Carlo Methods and Applications. 2024;30(2):137–148.
Arellano-Valle RB, Castro LM, Genton MG, Gómez HW. Bayesian inference for shape mixtures of skewed distributions, with application to regression analysis. 2008;.
Arellano-Valle RB, Ferreira CS, Genton MG. Scale and shape mixtures of multivariate skew-normal distributions. Journal of Multivariate Analysis. 2018;166:98–110.
Arellano-Valle RB, Genton MG, Loschi RH. Shape mixtures of multivariate skewnormal distributions. Journal of Multivariate Analysis. 2009;100(1):91–101.
Arellano-Valle RB, Gómez HW, Quintana FA. A new class of skew-normal distributions.Communications in statistics-Theory and Methods. 2004;33(7):1465–1480.
Arnold BC, Beaver RJ. The skew-Cauchy distribution. Statistics & probability letters. 2000;49(3):285–290.
Ashour SK, Abdel-hameed MA. Approximate skew normal distribution. Journal of Advanced Research. 2010;1(4):341–350.
Azzalini A. A class of distributions which includes the normal ones. Scandinavian journal of statistics. 1985;p. 171–178.
Azzalini A, Capitanio A. Statistical applications of the multivariate skew normal distribution. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 1999;61(3):579–602.
Azzalini A, Capitanio A. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t-distribution. Journal of the Royal Statistical Society Series B: Statistical Methodology. 2003;65(2):367–389.
Branco MD, Dey DK. A general class of multivariate skew-elliptical distributions. Journal of Multivariate Analysis. 2001;79(1):99–113.
Celeux G. The SEM algorithm: a probabilistic teacher algorithm derived from the EM algorithm for the mixture problem. Computational statistics quarterly. 1985;2:73–82.
Cook RD, Weisberg S. An introduction to regression graphics, vol. 405. John Wiley & Sons; 2009.
Delyon B, Lavielle M, Moulines E. Convergence of a stochastic approximation version of the EM algorithm. Annals of statistics. 1999;p. 94–128.
Dempster AP, Laird NM, Rubin DB. Maximum likelihood from incomplete data via the EM algorithm. Journal of the royal statistical society: series B (methodological). 1977;39(1):1–22.
Ferreira C. Diagnostics Analysis for Skew Normal Regression Models. Instituto de Matemática, Estatística e Computação Científica, Universidade; 2008.
Ferreira C, Bolfarine H, Lachos VH. Skew scale mixtures of normal distributions: Properties and estimation. Statistical Methodology. 2011;8(2):154–171.
Ferreira CS, Lachos VH, Bolfarine H. Inference and diagnostics in skew scale mixtures of normal regression models. Journal of Statistical Computation and Simulation. 2015;85(3):517–537.
Genton MG, Loperfido NM. Generalized skew-elliptical distributions and their quadratic forms. Annals of the Institute of Statistical Mathematics. 2005;57(2):389–401.
Gómez HW, Venegas O, Bolfarine H. Skew-symmetric distributions generated by the distribution function of the normal distribution. Environmetrics: The official journal of the International Environmetrics Society. 2007;18(4):395–407.
Gupta AK, Chang F, Huang WJ. Some skew-symmetric models. Random Operators & Stochastic Equations. 2002;10(2).
Kahrari F, Rezaei M, Yousefzadeh F, Arellano-Valle RB. On the multivariate skewnormal-Cauchy distribution. Statistics & Probability Letters. 2016;117:80–88.
Kuhn E, Lavielle M. Coupling a stochastic approximation version of EMwith an MCMC procedure. ESAIM: Probability and Statistics. 2004;8:115–131.
Louis TA. Finding the observed information matrix when using the EM algorithm. Journal of the Royal Statistical Society Series B: Statistical Methodology. 1982;44(2):226–233.
Mattos TdB, Garay AM, Lachos VH. Likelihood-based inference for censored linear regression models with scale mixtures of skew-normal distributions. Journal of Applied Statistics. 2018;45(11):2039–2066.
Meilijson I. A fast improvement to the EM algorithm on its own terms. Journal of the Royal Statistical Society Series B: Statistical Methodology. 1989;51(1):127–138.
Meng XL, Rubin DB. Maximum likelihood estimation via the ECM algorithm: A general framework. Biometrika. 1993;80(2):267–278.
Nadarajah S, Kotz S. Skewed distributions generated by the normal kernel. Statistics & probability letters. 2003;65(3):269–277.
Nematollahi N, Farnoosh R, Rahnamaei Z. Location-scale mixture of skew-elliptical distributions: Looking at the robust modeling. Statistical Methodology. 2016;32:131–146.
Team RC. R: A language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria. 2021;.
Sahu SK, Dey DK, Branco MD. A new class of multivariate skew distributions with applications to Bayesian regression models. Canadian Journal of Statistics. 2003;31(2):129–150.
Vilca F, Balakrishnan N, Zeller CB. Multivariate skew-normal generalized hyperbolic distribution and its properties. Journal of Multivariate Analysis. 2014;128:73–85.
Wang J, Genton MG. The multivariate skew-slash distribution. Journal of Statistical Planning and Inference. 2006;136(1):209–220.
Wang J, Boyer J, Genton MG. A skew-symmetric representation of multivariate distributions. Statistica Sinica. 2004;p. 1259–1270.
Watagoda L, Don HSRA, Sanqui JAT. A cosine approximation to the skew normal distribution. In: Int. Math. Forum, vol. 14; 2019. p. 253–261.
Wei GC, Tanner MA. A Monte Carlo implementation of the EM algorithm and the poor man’s data augmentation algorithms. Journal of the American statistical Association. 1990;85(411):699–704.

Articles in Press, Accepted Manuscript
Available Online from 02 August 2026

  • Receive Date 15 February 2026
  • Revise Date 10 July 2026
  • Accept Date 25 July 2026