대칭성 교란을 통한 분포 생성 및 다변량 스키우스 t 분포에 관한 연구

초록

본 논문은 약한 형태의 다변량 대칭성을 만족하는 밀도함수를 교란시켜 비대칭 밀도들의 전체 집합을 생성하는 일반적인 절차를 제시한다. 이 방법은 기존의 스키우스 정규분포와 관련된 여러 제안을 포괄하며, 특히 스키우스 타원형 밀도와 다변량 스키우스 t 분포에 초점을 맞춘다. 논문은 해당 분포들의 정의, 성질, 그리고 최대우도 추정 방법을 제시하고, 수치 예시를 통해 실용성을 검증한다.

상세 분석

논문은 먼저 “약한 대칭성”(weak symmetry)이라는 개념을 정의한다. 이는 원점에 대한 중앙대칭이 아니라, 어떤 고정된 방향벡터 α에 대해 f( x )=f( −x )·g(αᵀx) 형태로 표현될 수 있는 조건을 의미한다. 이러한 구조를 갖는 기본 밀도 f₀(x) 에 대해, 교란 함수 w(·)를 도입하여 새로운 밀도 f(x)=2 f₀(x) W(αᵀx) 를 만든다. 여기서 W 는 누적분포함수(CDF)이며, α는 스키우스 방향벡터이다. 이 일반식은 Azzalini와 Dalla Valle(1996)의 스키우스 정규분포를 포함하고, 이후 제안된 다양한 스키우스 타원형 모델들의 공통된 틀을 제공한다.

특히 타원형 기본밀도 f₀(x)=c |Σ|^{-1/2} g\big((x‑μ)ᵀΣ^{-1}(x‑μ)\big) 에 대해, g가 적절히 선택되면 스키우스 타원형 분포가 도출된다. 논문은 이때 α와 Σ가 독립적으로 파라미터화될 수 있음을 강조하며, 스키우스 방향과 스케일 구조가 분리되어 해석이 용이함을 보여준다.

다변량 스키우스 t 분포는 g(u)=\big(1+u/ν\big)^{-(ν+p)/2} 형태의 t-밀도에서 출발한다. 여기서 ν는 자유도, p는 차원이다. 교란 함수 W 를 표준 정규분포의 CDF인 Φ(·)로 두면, 최종 밀도는

f(x)=2 t_p(x;μ,Σ,ν) Φ\big(αᵀΣ^{-1/2}(x‑μ) \sqrt{(ν+p)/(ν+δ(x))}\big)

의 형태를 갖는다. δ(x)=(x‑μ)ᵀΣ^{-1}(x‑μ)이며, 이는 기존 스키우스 t 모델과 동일하지만, 자유도와 스케일 매개변수가 명시적으로 분리되어 추정이 직관적이다.

추정 측면에서 논문은 EM 알고리즘을 활용한 최대우도 추정 절차를 제시한다. 잠재 변수 Z 를 도입해 t-분포를 정규-가우시안 혼합 형태로 표현하고, E‑단계에서 조건부 기대값 E