[论文解读] Möbius transformation and a Cauchy family on the sphere

Shogo Kato, Peter McCullagh|arXiv (Cornell University)|Oct 26, 2015

Morphological variations and asymmetry参考文献 4被引用 4

一句话总结

本文在单位球面 $S^d$ 上引入了球面柯西分布，扩展了圆周上的Wrapped Cauchy分布。证明了该分布族在球面上的Möbius变换下保持封闭，且参数变换由扩展的Möbius映射给出，从而可通过球极投影映射到多元 $t$-分布族，实现高效的推断。主要贡献包括闭式估计量、基于似然的高效推断，以及用于统计检验的枢轴分布。

ABSTRACT

We present some properties of a Cauchy family of distributions on the sphere, which is a spherical extension of the wrapped Cauchy family on the circle. The spherical Cauchy family is closed under the Möbius transformation on the sphere and there is a similar induced transformation on the parameter space. Stereographic projection transforms the the spherical Cauchy family into a multivariate $t$-family with a certain degree of freedom on Euclidean space. Many tractable properties of the spherical Cauchy are derived using the Möbius transformation and stereographic projection. A method of moments estimator and an asymptotically efficient estimator are expressed in closed form. The maximum likelihood estimation is also straightforward.

研究动机与目标

• 将圆周上的Wrapped Cauchy分布扩展至 $S^d$ 上的球面柯西分布族，用于方向数据分析。
• 建立球面柯西分布族在球面上Möbius变换下的封闭性，并推导参数空间上的诱导变换。
• 利用球极投影将球面柯西分布与欧氏空间上的多元 $t$-分布族联系起来，实现可 tractable 的推断。
• 推导闭式估计量——矩法估计量与渐近高效估计量——并实现简便的最大似然估计。
• 从分布特性与统计实用性两方面，比较球面柯西分布与von Mises–Fisher分布。

提出的方法

• 在 $S^d$ 上定义球面柯西密度为 $f(y;\mu,\rho) = \frac{\Gamma\{(d+1)/2\}}{2\pi^{(d+1)/2}} \left( \frac{1-\rho^2}{1+\rho^2 - 2\rho \mu^T y} \right)^d$，其中位置参数 $\mu \in S^d$，集中度参数 $\rho \in [0,1)$。
• 使用定义在 $S^d$ 上的Möbius变换：$\tilde{\cal M}_{R,\psi}(y) = R\left\{ \frac{1-\|\psi\|^2}{\|y+\psi\|^2}(y+\psi) + \psi \right\}$，其中 $\psi \in \mathbb{R}^{d+1} \setminus S^d$ 且 $R \in SO(d+1)$。
• 证明球面柯西分布族在Möbius变换下保持不变，且参数通过 $\mathbb{R}^{d+1}$ 上的扩展Möbius映射进行变换。
• 应用球极投影，将球面柯西分布映射为 $\mathbb{R}^d$ 上自由度为 $d$ 的多元 $t$-分布，实现似然等价性与可 tractable 推断。
• 推导矩法估计量与渐近高效估计量的闭式表达式，模拟研究显示后者在某些情形下表现更优。
• 证明费雪信息矩阵与最大似然估计均可直接计算，且在MLE处对数似然函数的黑塞矩阵为负定，确保局部极大性。

实验结果

研究问题

• RQ1 球面柯西分布族在 $S^d$ 上是否在Möbius变换下保持封闭？参数变换如何表征？
• RQ2 球极投影如何将球面柯西分布与欧氏空间上的多元 $t$-分布族关联起来？
• RQ3 能否为球面柯西分布推导出闭式估计量？其在有限样本中的表现如何？
• RQ4 球面柯西分布在变换不变性与尾部行为方面相比von Mises–Fisher分布有何统计优势？
• RQ5 球面柯西分布的最大似然估计是否计算上可行？是否达到渐近效率？

主要发现

• 球面柯西分布族在 $S^d$ 上的Möbius变换下保持封闭，参数变换由 $\mathbb{R}^{d+1}$ 上的扩展Möbius映射给出，支持基于不变性的推断。
• 球极投影将球面柯西分布转化为 $\mathbb{R}^d$ 上自由度为 $d$ 的多元 $t$-分布，似然函数等价，支持高效推断。
• 矩法估计量与渐近高效估计量均可表示为闭式，模拟结果表明后者在某些情形下表现更优。
• 最大似然估计计算直接，对数似然函数在MLE处的黑塞矩阵为负定，确保局部极大性与参数收敛。
• 球面柯西密度在众数与反模式处分配更高质量，而在赤道区域质量更低，表现出更重的尾部。
• 费雪信息矩阵可显式计算，且似然函数等价于多元 $t$-分布，简化了推断与假设检验。