[논문 리뷰] Moments Calculation For the Doubly Truncated Multivariate Normal Density

B.G. Manjunath, Stefan Wilhelm|arXiv (Cornell University)|2012. 06. 23.

Statistical Distribution Estimation and Applications참고 문헌 16인용 수 28

한 줄 요약

이 논문은 모멘트 생성 함수 접근법을 사용하여 임의의 직사각형 이중 절단 하에서 다변량 정규분포의 잘라낸 평균과 분산에 대한 명시적 공식을 유도한다. 주요 기여는 모멘트에 대한 폐쇄형 해와 절단 이후 정밀도 행렬에서의 불변성 성질이며, 통계 모델링 및 그래픽 모델링에 실용적으로 활용할 수 있도록 R에 구현되어 있다.

ABSTRACT

In the present article we derive an explicit expression for the trun- cated mean and variance for the multivariate normal distribution with ar- bitrary rectangular double truncation. We use the moment generating ap- proach of Tallis (1961) and extend it to general μ, Σ and all combinations of truncation. As part of the solution we also give a formula for the bivari- ate marginal density of truncated multinormal variates. We also prove an invariance property of some elements of the inverse covariance after trunca- tion. Computer algorithms for computing the truncated mean, variance and the bivariate marginal probabilities for doubly truncated multivariate normal variates have been written in R and are presented along with three examples.

연구 동기 및 목표

• 임의의 직사각형 이중 절단 하에서 다변량 정규분포의 1차 및 2차 모멘트를 위한 일반적 해를 개발하는 것.
• Tallis(1961)의 모멘트 생성 함수 방법을 일반적인 평균 벡터와 공분산 행렬을 갖는 경우에 대해 전체 이중 절단을 다룰 수 있도록 확장하는 것.
• 절단된 다변량 정규 변수의 이변량 근접 밀도에 대한 공식을 제공하는 것.
• 특히 절단되지 않은 변수에 대해 절단 후 정밀도 행렬 원소에 대한 불변성 성질을 증명하는 것.
• 금융 및 그래픽 모델링 분야의 응용을 지원하기 위해 R에 절단된 모멘트와 근접 확률을 위한 계산 알고리즘을 구현하는 것.

제안 방법

• Tallis(1961)의 모멘트 생성 함수 접근법을 변형하여 일반적인 이중 절단 하에서 절단된 평균과 분산에 대한 명시적 표현을 도출하는 것.
• 분할된 행렬 역행렬과 Johnson-Kotz 공식을 사용하여 절단된 공분산 및 정밀도 행렬을 표현하는 것.
• 해결 프레임워크의 일부로 절단된 다변량 정규 변수의 이변량 근접 밀도 함수를 도출하는 것.
• 절단되지 않은 변수에 대해 정밀도 행렬의 비대각 원소가 절단 후에도 그대로 유지됨을 증명하는 것.
• R에 절단된 모멘트, 근접 확률, 정밀도 행렬 원소를 계산하기 위한 수치 알고리즘을 구현하는 것.
• Whittaker(1990)의 빛나는 그래픽 모델을 포함한 세 가지 예시를 통해 결과를 검증하는 것.

실험 결과

연구 질문

• RQ1 임의의 직사각형 이중 절단 하에서 다변량 정규분포의 1차 및 2차 모멘트에 대한 폐쇄형 표현은 무엇인가?
• RQ2 이중 절단 다변량 정규 변수의 이변량 근접 밀도를 명시적으로 유도할 수 있는가?
• RQ3 이중 절단 후 어떤 정밀도 행렬 원소가 불변성을 유지하며, 그 조건은 무엇인가?
• RQ4 모멘트 생성 함수 방법을 어떻게 확장하여 이중 절단 경우의 일반적인 평균 벡터와 공분산 행렬을 다룰 수 있는가?
• RQ5 불변성 성질이 선택 후 그래픽 모델에서 부분 상관계수 추정에 미치는 계산적 영향은 무엇인가?

주요 결과

• 임의의 직사각형 이중 절단 하에서 다변량 정규분포의 잘라낸 평균과 분산에 대한 명시적 폐쇄형 표현이 도출되었다.
• 절단된 다변량 정규 변수의 이변량 근접 밀도가 해의 일부로 도출되어 공동 확률 계산이 가능해졌다.
• 불변성 성질이 증명되었으며, 절단되지 않은 변수에 대해 정밀도 행렬의 비대각 원소는 절단 후에도 그대로 유지된다.
• 절단된 변수의 정밀도 행렬 대각 원소는 갱신되며, 이는 선택으로 인한 부분 분산의 변화를 반영한다.
• 유도된 공식은 R 패키지 `tmvtnorm`에 구현되어 있어 절단된 모멘트와 근접 확률의 효율적 계산을 가능하게 한다.
• 이 방법은 금융 모델링(예: 자동 콜러블 파생상품) 및 그래픽 모델링, 특히 조건부 독립성 구조에서의 응용을 지원한다.