[論文レビュー] 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）を計算すること。
• 各ネットの逆曲面 ℙℒ⁻¹ を、ℙ⁵内に埋め込まれたヴェロネーゼ曲面の射影として特徴付けること。
• MLベース点集合を逆曲面と極的ネット ℙℒ⊥ の交点として定義し、それが共通する円錐曲線の零点の幾何的性質とどのように関係するかを明らかにすること。
• pencil（線型空間）からネットへの一般化において、恒等式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 がすべてのネットに成り立つかどうかを調査すること。
• 一般ネット（タイプA）が rmld = 7 かつ mld = 7 を満たすことを確立し、逆最大尤度次数が7であるタイプAネットが、逆最大尤度次数に関して一般であると仮説を立てる。

提案手法

• Wall（1977）の分類に従い、判別三次曲線に基づいて13種類の幾何的タイプ（A–H、およびその双対）に、すべての正則円錐曲線のネットを分類する。
• 逆曲面 ℙℒ⁻¹ を、ℙ⁵内における {A⁻¹ | A ∈ ℒ, det(A) ≠ 0} のZariski閉包として定義し、これがヴェロネーゼ曲面の射影であることを証明する。
• トレースペアリングを用いて極的ネット ℙℒ⊥ を定義し、MLベース点集合を ℙℒ⁻¹ ∩ ℙℒ⊥ として計算する。
• 代数幾何学的手法を適用：次数の計算、特異点解析（例：ランク1点での滑らかさ）、および逆曲面がベース点集合で滑らかである場合に mld(ℒ) = deg(ℙℒ⁻¹) − deg(base locus) が成り立つ公式を用いる。
• Macaulay2を用いて非一般タイプの逆最大尤度次数を計算し、タイプAの複数の例をサンプリングして rmld = 7 であることを確認する。
• 合同作用による不変性を活用し、各軌道タイプごとに結果を計算する。

実験結果

リサーチクエスチョン

• RQ1 ℙ⁵内のすべての正則円錐曲線のネットの最大尤度次数（MLD）は何か？
• RQ2 円錐曲線のネットの逆曲面 ℙℒ⁻¹ は、幾何学的にヴェロネーゼ曲面とどのように関係するか？
• RQ3 恒等式 rmld(ℒ) = deg(ℙℒ⁻¹) + mld(ℒ) − 1 はすべてのネットに成り立つか、それとも一般ネットにしか成り立たないのか？
• RQ4 MLベース点集合 ℙℒ⁻¹ ∩ ℙℒ⊥ の構造は何か？ また、それがネットに含まれる円錐曲線の共通零点とどのように関係するか？
• RQ5 タイプA（一般ネット）のネットが唯一 rmld = 7 を満たすものであり、かつすべての3次元対称行列線形空間に対して一般の逆最大尤度次数が7であると推測できるか？

主な発見

• すべての正則円錐曲線のネットの逆曲面 ℙℒ⁻¹ は、ℙ⁵内に埋め込まれたヴェロネーゼ曲面の射影であり、これはネットのタイプ（合同変換を除いて）にのみ依存する。
• タイプ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は2次交差とその点での滑らかさを用いて計算され、0である。
• タイプAの逆最大尤度次数は、サンプリングと[9, Table 1]との整合性に基づき、すべての同様のネットで7であると仮説を立てた。この値は他のタイプ（例：B*で6、D*で5、E*で4）よりも高い。