[論文レビュー] Maximum likelihood for high-noise group orbit estimation and single-particle cryo-EM

Zhou Fan, Roy R. Lederman+3 more|arXiv (Cornell University)|Jul 2, 2021

Advanced Electron Microscopy Techniques and Applications参考文献 36被引用数 4

ひとこと要約

本稿は、単粒子低温電子線シャドウ・エレクトロニクス（cryo-EM）における高ノイズ下でのグループ軌道推定のための最尤推定フレームワークを開発し、代数的不変量によって分類されたフィッシャー情報固有値を用いて対数尤度関数の形状を分析する。理論的予測を低次元のタンパク質ポテンシャルマップにおけるノイズスケーリングの実証的検証とともに、SO(2)およびSO(3)モデルにおいて第三階モーメントが回転を除いて局所的に信号を同定可能であることを証明する。

ABSTRACT

Motivated by applications to single-particle cryo-electron microscopy (cryo-EM), we study several problems of function estimation in a high noise regime, where samples are observed after random rotation and possible linear projection of the function domain. We describe a stratification of the Fisher information eigenvalues according to transcendence degrees of graded pieces of the algebra of group invariants, and we relate critical points of the log-likelihood landscape to a sequence of moment optimization problems, extending previous results for a discrete rotation group without projections. We then compute the transcendence degrees and forms of these optimization problems for several examples of function estimation under $SO(2)$ and $SO(3)$ rotations, including a simplified model of cryo-EM as introduced by Bandeira, Blum-Smith, Kileel, Perry, Weed, and Wein. We affirmatively resolve conjectures that $3^ ext{rd}$-order moments are sufficient to locally identify a generic signal up to its rotational orbit in these examples. For low-dimensional approximations of the electric potential maps of two small protein molecules, we empirically verify that the noise-scalings of the Fisher information eigenvalues conform with our theoretical predictions over a range of SNR, in a model of $SO(3)$ rotations without projections.

研究の動機と目的

• 高ノイズ下でのグループ軌道回復における最尤推定の統計的・幾何的構造を理解すること。
• 統計的推定と代数幾何学を結びつけるために、群不変量の超越次数に基づいてフィッシャー情報行列を特徴づけること。
• SO(2)およびSO(3)回転群における局所的信号同定に第三階モーメントが十分であるという予想を解消すること。
• 低次元のタンパク質ポテンシャルマップにおける数値シミュレーションを通じて、フィッシャー情報固有値のノイズスケーリングに関する理論的予測を検証すること。

提案手法

• SO(2)およびSO(3)群の下での群不変量の代数の各次数付き部分の超越次数に応じて、フィッシャー情報固有値を階層化する。
• 対数尤度関数の臨界点を段階的なモーメント最適化問題に再定式化する。
• SO(2)およびSO(3)群の下での関数推定における超越次数およびモーメント形式を計算し、簡略化されたcryo-EMモデルを含む。
• 回転推定の解析的ツールとして、直交プロクラステス整合と連続的多次元参照アラインメントを用いる。
• 信号対ノイズ比の変動に応じて、フィッシャー情報固有値を数値的シミュレーションで評価し、理論的予測と比較する。
• フーリエ断層撮影定理およびトモグラフィック投影モデルを用いて、cryo-EM像と2次元投影からの3次元ボリューム再構成を関連付ける。

実験結果

リサーチクエスチョン

• RQ1 SO(2)およびSO(3)群軌道モデルにおいて、第三階モーメントは一般の信号をその回転軌道を除いて局所的に同定可能か？
• RQ2 フィッシャー情報行列の構造は、回転群の代数的不変量とどのように関係するか？
• RQ3 信号対ノイズ比が変化するcryo-EMに類似したモデルにおいて、フィッシャー情報固有値のノイズスケーリング挙動はいかなるものか？
• RQ4 低次元のタンパク質ポテンシャルマップにおける理論的予測と実証的観測の間には、フィッシャー情報固有値にどのような差異があるか？
• RQ5 高ノイズ・回転不変推定モデルにおいて、対数尤度関数の形状はモーメント最適化問題として特徴づけられるか？

主な発見

• 本稿は、SO(2)およびSO(3)モデルにおいて、第三階モーメントが一般の信号をその回転軌道を除いて局所的に同定可能であることを確認し、先行研究からの予想を解決する。
• フィッシャー情報固有値が、群不変量の超越次数から導かれた理論的予測と整合したノイズ依存性を示し、数値的シミュレーションでさまざまなSNR範囲で検証された。
• 群不変量の次数付き部分の超越次数に応じてフィッシャー情報固有値を階層化することで、回転不変モデルにおける推定の難易度を理解するための新しい代数的幾何的枠組みが得られた。
• 2つの小分子タンパク質の低次元近似において、フィッシャー情報固有値のノイズスケーリングが理論的予測と一致し、実用的状況下でのモデルの正確性を確認した。
• 対数尤度関数の臨界点が、離散群に限らない段階的なモーメント最適化問題の解に対応することが示された。
• 理論的および数値的結果は、第二階モーメント手法が失敗する高ノイズ領域においても、高階モーメントを用いた最尤推定の有効性を支持する。