互いに独立な指数分布の比の分布の導出

公開日: 更新日:

【2023年3月5週】 【B000】数理統計学 【B040】連続型の確率分布

この記事をシェアする
  • B!
サムネイル画像

本稿では、互いに独立な指数分布の比の分布が連続一様分布になることを証明しています。変数変換後の周辺確率密度関数を求める問題としては手ごろな問題です。

なお、閲覧にあたっては、以下の点にご注意ください。

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。

【命題】互いに独立な指数分布の比の分布

【命題】
互いに独立な指数分布の比の分布
A Ratio of Independent Exponential Random Variables of Exponential Distribution

確率変数 $X,Y$ がそれぞれ独立に、指数分布 \begin{align} X \sim \mathrm{Ex} \left(\lambda\right) \quad Y \sim \mathrm{Ex} \left(\lambda\right) \end{align} に従うとき、 2変数を用いてできる新たな確率変数 \begin{align} U=\frac{X}{X+Y} \end{align} は、 連続一様分布 \begin{align} U \left(0,1\right) \end{align} に従う。

証明

証明

確率変数の独立性の定義 $f \left(x,y\right)=g \left(x\right) \cdot h \left(y\right)$ より、同時確率密度関数は、 \begin{align} f \left(x,y\right)&=\lambda e^{-\lambda x} \cdot \lambda e^{-\lambda y}\\ &=\lambda^2e^{-\lambda \left(x+y\right)} \end{align} ここで、以下のように変数変換を行うと、 \begin{align} \left\{\begin{matrix}u=\frac{x}{x+y}\\v=x+y\\\end{matrix}\right.\Leftrightarrow \left\{\begin{matrix}u=x=uv\\y=v-uv\\\end{matrix}\right. \end{align} 変数変換のヤコビアンは、 \begin{align} \left|J\right|&= \left|\begin{matrix}\frac{\partial x}{\partial u}&\frac{\partial y}{\partial u}\\\frac{\partial x}{\partial v}&\frac{\partial y}{\partial v}\\\end{matrix}\right|\\ &= \left|\begin{matrix}v&-v\\u&1-u\\\end{matrix}\right|\\ &= \left|v \left(1-u\right)+uv\right|\\ &=v \end{align} 変数変換後の同時確率密度関数の公式 $k \left(u,v\right)=f \left\{x \left(u,v\right),y \left(u,v\right)\right\} \left|J\right|$ より、 \begin{align} g \left(u,v\right)&=f \left(uv,v-uv\right) \cdot v\\ &=\lambda^2e^{-\lambda v} \cdot v \end{align} 確率変数 $U$ の周辺確率密度関数を求めると、 \begin{align} h \left(u\right)&=\int_{0}^{\infty}{\lambda^2 \cdot v e^{-\lambda v}dv}\\ &=\lambda\int_{0}^{\infty}{\lambda v \cdot e^{-\lambda v}dv} \end{align} 部分積分法により、 \begin{align} \int_{0}^{\infty}{\lambda v \cdot e^{-\lambda v}dv}&= \left[-\lambda v \cdot e^{-\lambda v}\right]_0^\infty-\int_{0}^{\infty}{-e^{-\lambda v}dv}\\ &=-\lim_{x\rightarrow\infty}{v \cdot e^{-\lambda v}}+\int_{0}^{\infty}{e^{-\lambda v}dv} \end{align} ロピタルの定理より、 \begin{align} \lim_{x\rightarrow\infty}{xe^{-\lambda x}}=\lim_{x\rightarrow\infty}{\frac{x}{e^{\lambda x}}}=\lim_{x\rightarrow\infty}{\frac{1}{\lambda e^{\lambda x}}}=0 \end{align} したがって、 \begin{align} \int_{0}^{\infty}{\lambda v \cdot e^{-\lambda v}dv}&=\int_{0}^{\infty}{e^{-\lambda x}dx}\\ &= \left[-\frac{1}{\lambda}e^{-\lambda x}\right]_0^\infty\\ &=-\frac{1}{\lambda} \left(\lim_{x\rightarrow\infty}{e^{-\lambda x}}-e^0\right)\\ &=\frac{1}{\lambda} \end{align} したがって、 \begin{align} h \left(u\right)=\lambda \cdot \frac{1}{\lambda}=1 \end{align} これは、連続一様分布 \begin{align} U \left(0,1\right) \end{align} の確率密度関数なので、 確率密度関数の一意性から、確率変数 $U$ は連続一様分布に従う。 $\blacksquare$

参考文献

  • 野田 一雄, 宮岡 悦良 著. 入門・演習数理統計. 共立出版, 1990, p.143-144 練習問題 ex.3.8.3

関連記事

自己紹介

自分の写真

yama

大学時代に読書の面白さに気づいて以来、読書や勉強を通じて、興味をもったことや新しいことを学ぶことが生きる原動力。そんな人間が、その時々に学んだことを備忘録兼人生の軌跡として記録しているブログです。

このブログを検索

ブログ アーカイブ

QooQ