[论文解读] Existence and non-existence of the non-central Wishart distributions

Gérard Letac, Hélène Massam|arXiv (Cornell University)|Aug 14, 2011

Point processes and geometric inequalities被引用 6

一句话总结

本文通过分析形如 $(\det s)^{-p} \exp(\mathrm{tr}(s^{-1}w))$ 的拉普拉斯变换，解决了正半定矩阵锥上非中心 Wishart 分布的存在性条件，其中 $w$ 为秩 $k$ 的对角矩阵。证明了 $\mathcal{NCW}(2p,w,\Sigma)$ 存在当且仅当 $p \in \Lambda_d$ 且当 $p \in \{1/2, 1, \dots, (d-2)/2\}$ 时 $k \leq 2p$；而当 $p \geq (d-1)/2$ 时，对所有 $w$ 均存在，且显式计算了 $m(d-1,d,d)$，并证明了 $m(d-2,d,d)$ 和 $m(d-2,d-1,d)$ 不存在。结果解决了 Mayerhofer 长期以来的猜想，并弥补了该领域先前文献中的一个漏洞。

ABSTRACT

The problem considered in this paper is to find when the non-central Wishart distribution, defined on the cone $\bar{\mathcal{P}_d}$ of semi positive definite matrices of order $d$ and with a real valued shape parameter, exists. We reduce this problem to the problem of existence of the measures $m(n,k,d)$ defined on $\bar{\mathcal{P}_d}$ and with Laplace transform $(\det s)^{-n/2}\exp r(s^{-1}w)$ where $n$ is an integer and where $w=\mathrm{diag}(0,...,0,1,...,1)$ has order $d$ and rank $k.$ We compute $m(d-1,d,d)$ and we show that neither $m(d-2,d,d)$ nor $m(d-2,d-1,d)$ exist. This proves a conjecture of E. Mayerhofer.

研究动机与目标

• 确定非中心 Wishart 分布 $\mathcal{NCW}(2p,w,\Sigma)$ 在正半定矩阵锥上存在的精确条件。
• 通过证明当 $p \in \{1/2, 1, \dots, (d-2)/2\}$ 时 $\mathcal{NCW}(2p,w,\Sigma)$ 仅在 $k \leq 2p$ 时存在，从而解决文献中长期存在的漏洞，纠正先前工作中错误的断言。
• 显式计算测度 $m(d-1,d,d)$ 并确立 $m(d-2,d,d)$ 和 $m(d-2,d-1,d)$ 的非存在性，从而完成此类分布的分类。
• 通过 zonal 多项式方法与线性代数，为非中心 Wishart 分布的存在性提供严格的理论基础，尤其在临界情形 $p = (d-1)/2$ 下。

提出的方法

• 本文分析了 $s \in \mathcal{P}_d$ 的拉普拉斯变换 $\mathbb{E}[e^{-\mathrm{tr}(sX)}] = (\det s)^{-p} \exp(\mathrm{tr}(s^{-1}w))$，其中 $w = \mathrm{diag}(0,\dots,0,1,\dots,1)$ 秩为 $k$。
• 定义了在 $\overline{\mathcal{P}_d}$ 上的测度 $m(n,k,d)$，其具有上述拉普拉斯变换，从而将存在性问题转化为对这些测度的研究。
• 作者利用 zonal 多项式展开及欧几里得 Jordan 代数的性质，分析了这些测度的结构，尤其关注低维情形。
• 通过极限论证表明 $\mathcal{NCW}(d-1,w,\Sigma)$ 是 $\mathcal{NCW}(2p,w,\Sigma)$ 当 $p \downarrow (d-1)/2$ 时的极限，从而在临界情形下确立了存在性。
• 一个关键技术工具是使用变换 $G(c_2,\dots,c_d)$ 将 $\mathbb{R}^{d-1}$ 中的向量映射到对称矩阵，用于推导秩约束，并通过反证法证明非存在性。
• 对 $m(d-2,d,d)$ 和 $m(d-2,d-1,d)$ 的非存在性证明依赖于：该测度在 $G$ 下的像无法具有满秩 $d-1$，与必须存在绝对连续部分的要求矛盾。

实验结果

研究问题

• RQ1 对于给定的 $p > 0$ 和秩为 $k$ 的 $w \in \overline{\mathcal{P}_d}$，非中心 Wishart 分布 $\mathcal{NCW}(2p,w,\Sigma)$ 在什么条件下存在？
• RQ2 当 $p \in \{1/2, 1, \dots, (d-2)/2\}$ 时，是否存在必须满足的秩约束 $k \leq 2p$ 以保证 $\mathcal{NCW}(2p,w,\Sigma)$ 的存在性？
• RQ3 当 $p \geq (d-1)/2$ 时，包括临界情形 $p = (d-1)/2$，$\mathcal{NCW}(2p,w,\Sigma)$ 是否对所有 $w$ 都存在？
• RQ4 测度 $m(d-1,d,d)$ 的显式形式是什么？为何 $m(d-2,d,d)$ 和 $m(d-2,d-1,d)$ 不存在？
• RQ5 Mayerhofer 关于 $\mathrm{rank}(w) \leq 2p$ 是存在性的必要条件的猜想能否被严格证明？

主要发现

• 测度 $m(d-1,d,d)$ 存在且已显式计算，为 $p = (d-1)/2$ 且 $w$ 满秩时的非中心 Wishart 分布提供了完整描述。
• $m(d-2,d,d)$ 和 $m(d-2,d-1,d)$ 不存在，证明了当 $p \in \{1/2, 1, \dots, (d-2)/2\}$ 时，秩约束 $k \leq 2p$ 是必要的。
• 当 $p \in \{1/2, 1, \dots, (d-2)/2\}$ 时，非中心 Wishart 分布 $\mathcal{NCW}(2p,w,\Sigma)$ 存在当且仅当 $\mathrm{rank}(w) \leq 2p$。
• 当 $p \geq (d-1)/2$ 时，$\mathcal{NCW}(2p,w,\Sigma)$ 对所有 $w \in \overline{\mathcal{P}_d}$ 均存在，包括临界情形 $p = (d-1)/2$ 作为极限。
• 本文证实了 Mayerhofer 的猜想，即 $\mathrm{rank}(w) \leq 2p$ 是存在性的必要条件，解决了长期悬而未决的问题。
• 该证明方法基于线性代数与测度变换，提供了一种新颖且自包含的途径，与此前使用的随机过程方法截然不同。