QCプラネッツ 品質のプロフェッショナルを育成するサイト

2変数の確率変数の変換がよくわかる(事例1)

統計学

「確率変数の変換が、わからない、解けない?」、「t分布、F分布の確率密度関数への導出がわからない」と困っていませんか?

こういう疑問に答えます。

本記事のテーマ

2変数の確率変数の変換がよくわかる(事例1)
  • ①2変数の確率変数の変換の基本をマスターする
  • ➁実例を使って理解する!

①2変数の確率変数の変換の基本をマスターする

確率変数の変換は高校数学でほぼイケます!大丈夫!

確率変数の変換は難しいけど、
理解しないと、正規分布、t分布、χ2乗分布、F分布との関係が理解できないから困っている!
確率変数の変換は高校数学でほぼイケます!大丈夫!
1つ条件があります!

それは、

公式暗記より、実演でマスターした方が速い!
1つ解法で解ける解法で、たくさんの例題を見る方がマスターは速い!

慣れてきたら、公式を見ましょう。

2変数の確率変数の変換の基本をマスターする

関連記事に2変数の確率変数の変換の求め方をわかりやすく解説しています。

【まとめ】2変数の確率変数の変換がよくわかる
2変数の確率変数の変換が計算できますか?本記事では,理解が難しい公式をそのまま使わずに,高校数学で十分解ける解法を解説します。今回は変換したいパターンをすべてを解説!教科書よりわかりやすく、 ほぼ高校数学でイケる方法で解説! t分布、F分布の確率密度関数を導出したい方は必読な記事です。

同じ1つの解法でイケますので、ご安心ください。

2変数の確率変数の変換の求め方

1変数の確率変数の変換方法と同様に決まった解法があります。

変数\(x,y\)を変数\(z,w\)に変換するとします。

  1. \(x=x(z,w),y=(z,w)\)の式を\(z=z(x,y),w=w(x,y)\)の式に直す
  2. \(f(x,y)dxdy\)=\(f(x(z,w),y(z,w)|det J| dzdw\)に変換する
  3. 2変数\(z,w\)の同時確率密度関数\(g(z,w)\)は
    \(g(z,w)=f(x(z,w),y(z,w)|det J| \)で求まる。
  4. 実際は\(z,w\)のうち、どちらかは不要な変数なので、片方の変数で積分して、残りの変数についての周囲確率密度関数
    (例えば \(g(z)= \displaystyle \int_{w_1}^{w_2} g(z,w)dw \))
    を計算する。

ここで、注意点があります。
Jは
J=\(\begin{pmatrix}
\frac{\partial x}{ \partial z} & \frac{\partial x}{\partial w} \\
\frac{\partial y}{\partial z} & \frac{\partial y}{\partial w}
\end{pmatrix}\)

また、\(det J\)は行列式ヤコビアンといいますね。

A=\(\begin{pmatrix}
a & b \\
c & d
\end{pmatrix}\)
のとき、行列式ヤコビアン\(det A\)は、
\(det A=ad-bc\)
で計算できます。

計算力が求められる場合がありますが、基本は高校数学でイケます!

では、実践編に入ります。最初は簡単な式から行きます!

➁実例を使って理解する!

QCプラネッツでは、5つの事例を関連記事で紹介していきます。ご確認ください。

  1. 簡単な関数の変換事例
  2. t分布の確率密度関数の導出
  3. F分布の確率密度関数の導出
  4. 1変数でZ=XY(積)の場合の変換方法
  5. 1変数でZ=X/Y(商)の場合の変換方法

今回は、その1「簡単な関数の変換事例」です。

(1) 簡単な関数の変換事例

【例題】
確率変数(\(X,Y\))の同時確率密度関数\(f(x,y)\)が
\(f(x,y)\)=\(\frac{1}{2}xy^2\) (0 ≤ \(x\) ≤ 2, 0 ≤ \(y\) ≤ 1)
に対して、確率変数\(Z,W\)を
\(Z=2X+Y\)
\(W=X-2Y\)
と定義した場合の、確率変数\(Z,W\)についての同時確率密度関数\(g(z,w)\)を求めよ。

解いていきましょう。解法は、

  1. \(x=x(z,w),y=(z,w)\)の式を\(z=z(x,y),w=w(x,y)\)の式に直す
  2. \(f(x,y)dxdy\)=\(f(x(z,w),y(z,w)|det J| dzdw\)に変換する
  3. 2変数\(z,w\)の同時確率密度関数\(g(z,w)\)は
    \(g(z,w)=f(x(z,w),y(z,w)|det J| \)で求まる。
  4. 実際は\(z,w\)のうち、どちらかは不要な変数なので、片方の変数で積分して、残りの変数についての周囲確率密度関数
    (例えば \(g(z)= \displaystyle \int_{w_1}^{w_2} g(z,w)dw \))
    を計算する。

ですから、1つずつ行きましょう。

(i)\(x=x(z,w),y=(z,w)\)の式を\(z=z(x,y),w=w(x,y)\)の式に直す

\(x=x(z,w),y=y(z,w)\)に直します。
\(x\)=\(\frac{1}{5}(2z+w)\)
\(y\)=\(\frac{1}{5}(z-2w)\)
連立方程式から求められます。

(ii)\(f(x,y)dxdy\)=\(f(x(z,w),y(z,w)|det J| dzdw\)に変換する

次に、ヤコビアン行列から行列式ヤコビアンを求めます。

\(det J\)=\(\frac{2}{5}・(-\frac{2}{5})\)-\(\frac{1}{5}・\frac{1}{5}\)
=\(-\frac{1}{5}\)
で計算できます。

(iii)2変数\(z,w\)の同時確率密度関数\(g(z,w)\)を導出

代入すると、

\(f(x,y)dxdy\)=\(f(x(z,w),y(z,w)|det J| dzdw\)
=\(\frac{1}{2} × \frac{1}{5}(2z+w) ×\frac{1}{25}(z-2w)^2 ×|-\frac{1}{5}| dzdw\)
=\(\frac{1}{1250} (2z+w) (z-2w)^2 dzdw\)

よって、同時確率密度関数\(g(z,w)\)は
\(g(z,w)= \frac{1}{1250} (2z+w) (z-2w)^2 \)
と計算できます。

なお、ここから\(z\)または、\(w\)だけの周辺確率分布関数が必要なら、不要な変数について積分が必要となります。

今回の\(f(x,y)\)は簡単な式ですが、正規分布やχ2乗分布の確率分布関数でも同様の解法で変換していきます!
簡単だ! 解けるぞ!と自信持ってください!

いろいろな関数を使って、確率変数の変換を見て慣れていきましょう!

本記事の内容は、ほぼ高校数学で解けましたね!

まとめ

「2変数の確率変数の変換がよくわかる(事例1)」を解説しました。

  • ①2変数の確率変数の変換の基本をマスターする
  • ➁実例を使って理解する!


Warning: count(): Parameter must be an array or an object that implements Countable in /home/qcplanets/qcplanets.com/public_html/wp-content/themes/m_theme/sns.php on line 119

    Warning: Invalid argument supplied for foreach() in /home/qcplanets/qcplanets.com/public_html/wp-content/themes/m_theme/sns.php on line 122
error: Content is protected !!