mnorm: An R Package for Calculation and Differentiation of Conditional Multivariate Normal Densities and Probabilities

あなたは天気を予測しようとしていると想像してください。しかし、単に雨か晴れかといったことではなく、気温、湿度、風速、気圧、さらには地元のリスの機嫌まで、十数もの変数を同時に操りこなしているのです。統計学の世界では、この混沌としたジャグリングのような行為は「多変量正規分布」と呼ばれます。これは、さまざまな事象がどのように連動して動くかを記述するための数学的な手法です。時には、それらは同期して動くこともあります（風と雨のように）。また、時にはそうでないこともあります。

さて、外を覗いてみて、風が時速20マイルで吹いているのが見えたとしましょう。すると突然、あなたのリスの機嫌に対する予測が変わり、湿度に対する予想も変化します。これは「条件付き確率」と呼ばれます。つまり、ある一部についてすでに知っている事実がある場合に、システムの他の部分がどのように動いているかを突き止めることです。

なぜこれほどまでに重要なのでしょうか？ それは、現実世界がこうした複雑に絡み合った網で満ちているからです。経済学者は、金利の変化がどのように株式市場に波及するかを理解するためにこれらを用いますし、学生の数学のテストの点数が歴史の点数をどのように予測するかを知るためにも用います。これらの予測を正確にするために、科学者たちは「微分」と呼ばれる作業を行う必要があります。微分とは、「丘の傾斜」を見つけることだと考えてください。もしあなたがキャンプを設営するのに最適な場所（最尤値）を見つけようとしているなら、どちらの方向が「上」で、どれくらい急な斜面なのかを知る必要があります。もし傾斜を完璧に計算できれば、真っ直ぐ頂上へと進むことができます。もし推測に頼るようなことがあれば、長い間、円を描くように彷徨うことになるでしょう。

ここで、新しいツールである「mnorm」パッケージが登場します。これは、これらの統計的な山の頂上を目指すための究極のGPSとなるR言語用のソフトウェアです。

---

mnorm パッケージ：統計的な山々へのスピーディなガイド

この論文では、これらの複雑な統計的な山々に対する高速なマルチツール・計算機として機能する新しいソフトウェアパッケージ、「mnorm」を紹介しています。mnormが登場する前、研究者たちは、処理が遅かったり、「条件付き」の状況（すでにいくつかの事実を知っている場合）を扱えなかったり、あるいは最善の答えを素早く見つけるために必要な「傾斜」（勾配）を計算できなかったりする他のツールを使わなければなりませんでした。

確率を計算することを、奇妙な形をした目に見えない雲の体積を測ろうとすることだと考えてみてください。その雲が2次元や3次元であれば、それは風船の体積を測るようなものです。しかし、もしその雲が10次元や100次元であったなら、それは追加の次元が存在する宇宙に存在する雲を測ろうとするようなものです。mnormパッケージが特別なのは、単に雲の大きさを測るだけでなく、風や温度、あるいは出発点を少し調整したときに、雲の形がどのように変化するかを正確に伝えることができる点にあります。

強力な武器：解析的勾配
mnormの最大の秘策は、「解析的勾配」を計算できる能力です。平たく言えば、これは、数値をわずかに動かして様子を見るという「推測」による方法ではなく、正確な数学公式を用いて曲線の「傾斜」を導き出せることを意味します。

著者らはシミュレーションを行い、このテストを行いました。彼らは、これらの正確な公式を使用することが、推測による方法よりも劇的に速いことを発見しました。

• 単純な密度計算（雲の大きさを測るだけ）において、mnormは高速でした。
• 勾配（傾斜）を計算する場合、その差は歴然でした。5次元のテストでは、傾斜を解析的に計算するのに約0.09秒しかかかりませんでしたが、「推測」による方法（数値微分）では、最低限必要な回数の推測を行うだけで少なくとも0.54秒はかかっていました。
• 10次元のような高次元では、解析的手法は3.47秒で傾斜を計算しましたが、推測による方法は少なくとも4.51秒を要しました。

この論文は、このスピードアップが、モデル構築のためにこれらの計算を数百万回繰り返す必要がある経済学者にとって、ゲームチェンジャーであることを示唆しています。それは、山を歩いて登るのか、ヘリコプターで移動するのかほどの違いがあります。

「与えられた」部分の扱い
最も有用な機能の一つは、mnormが「条件付き」データをどのように扱うかという点です。例えば、ある学生の最終成績を予測しようとしているとします。もし、その学生の中間テストの点数が分からなければ、あなたは推測するしかありません。しかし、もし中間テストの点数（「与えられた」部分）を知っているならば、予測は変わります。

mnormは、その特定の中間テストの点数が「与えられた」状態での確率と傾斜を計算できます。論文によれば、これらの傾斜を計算するのが非常に速いのは、条件付きの部分の数学が、他の部分の数学の単純な再構成に過ぎないからです。5つの変数があり、そのうち2つが「与えられた」ものとして扱われるテストにおいて、パッケージはほぼ瞬時に傾斜を計算しましたが、推測による方法ははるかに時間がかかったはずです。

小規模から大規模まで
このパッケージは、以下の2つの領域でも優れた性能を発揮します。

1. 小さな雲（2次元および3次元）： これらに対して、mnormはガウス・レジェンド積分（Gauss-Legendre quadrature）という巧妙な数学的トリックを使用します。著者らが他の一般的なツールと比較したところ、mnormは大幅に高速であることが分かりました。例えば、2次元の雲の確率を計算する場合、100万回の計算に対してmnormは0.57秒かかりましたが、競合するツールは0.75秒かかりました。さらに驚くべきことに、「マルチコア」機能（8つのコンピュータの頭脳を同時に使用する機能）をオンにすると、mnormは0.16秒まで短縮されましたが、競合するツールは依然として1.11秒でした。
2. 巨大な雲（高次元）： 16、64、あるいは512次元といった雲に対しても、mnormは利用可能な最高のツールと互角に渡り合います。512次元のテストにおいて、mnormは2.54秒で確率を計算しました。これは他のトップクラスのツールと同等の性能であり、mnormが小さな問題だけでなく、大きな問題にも対応できることを証明しています。

非正規分布の仲間たち
現実の世界は、常に完璧に「正規」であるとは限りません。データが歪んでいたり、奇妙な裾（テイル）を持っていたりすることもあります。mnormは、「ガウス・コピュラ（Gaussian copula）」を用いることでこれに対処できます。これは「翻訳機」のようなものです。奇妙で非正規なデータ（例えば、リスが餌箱を訪れる回数など）を取り込み、数学が機能するように正規分布の言語へと翻訳します。このパッケージは、これらの奇妙な分布のパラメータに対する傾斜を計算することも可能であり、これは他の多くのツールが持っていない機能です。

結論
この論文は、宇宙の新しい法則を発見したと主張しているわけではありません。代わりに、実用的で極めて効率的なツールを提供しています。著者らは、mnormを使用することで、研究者がより優れた、より高速な経済モデルを構築できることを示唆しています。彼らはシミュレーションを通じて、彼らの手法が代替案よりも高速であることを証明し、誰でもその作業を検証できるように正確な数学的公式（「秘伝のソース」）を提供しました。

端的に言えば、mnormは統計家に、自動車のターボエンジンを与えてくれるようなものです。目的地を変えるわけではありませんが、地形の明確な地図を携えながら、記録的な速さで目的地へと到達させてくれます。2つの変数の単純な問題であっても、512次元の巨大なパズルであっても、このパッケージは、それが答えを得るための最も速い方法であることを示しています。