Robust mixture modeling using multivariate skew t distributions

被引:179
作者
Lin, Tsung-I [1 ,2 ]
机构
[1] Natl Chung Hsing Univ, Dept Appl Math, Taichung 402, Taiwan
[2] Natl Chung Hsing Univ, Inst Stat, Taichung 402, Taiwan
关键词
MCEM-type algorithms; MSN; MST; Multivariate truncated normal; Multivariate truncated t; Outliers; MAXIMUM-LIKELIHOOD-ESTIMATION; UNKNOWN NUMBER; EM ALGORITHM; INFERENCE; ECM;
D O I
10.1007/s11222-009-9128-9
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
This paper presents a robust mixture modeling framework using the multivariate skew t distributions, an extension of the multivariate Student's t family with additional shape parameters to regulate skewness. The proposed model results in a very complicated likelihood. Two variants of Monte Carlo EM algorithms are developed to carry out maximum likelihood estimation of mixture parameters. In addition, we offer a general information-based method for obtaining the asymptotic covariance matrix of maximum likelihood estimates. Some practical issues including the selection of starting values as well as the stopping criterion are also discussed. The proposed methodology is applied to a subset of the Australian Institute of Sport data for illustration.
引用
收藏
页码:343 / 356
页数:14
相关论文
共 41 条
[1]  
[Anonymous], 2000, Sankhya Ser. A, DOI DOI 10.2307/25051289
[2]  
[Anonymous], 2008, EM ALGORITHM EXTENSI
[3]   Bayesian inference for skew-normal linear mixed models [J].
Arellano-Valle, R. B. ;
Bolfarine, H. ;
Lachos, V. H. .
JOURNAL OF APPLIED STATISTICS, 2007, 34 (06) :663-682
[4]   Statistical applications of the multivariate skew normal distribution [J].
Azzalini, A ;
Capitanio, A .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY, 1999, 61 :579-602
[5]   The multivariate skew-normal distribution [J].
Azzalini, A ;
DallaValle, A .
BIOMETRIKA, 1996, 83 (04) :715-726
[6]   The skew-normal distribution and related multivariate families [J].
Azzalini, A .
SCANDINAVIAN JOURNAL OF STATISTICS, 2005, 32 (02) :159-188
[7]   Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t-distribution [J].
Azzalini, A ;
Capitanio, A .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY, 2003, 65 :367-389
[8]  
Basford KE, 1997, COMPUTATION STAT, V12, P1
[9]   Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm [J].
Booth, JG ;
Hobert, JP .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY, 1999, 61 :265-285
[10]  
Cook RD., 1994, INTRO REGRESSION GRA