다차원 정규분포 계산을 위한 교정된 추상 튜브와 사전 교정 다면체 이론

초록

본 논문은 선형 부등식으로 정의되는 폐합 볼록 다면체 K에 대해, 외부 방향의 사전 교정(lexicographic perturbation)을 적용한 추상 튜브(abstract tube) 이론을 정교화한다. 선형계획법을 이용해 교정된 부등식 시스템으로부터 추상 튜브의 복합체 F(0⁺)를 효율적으로 구성하는 알고리즘을 제시하고, 이를 Miwa‑Hayter‑Kuriki(2003)의 재귀적 적분 기법과 결합해 다차원 정규 확률 P(K)를 정확히 계산한다. 또한 학생화 범위 통계량의 분포 함수 계산을 사례로 들어 실험 결과를 보여준다.

상세 분석

이 논문은 Naiman과 Wynn(1997)이 제시한 “추상 튜브” 개념을 다면체의 부등식 시스템이 일반 위치(general position)에 있지 않을 때 발생하는 문제점을 해결하기 위해 확장한다. 일반 위치란 어떤 J⊆{1,…,m}에 대해 교차된 초평면 ∩_{i∈J}∂H_i가 차원 max{n+1−|J|,0}인 아핀 부분공간을 포함하지 않는 경우를 말한다. 기존 정의는 약했으며, 저자들은 보다 강력한 정의를 채택해 교정된 시스템이 언제든 일반 위치에 놓이도록 보장한다.

핵심 교정은 (4)식에 나타난 사전 교정(b(ε)=b+ε·(1,2,…,m)ᵀ)이다. ε>0를 충분히 작게 잡으면 각 초평면 ∂H_i(ε)는 원래 위치에서 미세하게 외부로 이동하고, 이때 발생하는 다면체 K(ε)는 일반 위치에 있다. Lemma 2.1은 ε가 충분히 작을 경우(1) 모든 초평면 집합이 일반 위치에 놓이며, (2) 복합체 F(ε)가 ε에 의존하지 않음을 증명한다. 즉, F(0⁺)=F(ε)는 ε→0⁺ 한계에서 고정된 단순 복합체가 된다.

다음으로 저자들은 F(0⁺)가 원래 복합체 F보다 작거나 같으며(|F|≥|F(0⁺)|)임을 Lemma 2.2와 Theorem 2.1을 통해 보인다. 특히, F(0⁺)에 포함된 모든 J에 대해 {a_i | i∈J}는 선형 독립이며, |J|는 행렬 A의 랭크를 초과하지 않는다. 따라서 ∩_{i∈J}H_i^c는 n‑차원 단순 원뿔(simple cone) 혹은 단순 원뿔과 선형 부분공간의 직합 형태가 된다. 이는 Miwa‑Hayter‑Kuriki(2003)의 재귀 적분 알고리즘이 직접 적용 가능함을 의미한다.

추상 튜브의 실제 구성은 선형계획법(LP)으로 변환된다. Problem 3.1은 주어진 J에 대해 시스템 a_iᵀx = b_i(ε) (i∈J), a_iᵀx ≤ b_i(ε) (i∉J)의 실현 가능성을 검사한다. 이를 표준 형태의 tableau로 만든 뒤, Irie(1973)의 알고리즘 3.1을 적용해 feasibility를 판정한다. 이 절차는 모든 J⊆{1,…,m}에 대해 반복되며, 결과적으로 F(0⁺)를 완전하게 식별한다.

마지막으로, 식 (6)의 기대값을 취하면
1−P(K)=∑{J∈F(0⁺)} (−1)^{|J|−1} P(∩{i∈J}H_i^c)
가 된다. 각 교차 원뿔에 대한 정규 확률은 앞서 언급한 재귀 적분법으로 효율적으로 계산 가능하므로, 전체 P(K)를 고차원에서도 실용적으로 구할 수 있다. 논문은 3‑차원 피라미드와 중복 부등식 예시, 그리고 Tukey의 학생화 범위 통계량에 대한 계산을 통해 알고리즘의 정확도와 계산량 감소 효과를 실증한다.