連続型確率変数のモード
確率空間\(\left( \Omega ,\mathcal{F},P\right) \)に加えて絶対連続型の確率変数\begin{equation*}X:\Omega \rightarrow \mathbb{R} \end{equation*}が与えられているものとします。加えて、確率変数\(X\)の確率分布が確率密度関数\begin{equation*}f_{X}:\mathbb{R} \rightarrow \mathbb{R} \end{equation*}によって記述されているものとします。つまり、\begin{equation*}
\forall B\in \mathcal{B}\left( \mathbb{R} \right) :P\left( X\in B\right) =\int_{B}f_{X}\left( x\right) dx
\end{equation*}が成り立つということです。ただし、\(\mathcal{B}\left( \mathbb{R} \right) \)は\(\mathbb{R} \)上のボレル集合族です。確率変数\(X\)の分布関数が、\begin{equation*}F_{X}:\mathbb{R} \rightarrow \mathbb{R} \end{equation*}である場合、以下の関係\begin{equation*}
\forall x\in \mathbb{R} :F_{X}\left( x\right) =\int_{-\infty }^{x}f_{X}\left( t\right) dt
\end{equation*}もまた成立します。
離散型の確率変数\(X\)に対しては、実数\(m\in \mathbb{R} \)が\(X\)のモードであることを、\begin{equation*}\forall x\in \mathbb{R} :P\left( X=m\right) \geq P\left( X=x\right)
\end{equation*}が成り立つこととして定義しました。一方、絶対連続型の確率変数\(X\)に関しては、それが特定の値をとる確率は常にゼロであるため、上の命題は恒真式\begin{equation*}\forall x\in \mathbb{R} :0\geq 0
\end{equation*}になってしまいます。このような事情を踏まえると、絶対連続型の確率変数のモードについては、それを別の形で定義する必要があります。
絶対連続型の確率変数\(X:\Omega \rightarrow \mathbb{R} \)の確率密度関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)が与えられたとき、実数\(m\in \mathbb{R} \)が以下の条件\begin{equation*}\forall x\in \mathbb{R} :f_{X}\left( m\right) \geq f_{X}\left( x\right)
\end{equation*}を満たす場合には、つまり、\(m\)が\(f_{X}\)の大域的最大点である場合には、これを\(X\)のモード(mode)と呼びます。
\begin{array}{cl}
x+\frac{1}{2} & \left( if\ x\in X\left( \Omega \right) \right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。\(f_{X}\left( x\right) \)は\(X\left( \Omega \right) \)上において狭義単調増加であるため、\(X\)のモードは\(1\)です。
微分を用いたモードの特定方法
絶対連続型確率変数\(X:\Omega \rightarrow \mathbb{R} \)の確率密度関数\(f_{X}:\mathbb{R} \rightarrow \mathbb{R} \)が与えられた状況を想定します。確率密度関数の定義より、\begin{eqnarray*}\forall x &\in &X\left( \Omega \right) :f_{X}\left( x\right) \geq 0 \\
\forall x &\not\in &X\left( \Omega \right) :f_{X}\left( x\right) =0
\end{eqnarray*}がともに成り立つため、\(X\)のモードは\(X\left( \Omega\right) \)上の点です。モードの定義より、\(X\)のモードは最大化問題\begin{equation*}\max_{x\in X\left( \Omega \right) }\ f_{X}\left( x\right)
\end{equation*}の解です。特に、\(X\left(\Omega \right) \)が\(\mathbb{R} \)上のコンパクト集合であるとともに\(f_{X}\)が\(X\left( \Omega \right) \)上の連続関数である場合には、最大値・最小値の定理より\(f_{X}\)は\(X\left( \Omega \right) \)上に最大点を持つため、\(X\)のモードが存在することが保証されます。さらに、以下の手順を通じて\(X\)のモードを求めることができます。
- \(X\left( \Omega \right) \)の内点の中でも、局所最大化のための必要条件\begin{eqnarray*}&&\left( a\right) \ f_{X}^{\prime }\left( x\right) =0 \\&&\left( b\right) \ f_{X}^{\prime \prime }\left( x\right) <0
\end{eqnarray*}を満たす点\(x\)をすべて特定する。 - \(X\left( \Omega \right) \)の内点の中でも、\(f_{X}\)が微分可能ではない点をすべて特定する。
- \(X\left( \Omega \right) \)の境界点をすべて特定する。
- 以上の点の中でも\(f_{X}\)の値を最大化する点が最大点であり、したがってそれが\(X\)のモードである。
\begin{array}{cl}
\frac{3}{64}x^{2}\left( 4-x\right) & \left( if\ x\in X\left( \Omega \right)
\right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。\(X\left( \Omega \right) \)は有界閉区間であり\(f_{X}\)は\(X\left( \Omega\right) \)上の連続関数であるため、\(X\)のモードが存在します。\(f_{X}\)は\(X\left(\Omega \right) \)の内部\(\left( 0,4\right) \)において微分可能であり、導関数\(f_{X}^{\prime }:\mathbb{R} \supset \left( 0,4\right) \rightarrow \mathbb{R} \)はそれぞれの\(x\in \left( 0,4\right) \)に対して、\begin{eqnarray*}f_{X}^{\prime }\left( x\right) &=&\frac{d}{dx}\frac{3}{64}x^{2}\left(
4-x\right) \\
&=&\frac{3}{64}\frac{d}{dx}x^{2}\left( 4-x\right) \\
&=&\frac{3}{64}\left( 8x-3x^{2}\right)
\end{eqnarray*}を定めます。また、\(f_{X}\)は\(X\left( \Omega \right) \)の内部\(\left(0,4\right) \)において2階微分可能であり、2階導関数\(f_{X}^{\prime \prime }:\mathbb{R} \supset \left( 0,4\right) \rightarrow \mathbb{R} \)はそれぞれの\(x\in \left( 0,4\right) \)に対して、\begin{eqnarray*}f_{X}^{\prime \prime }\left( x\right) &=&\frac{d}{dx}\frac{3}{64}\left(
8x-3x^{2}\right) \\
&=&\frac{3}{64}\left( 8-6x\right)
\end{eqnarray*}を定めます。局所最大化のための1階の必要条件は、\begin{equation*}
\frac{3}{64}\left( 8x-3x^{2}\right) =0
\end{equation*}であり、これを解くと、\begin{equation*}
x=0,\frac{8}{3}
\end{equation*}を得ます。これらの中でも\(\left( 0,4\right) \)上の点は\(\frac{8}{3}\)です。さらに、\begin{eqnarray*}f^{\prime \prime }\left( \frac{8}{3}\right) &=&\frac{3}{64}\left( 8-6\cdot
\frac{8}{3}\right) \\
&=&-\frac{3}{8} \\
&<&0
\end{eqnarray*}が成り立つため、点\(\frac{8}{3}\)は局所最大化のための2階の必要条件を満たします。したがって、点\(\frac{8}{3}\)は最大点の候補であるとともに、\begin{equation*}f_{X}\left( \frac{8}{3}\right) =\frac{3}{64}\left( \frac{8}{3}\right)
^{2}\left( 4-\frac{8}{3}\right) =\frac{4}{9}
\end{equation*}です。\(X\left( \Omega \right) \)の境界点である\(0\)と\(4\)においては、\begin{eqnarray*}f_{X}\left( 0\right) &=&\frac{3}{64}\cdot 0^{2}\left( 4-0\right) =0 \\
f_{X}\left( 4\right) &=&\frac{3}{64}\cdot 4^{2}\left( 4-4\right) =0
\end{eqnarray*}となります。したがって、最大点すなわち\(X\)のモードは\(\frac{8}{3}\)であることが明らかになりました。
連続型確率変数のモードの個数
絶対連続型確率変数のモードは一意的に定まるとは限りません。以下の例より明らかです。
\begin{array}{cl}
1 & \left( if\ x\in X\left( \Omega \right) \right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。\(f_{X}\)は\(\left[ 0,1\right] \)上において定数関数であるため、\(\left[ 0,1\right] \)上の任意の値が\(X\)のモードです。
演習問題
\begin{array}{cl}
\frac{3}{4}x^{2}\left( 2-x\right) & \left( if\ x\in X\left( \Omega \right)
\right) \\
0 & \left( if\ x\not\in X\left( \Omega \right) \right)
\end{array}\right.
\end{equation*}を定めるものとします。\(X\)のモードを求めてください。
\begin{array}{cl}
\frac{2}{25}x\left( 4-x\right) & \left( if\ 0\leq x<3\right) \\
\frac{2}{25}x & \left( if\ 3\leq x\leq 4\right) \\
0 & \left( otherwise\right)
\end{array}\right.
\end{equation*}を定めるものとします。\(X\)のモードを求めてください。
プレミアム会員専用コンテンツです
【ログイン】【会員登録】