[論文レビュー] A Toolbox for the Radial and Angular Marginalization of Bivariate Normal Distributions

Emily A. Cooper, Hany Farid|arXiv (Cornell University)|May 19, 2020

Spectroscopy and Chemometric Analyses被引用数 5

ひとこと要約

本論文は、一般の条件下—非ゼロ平均、非等方的分散、非対角の共分散—における2変量正規分布の半径および角度の周辺化を計算するための包括的なMATLABおよびRリソース、RAMBiNo Toolboxを提示する。極座標における解析的積分に基づく閉形式および数値解法を提供し、半径および角度の周辺分布を効率的に計算可能である。

ABSTRACT

Bivariate normal distributions are often used to describe the joint probability density of a pair of random variables. These distributions arise across many domains, from telecommunications, to meteorology, ballistics, and computational neuroscience. In these applications, it is often useful to radially and angularly marginalize (i.e.,~under a polar transformation) the joint probability distribution relative to the coordinate system's origin. This marginalization is trivial for a zero-mean, isotropic distribution, but is non-trivial for the most general case of a non-zero-mean, anisotropic distribution with a non-diagonal covariance matrix. Across domains, a range of solutions with varying degrees of generality have been derived. Here, we provide a concise summary of analytic solutions for the polar marginalization of bivariate normal distributions. This report accompanies a Matlab (Mathworks, Inc.) and R toolbox that provides closed-form and numeric implementations for the marginalizations described herein.

研究の動機と目的

• 非ゼロ平均、非等方的分散、非対角共分散を含む一般の2変量正規分布の極座標における周辺化に関する散在する解析的解を統合し、統一的かつアクセス可能なリソースとして提供すること。
• 一般の2変量正規分布における半径および角度の周辺分布を計算するための中央集権的かつ十分に文書化された実装が不足している問題を解決すること。
• 神経科学、気象学、通信工学など、多様な分野の研究者に対して、閉形式および数値実装の両方を提供すること。
• 非ゼロ平均、非等方的性、非対角共分散という非自明な条件下でも、角度および半径の周辺化を効率的かつ正確に計算できることを実現すること。
• コード例およびすべてのケースの導出を含む完全な文書化済みツールボックスを公開することで、再現可能性の高い研究を支援すること。

提案手法

• 修正ベッセル関数を用いた解析的解法を導出し、特に非ゼロ平均・等方的ケースではライシアン分布が得られることを示す。
• 座標回転と変換技術を用いて共分散行列を対角化し、一般の非等方的ケースを標準形に還元する。
• 非ゼロ平均および非対角共分散ケースに対応するため、第一種修正ベッセル関数を用いた級数展開を採用する。
• 離散的な角度または半径のサンプリングを用いた円形経路に沿った数値積分により、数値的周辺化を実装する。
• パラメトリックアプローチを用い、半径および角度の点における同時密度を計算した後、円形経路に沿って統合して周辺密度を求める。
• MATLABおよびR用のモジュラーなコードを提供し、解析的および数値的周辺化の両方の関数を実装。また、切り捨てられた級数近似に対応する。

実験結果

リサーチクエスチョン

• RQ1 非ゼロ平均および非対角共分散を持つ2変量正規分布の半径および角度の周辺化に対する閉形式解は何か？
• RQ2 一般の2変量正規分布に対して、解析的および数値的手法を用いて半径の周辺分布をどのように効率的に計算できるか？
• RQ3 非等方的性およびオフセット平均が存在する状況において、角度の周辺分布を最も正確かつ計算効率的に近似する方法は何か？
• RQ4 通信工学、気象学など、異なる分野から得られた既存の解法は、数学的定式化および計算可能性の観点からどのように比較できるか？
• RQ5 科学的分野のあらゆる分野で極座標周辺化の計算を標準化・簡素化できる統一されたオープンソースツールボックスを開発可能か？

主な発見

• 本論文は、非ゼロ平均・等方的2変量正規分布に対して、ライシアン分布を含むすべての標準ケースにおける角度周辺化の閉形式解を導出し、実装した。
• 非ゼロ平均、非等方的、非対角共分散ケースでは、修正ベッセル関数を用いた級数展開に依拠し、数値的安定性と精度を確保するため100項まで切り捨てられている。
• 数値実装は、円形経路のサンプリングと同時密度の数値的積分により高い精度を達成しており、結果は単位積分を満たすように正規化されている。
• ツールボックスは、平均（ゼロ／非ゼロ）、分散（等方的／非等方的）、共分散（対角／非対角）のすべての組み合わせをカバーし、2変量正規分布の完全な一般性を網羅している。
• 解析的解法は、現代の科学的計算環境で利用可能な最適化されたベッセル関数評価を活用することで、計算効率的かつ数値的に安定していることが示された。
• ツールボックスはGitHubに公開されており、すべてのケースのコード例が含まれており、研究応用における即時利用と再現可能性を可能にしている。