標準正規分布に従う2個の確率変数の商が従う確率密度関数を計算してみた。

はじめに

前の記事で、以下の確率密度関数について計算してみました。

和、差、積の次は商だと思います(※個人の感想です。)ので、標準正規分布に従う2個の確率変数の商が従う確率密度関数を計算してみることにしました。

問題の定義

まず、解くべき問題を以下のように定義します。

標準正規分布に従う2個の確率密度変数がある。このとき、によって定義される確率変数が従う確率密度関数を求めよ。
ただし、とする。

ということで積の計算に続きまして、が負の場合はとりあえず考えないことにします。(´・ω・｀)

サクサクと計算。

の符号による場合分け。

問題を定義したところで、サクサクと計算します。

まず、積分の領域を定めます。

の場合、積分を実行するの範囲を2次元平面上にプロットすると、以下の図の水色で示した部分になります(の場合については後述)。

となるについて積分を実行していきますが、の符号によって場合分けをする必要がありそうです。また、の場合はになると考えることにします。(´・ω・｀)

が負の場合はとりあえず考えないことにしましたが、の場合のことが気にかかるので、をで割った商で定義される確率変数のとる値が以下となる確率をの符号によって場合分けすると、( $1$ )式のように表すことができます。

(1)

ちなみに( $1$ )式にを代入すると、

(2)

となります。( $2$ )式の積分を実行するの範囲を2次元平面上にプロットすると、以下の図の水色で示した部分になりますので、の場合も( $1$ )式でまとめて考えることができそうです。

( $1$ )式の右辺の第1項は( $3$ )式のように計算できます。上図の(a)の向きに最初に積分するイメージになります。なお、( $3$ )式の右辺ではをと書き換えています。

(3)

さらに、( $3$ )式でと置きます。すると、( $3$ )式の右辺は、

(4)

と計算できます。

次に、( $1$ )式の右辺の第2項を計算してみます。

( $1$ )式の右辺の第2項は、

(5)

と計算できます。上図の(b)の向きに最初に積分するイメージになります。これを(

4

)式と似たような式に変形できないか考えてみます。

まず、は偶関数ですので、( $5$ )式の右辺は( $6$ )式のように書き換えることができます。

(6)

次に、

と置きます。すると、(

6

)式の右辺は、

(7)

となりますので、さらに

と置きます。すると、(

7

)式の右辺は、

(8)

となります。

( $8$ )式のをに置き換えると( $4$ )式の右辺と同じ式が得られます。

したがって、は( $9$ )式で表されます(をに置き換えています)。

(9)

( $9$ )式の及びについての積分の順序は入れ替えることができます。そこで、積分の順序を入れ替えてからで微分すると、における確率密度関数が以下のように計算できます。

(10)

最後の仕上げ。

和・差及び積の場合とは異なり、( $10$ )式が初等的な積分計算ができそうな式になりました。そこで、( $10$ )式でとおきます。すると、となりますので( $11$ )式のように計算できます。

(11)

( $11$ )式は標準コーシー分布(Standard Cauchy Distribution)という確率分布です。一般的なコーシー分布は位置母数(location parameter)及び四分位範囲を用いて、( $12$ )式で表されます。

(12)

( $11$ )式は、( $12$ )式でと置いたものになります。

ToDoまたは検討が必要かもしれない事項

の場合はこちらに書きましたので、ご覧いただけると幸いです。

また、の場合については上図の軸上でのについての区間における定積分を行うことになります。( $1$ )式右辺の第1項にの区間が、第2項にの区間がそれぞれ含まれていますが、これで正しいのかどうかついては検証できていません。

まとめ

商の計算はのときにちょっと歯切れが悪くなってしまうのが気になるところではありますが、計算自体は積の場合と比較すると難しくはないと思います。

コーシー分布は式が簡単な割には期待値が計算できなかったりする等、かなり面白い性質を持つ関数だと思います。詳細につきましてはWikipediaの記事が詳しいので、参考にしていただけると幸いです。

この記事は以上です。

References / 参考文献

Cauchy distribution(英語版Wikipedia, 日本語版Wikipediaのページはこちら。)

おまけ

( $1$ )式は最初は、

(13)

のように書いていましたが、のケースが考慮できていないような気がしたので、書き換えました。