[论文解读] Maximum Likelihood Estimation for Nets of Conics

Stefan Dye, Kathlén Kohn+2 more|arXiv (Cornell University)|Nov 17, 2020

Topological and Geometric Data Analysis参考文献 3被引用 5

一句话总结

本文研究了在代数统计背景下，3×3对称矩阵的3维线性空间（即二次曲线网）的最大似然估计（MLE）。通过分析每张网的对偶曲面（ℙ⁵中的有理曲面），作者计算了Wall（1977）分类的13种正则网的所有类型的最大似然次数（MLD）和对偶最大似然次数（rmld），发现关系式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 并非普遍成立，尽管它对类型A的通用网成立，此时 mld = 7 且 rmld = 7。

ABSTRACT

We study the problem of maximum likelihood estimation for $3$-dimensional linear spaces of $3 imes 3$ symmetric matrices from the point of view of algebraic statistics where we view these nets of conics as linear concentration or linear covariance models of Gaussian distributions on $\mathbb{R}^3$. In particular, we study the reciprocal surfaces of nets of conics which are rational surfaces in $\mathbb{P}^5$. We show that the reciprocal surfaces are projections from the Veronese surface and determine their intersection with the polar nets. This geometry explains the maximum likelihood degrees of these linear models. We compute the reciprocal maximum likelihood degrees. This work is based on Wall's classification of nets of conics from 1977.

研究动机与目标

• 计算ℙ⁵中所有13种正则二次曲线网的最大似然次数（MLD）和对偶最大似然次数（rmld）。
• 将每张网的对偶曲面ℙℒ⁻¹表征为ℙ⁵中Veronese曲面的投影。
• 将ML基点集确定为对偶曲面与极性网ℙℒ⊥的交集，并将其与网中公共二次曲线零点的几何结构联系起来。
• 探究恒等式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 是否可从线性系推广至二次曲线网。
• 确立通用型A类网的 mld = 7 且 rmld = 7，并推测A类网在对偶最大似然次数方面具有通用性。

提出的方法

• 基于判别三次曲线，使用Wall（1977）的分类方法，将所有正则二次曲线网划分为13种几何类型（A–H，含对偶类型）。
• 将对偶曲面ℙℒ⁻¹定义为集合{A⁻¹ | A ∈ ℒ, det(A) ≠ 0}在ℙ⁵中的Zariski闭包，并证明其为Veronese曲面的投影。
• 利用迹配对定义极性网ℙℒ⊥，并将ML基点集计算为ℙℒ⁻¹ ∩ ℙℒ⊥。
• 应用代数几何工具：次数计算、奇点分析（如在秩1点处的光滑性），以及当对偶曲面在基点集处光滑时的公式 mld(ℒ) = deg(ℙℒ⁻¹) − deg(base locus)。
• 使用Macaulay2计算非通用类型的对偶最大似然次数，并对A类类型进行多次采样以确认 rmld = 7。
• 利用对偶最大似然次数在合同变换下的不变性，按轨道类型计算结果。

实验结果

研究问题

• RQ1 ℙ⁵中每张正则二次曲线网的最大似然次数（MLD）是多少？
• RQ2 二次曲线网的对偶曲面ℙℒ⁻¹与Veronese曲面在几何上如何关联？
• RQ3 恒等式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 是否对所有二次曲线网成立，还是仅对通用网成立？
• RQ4 ML基点集ℙℒ⁻¹ ∩ ℙℒ⊥的结构是什么？它与网中二次曲线的公共零点有何关联？
• RQ5 是否只有A类网（通用网）的对偶最大似然次数为7？该值是否为所有3维对称矩阵线性空间的通用对偶最大似然次数？

主要发现

• 每张正则二次曲线网的对偶曲面ℙℒ⁻¹均为ℙ⁵中Veronese曲面的投影，且该投影仅依赖于网的类型（模合同变换）。
• 对于A类网（通用网），最大似然次数为4，对偶最大似然次数为7，确认A类网在对偶最大似然次数方面具有通用性。
• 恒等式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 并非对所有网成立：对F*类（rmld = 1 < 2）和G*类（rmld = 0 < 1）严格不成立，但对A–E类和H类成立。
• A、B、D、E、F*、G*、H类的ML基点集为空集；其余类型（如B*类的{(0:0:1)}和D*类的{(1:0:0),(0:1:0)}）对应实或复点（含重数）。
• 对于G*类，ML基点集为一个秩1的重合点，MLD为0，通过2次交点与基点处的光滑性计算得出。
• 基于采样与[9, Table 1]的一致性，推测A类网的对偶最大似然次数恒为7，且该值高于其他类型（如B*类为6，D*类为5，E*类为4）。