[논문 리뷰] Erdelyi-Kober Fractional Integral Operators from a Statistical Perspective -II

A. M. Mathai, H. J. Haubold|arXiv (Cornell University)|2013. 03. 16.

Numerical methods in inverse problems참고 문헌 5인용 수 3

한 줄 요약

이 논문은 에르데리-코버 분수적 적분 연산자의 통계적 해석을 행렬변량 사례로 확장하여, 독립적인 행렬변량 랜덤 변수들의 곱과 비율의 확률밀도함수를 유도한다. 첫 번째 및 두 번째 종류의 행렬변량 코버 연산자를 도입하고, 경로 모델과 행렬변량 하이퍼기하급수를 통해 이를 일반화하며, 사이고 연산자와 다변량 통계에서의 베타 분포와의 연결 고리를 설정한다.

ABSTRACT

In this article we examine the densities of a product and a ratio of two real positive definite matrix-variate random variables $X_1$ and $X_2$, which are statistically independently distributed, and we consider the density of the product $U_1=X_2^{1\over2}X_1X_2^{1\over2}$ as well as the density of the ratio $U_2=X_2^{1\over2}X_1^{-1}X_2^{1\over2}$. We define matrix-variate Kober fractional integral operators of the first and second kinds from a statistical perspective, making use of the derivation in the predecessor of this paper for the scalar variable case, by deriving the densities of products and ratios where one variable has a matrix-variate type-1 beta density and the other variable has an arbitrary density. Various types of generalizations are considered, by using pathway models, by appending matrix variate hypergeometric series etc. During this process matrix-variate Saigo operator and other operators are also defined and properties studied.

연구 동기 및 목표

• 독립적인 랜덤 행렬의 행렬변량 곱과 비율의 확률밀도함수를 유도하기 위해.
• 통계적 관점에서 첫 번째 및 두 번째 종류의 행렬변량 코버 분수적 적분 연산자를 정의하고 연구하기 위해.
• 경로 모델과 행렬변량 하이퍼기하급수를 사용하여 이러한 연산자를 일반화하기 위해.
• 유도된 연산자와 기존의 다변량 분포, 예를 들어 제1형 베타 분포와 감마 분포 사이의 연결 고리를 설정하기 위해.
• 극한 형태와 행렬 매개변수화, 예를 들어 매개변수 $ q \to 1_{-} $인 경로 모델을 포함하도록 프레임워크를 확장하기 위해.

제안 방법

• 행렬변량 제1형 베타 분포를 가진 $ X_1 $과 임의의 밀도를 가진 $ X_2 $를 고려하여, $ (X_1, X_2) $에서 $ (U_1, X_2) $로의 자코비안 변환을 통해 $ U_1 = X_2^{1/2} X_1 X_2^{1/2} $의 밀도를 유도한다.
• 동일한 변환을 통해 $ U_2 = X_2^{1/2} X_1^{-1} X_2^{1/2} $의 밀도를 유도함으로써 첫 번째 종류의 코버 연산자를 도출한다.
• 두 번째 종류의 코버 연산자를 $ K_X^{\zeta,\alpha}f(X) = \frac{|X|^\zeta}{\Gamma_p(\alpha)} \int_{T>X} |T-X|^{\alpha - \frac{p+1}{2}} |T|^{-\zeta - \alpha} f(T) dT $로 정의하며, $ \Re(\alpha) > \frac{p-1}{2} $를 조건으로 한다.
• 행렬변량 밀도의 경로 모델을 도입하여 제1형 베타 분포를 $ f_1(X_1) \propto |X_1|^{\gamma - \frac{p+1}{2}} |I - a(1-q)X_1|^{\frac{\eta}{1-q}} $로 대체함으로써, $ q \to 1_{-} $일 때 지수형 밀도로의 연속적 전이를 가능하게 한다.
• 경로 모델 하에서 $ U $의 주변 밀도가 제1종 일반화 코버 연산자로 이어지며, 이는 $ |U|^{-\gamma - (\frac{\eta}{1-q} + \frac{p+1}{2})} \int_{U > a(1-q)V} |U - a(1-q)V|^{\frac{\eta}{1-q}} |V|^\gamma f(V) dV $의 형태를 가진다.
• 극한 형태 $ q \to 1_{-} $를 유도하여, 지수 함수 $ \exp(-a\eta \cdot \mathrm{tr}(V^{1/2} U^{-1} V^{1/2})) $를 포함하는 밀도로 수렴함을 보이며, 감마 분포 및 일반화된 감마 분포와의 연결 고리를 확인한다.

실험 결과

연구 질문

• RQ1 스칼라에서 행렬변량 사례로의 코버 분수적 적분 연산자를 통계적으로 의미 있는 방식으로 어떻게 일반화할 수 있는가?
• RQ2 행렬변량 코버 연산자 제1종의 통계적 해석은 독립적인 랜덤 행렬의 비율 분포 측면에서 어떻게 해석될 수 있는가?
• RQ3 경로 모델은 어떻게 행렬변량 코버 연산자를 확장하며, 경로 매개변수 $ q \to 1_{-} $일 때의 극한 행동은 어떠한가?
• RQ4 행렬변량 하이퍼기하급수와 사이고 연산자는 이 일반화된 프레임워크에서 어떤 역할을 하는가?
• RQ5 유도된 연산자들은 기존의 다변량 분포, 예를 들어 제1형 베타 분포와 행렬변량 감마 분포와 어떻게 관련되어 있는가?

주요 결과

• $ X_1 \sim \text{Beta}_p(\zeta + \frac{p+1}{2}, \alpha) $이고 $ X_2 $가 임의의 밀도를 가지는 경우, $ U_1 = X_2^{1/2} X_1 X_2^{1/2} $의 밀도는 두 번째 종류의 코버 연산자로 이어진다.
• $ U_2 = X_2^{1/2} X_1^{-1} X_2^{1/2} $의 주변 밀도는 제1종 코버 연산자 $ I_X^{\zeta,\alpha}f(X) $에 비례하며, $ \Gamma_p(\zeta)/\Gamma_p(\zeta + \alpha) \cdot g_r(U) = I_X^{\zeta,\alpha}f(X) $로 표현된다.
• 제1종 코버 연산자의 경로 확장은 제1형 베타 밀도를 $ |I - a(1-q)X_1|^{\frac{\eta}{1-q}} $를 포함하는 일반화된 형태로 대체함으로써 도출되며, 더 넓은 종류의 연산자 클래스를 제공한다.
• $ q \to 1_{-} $일 때 경로 모델은 지수 함수 $ \exp(-a\eta \cdot \mathrm{tr}(V^{1/2} U^{-1} V^{1/2})) $를 포함하는 밀도로 수렴하며, 이는 일반화된 감마 분포와의 연결 고리를 형성한다.
• 경로 모델의 극한 형태는 $ \lim_{q \to 1_{-}} g_p(U) = \frac{(a\eta)^{p\gamma}}{\Gamma_p(\gamma)} |U|^{-\gamma - \frac{p+1}{2}} \int_{V>O} |V|^\gamma \exp(-a\eta \cdot \mathrm{tr}(V^{1/2} U^{-1} V^{1/2})) f(V) dV $로 표현되며, 감마 분포 및 역감마 분포와의 연결 고리를 확인한다.
• 이 프레임워크는 경로 모델에 정의된 행렬 $ A > O $를 포함할 수 있으며, 스칼라 매개변수 $ a $를 $ A > O $로 일반화할 수 있다. 이에 따라 유도된 밀도는 $ g_A(U) \propto |U|^{-\gamma - (\frac{\eta}{1-q} + \frac{p+1}{2})} \int_{U > (1-q)V^{1/2} A V^{1/2}} |U - (1-q)V^{1/2} A V^{1/2}|^{\frac{\eta}{1-q}} |V|^\gamma f(V) dV $의 형태로 표현된다.