出典(authority):フリー百科事典『ウィキペディア(Wikipedia)』「2017/06/09 10:39:45」(JST)
この記事は検証可能な参考文献や出典が全く示されていないか、不十分です。 出典を追加して記事の信頼性向上にご協力ください。(2011年11月) |
確率分布(かくりつぶんぷ, 英: probability distribution)は、確率変数の各々の値に対して、その起こりやすさを記述するものである。日本工業規格では、「確率変数がある値となる確率,又はある集合に属する確率を与える関数」と定義している[1]。
例えば、「サイコロを二つ振ったときの出た目の和」は確率変数であるが、その分布は次の表のように書くことができる。
出た目の和 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
---|---|---|---|---|---|---|---|---|---|---|---|
確率 | 1/36 | 2/36 | 3/36 | 4/36 | 5/36 | 6/36 | 5/36 | 4/36 | 3/36 | 2/36 | 1/36 |
すなわち、確率分布は値に確率を対応させる関数ということができる。確率変数がこのように離散的な値しかとらないときは上のような理解で十分である。しかし、例えば「次に電話がなるまでの時間」といった連続的な値をとる確率変数の分布はこのような形では表現できず、測度の概念が必要になる。
「次に電話がなるまでの時間」の分布の一部を表にしたとき、次のようになったとする。
次の電話 | 1時間以内 | 1 - 2時間後 | 2 - 3時間後 | 3 - 4時間後 | 4時間以上先 |
---|---|---|---|---|---|
確率 | 1/2 | 1/4 | 1/8 | 1/16 | 1/16 |
この表だけでは「次に電話がなるまでの時間」を正確に記述しているとはいえない。完全なものにするためには、次の電話がa - b時間後になる確率をすべて記述する必要がある。「次に電話がなるまでの時間」を X と書くことにすれば、この確率は
と書ける。累積分布関数(るいせきぶんぷかんすう cumulative distribution function, CDF) FX を
,&t<0\end{cases}}}</annotation>
</semantics></math>
で定めれば、
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo><</mo> <mi>X</mi> <mo>≤</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msub> <mi>F</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>−</mo> <msub> <mi>F</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>a</mi> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P(a<X\leq b)=F_{X}(b)-F_{X}(a)}</annotation> </semantics> </math>
のように、一変数関数で分布を表現できるので便利である。さらに、FX の導関数 fX は確率密度関数(frequency functionまたは probability density function(PDF)) と呼ばれ、確率は積分を用いて
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo><</mo> <mi>X</mi> <mo>≤</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msubsup> <mo>∫</mo> <mrow class="MJX-TeXAtom-ORD"> <mi>a</mi> </mrow> <mrow class="MJX-TeXAtom-ORD"> <mi>b</mi> </mrow> </msubsup> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>t</mi> <mo stretchy="false">)</mo> <mi>d</mi> <mi>t</mi> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P(a<X\leq b)=\int _{a}^{b}f_{X}(t)dt}</annotation> </semantics> </math>
と書ける。
通常、連続値をとる確率変数の分布は確率密度関数を用いて記述される。なぜかというと、確率密度関数は初等関数で書けるが、累積分布関数は書けない場合が多いからである。
公理主義的な確率論においては、d次元ベクトル値確率変数の確率分布とは、その確率変数の引き起こす像測度のことである。この測度は d次元ユークリッド空間上の確率測度であり、ユークリッド空間の部分集合に対して、確率変数の値がその集合に入る確率を与える関数となる。
単に確率分布というときは、d次元ユークリッド空間などのよく使われる可測空間上で定義された確率測度のことをいう。ただの確率測度と違って空間に散らばっている様子がグラフなどの目に見える形で表現できるので「分布」と呼ばれる。
確率論で、確率変数の分布を考えるのは、その変数だけを確率論的な議論の対象にしたい場合である。例えば、確率変数がある値を取る確率や、期待値、分散といった量は変数の分布が分かれば計算できる量である。 逆に分布を考えることによって隠れた変数ωと確率変数との対応関係は失われてしまい、他の確率変数との関連性も不明になる。例えば、確率変数XとYの分布がそれぞれPXとPYのように与えられたとしても、ふたつの変数の関連性は分からないので、X+Y がある値を取る確率や、積 X Y の期待値、X+Y の分散といった量は計算できない。このような量を計算したいときは、XとYの結合分布が必要となる。
よく使われる確率分布にはそれぞれ名前がついており性質がよく研究されている。このような分布をもつ確率変数に対して研究の結果を利用することができる。例えば、確率変数の分布が平均 0 分散 1 の正規分布だった場合、その変数が 2 以上の値を取る確率は数表から 2.28% である。
実数値確率変数 X の分布 PX は像測度であるから
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>X</mi> <mo>∈</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>,</mo> <mtext> </mtext> <mtext> </mtext> <mtext> </mtext> <mi>A</mi> <mo>∈</mo> <mi>B</mi> <mo stretchy="false">(</mo> <mrow class="MJX-TeXAtom-ORD"> <mi mathvariant="double-struck">R</mi> </mrow> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P_{X}(A)=P(X\in A),\ \ \ A\in B(\mathbb {R} )}</annotation> </semantics> </math>
で定義される確率測度のことである。
同様に Rd 値確率変数 X の分布 PX は
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>X</mi> <mo>∈</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>,</mo> <mtext> </mtext> <mtext> </mtext> <mtext> </mtext> <mi>A</mi> <mo>∈</mo> <mi>B</mi> <mo stretchy="false">(</mo> <msup> <mrow class="MJX-TeXAtom-ORD"> <mi mathvariant="double-struck">R</mi> </mrow> <mrow class="MJX-TeXAtom-ORD"> <mi>d</mi> </mrow> </msup> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P_{X}(A)=P(X\in A),\ \ \ A\in B(\mathbb {R} ^{d})}</annotation> </semantics> </math>
で定義される確率測度のことである。
確率変数 X の分布がある確率分布 μ に一致するとき、X は μ に従う確率変数であるという。例えば、「X は平均 0 分散 1 の正規分布に従う」のように使い、これを記号で
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>X</mi> <mo>∼</mo> <mi>N</mi> <mo stretchy="false">(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle X\sim N(0,1)}</annotation> </semantics> </math>
のように書く。
実数値確率変数 X の分布関数(ぶんぷかんすう, distribution function)あるいは、一次元確率分布 PX の分布関数とは
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>F</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>x</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>X</mi> <mo>≤</mo> <mi>x</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mo stretchy="false">(</mo> <mo>−</mo> <mi mathvariant="normal">∞</mi> <mo>,</mo> <mi>x</mi> <mo stretchy="false">]</mo> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle F_{X}(x)=P(X\leq x)=P_{X}((-\infty ,x])}</annotation> </semantics> </math>
で与えられる関数 FX のことである。累積分布関数 (るいせきぶんぷかんすう, cumulative distribution function) ともいう。
定義から分布関数は右連続であるが、左連続かどうかはわからない。これが連続であるときに確率分布は連続であるという。分布関数が高々可算個の値しかとらない場合は確率分布は離散であるという。
確率分布 PX が絶対連続であるというのは、任意の(ルベーグ測度に関しての)零集合 A にたいして、
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mn>0</mn> <mspace width="negativethinmathspace" /> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P_{X}(A)=0\!}</annotation> </semantics> </math>
が成り立つことを言う。これは測度の絶対連続性と同じである。 確率分布 PX が絶対連続のとき、測度論のラドン-ニコディムの定理によりラドン-ニコディム微分 fX が存在する。この fX のことを確率密度関数(frequency functionまたは probability density function(PDF)) と呼ぶ。 PX は確率密度関数を用いて
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>P</mi> <mo stretchy="false">(</mo> <mi>X</mi> <mo>∈</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>A</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msub> <mo>∫</mo> <mrow class="MJX-TeXAtom-ORD"> <mi>A</mi> </mrow> </msub> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>x</mi> <mo stretchy="false">)</mo> <mi>d</mi> <mi>x</mi> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P(X\in A)=P_{X}(A)=\int _{A}f_{X}(x)dx}</annotation> </semantics> </math>
と書くことができる。とくに A が区間の場合は
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo><</mo> <mi>X</mi> <mo><</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo>≤</mo> <mi>X</mi> <mo><</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo><</mo> <mi>X</mi> <mo>≤</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>a</mi> <mo>≤</mo> <mi>X</mi> <mo>≤</mo> <mi>b</mi> <mo stretchy="false">)</mo> <mo>=</mo> <msubsup> <mo>∫</mo> <mrow class="MJX-TeXAtom-ORD"> <mi>a</mi> </mrow> <mrow class="MJX-TeXAtom-ORD"> <mi>b</mi> </mrow> </msubsup> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>x</mi> <mo stretchy="false">)</mo> <mi>d</mi> <mi>x</mi> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle P(a<X<b)=P(a\leq X<b)=P(a<X\leq b)=P(a\leq X\leq b)=\int _{a}^{b}f_{X}(x)dx}</annotation> </semantics> </math>
となる。区間の端点はいれてもいれなくても確率は同じである。
確率分布 PX が離散のときに確率密度関数に対応する関数として確率質量関数(かくりつしつりょうかんすう, probability mass function)(単に確率関数(かくりつかんすう、probability functionまたはprobability mass function、random function、stochastic function)ともいう)を使うことができる。確率変数 X のとる値の集合が S = {x1, x2, ...} だとすると確率質量関数は
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>i</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mo>=</mo> <mi>P</mi> <mo stretchy="false">(</mo> <mi>X</mi> <mo>=</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>i</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mo>=</mo> <msub> <mi>P</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mo fence="false" stretchy="false">{</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>i</mi> </mrow> </msub> <mo fence="false" stretchy="false">}</mo> <mo stretchy="false">)</mo> <mspace width="negativethinmathspace" /> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle f_{X}(x_{i})=P(X=x_{i})=P_{X}(\{x_{i}\})\!}</annotation> </semantics> </math>
で定まる関数 fX のことである。
二つ以上の変数の確率分布関数を、多次元分布関数と呼ぶ。また、二つ以上の変数の組の確率分布のことを同時分布、同時確率分布 (joint probability distribution) という[2]。
二変数の確率分布関数を、二次元分布関数と呼ぶ[2]。
複数の確率変数の挙動を多次元の確率分布で表したものを同時分布(どうじぶんぷ, simultaneous distribution)という。同時分布から各変数の分布だけを取り出したものを周辺分布(しゅうへんぶんぷ, marginal distribution)と呼ぶ。日本工業規格では、周辺(確率)分布(しゅうへん(かくりつ)ぶんぷ marginal probability distribution)を、「k次元確率変数の部分集合である k1変数の同時分布」と定義している[3] 。
代表的な確率分布は離散型または絶対連続型のもののみを挙げる。他には離散でも連続でもないもの、連続であるが絶対連続ではないものなどが考えられるが通常現れる分布ではない。
サイコロを投げた時に出る目の数字など、確率変数が離散的な値をとる場合の確率分布は離散型確率分布である。
離散型の分布は母数と台 S と確率質量関数 f で特徴付けられる。台というのは確率変数のとる値の集合のことである。
ある地点での通行人の体重など、確率変数が連続的な場合の確率分布は連続型確率分布である。
絶対連続な分布は母数と台と確率密度関数 f で特徴付けられる
確率変数の確率分布が与えられると、その変数に関する確率・期待値・分散などが以下のように計算できる。
X は連続値をとる確率変数で密度関数は fX であるとする。 Y は離散値をとる確率変数で台は S = {y1, y2, ...} で質量関数は fY であるとする。
関数 g が与えられたときに g(X) と g(Y) の期待値は
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>g</mi>
<mo stretchy="false">(</mo>
<mi>X</mi>
<mo stretchy="false">)</mo>
<mo stretchy="false">]</mo>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<mi>g</mi>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo stretchy="false">)</mo>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo stretchy="false">)</mo>
<mi>d</mi>
<mi>x</mi>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle E[g(X)]=\int _{-\infty }^{+\infty }g(x)f_{X}(x)dx}</annotation>
</semantics>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>g</mi>
<mo stretchy="false">(</mo>
<mi>Y</mi>
<mo stretchy="false">)</mo>
<mo stretchy="false">]</mo>
<mo>=</mo>
<munder>
<mo>∑</mo>
<mrow class="MJX-TeXAtom-ORD">
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo>∈</mo>
<mi>S</mi>
</mrow>
</munder>
<mi>g</mi>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo stretchy="false">)</mo>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo stretchy="false">)</mo>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle E[g(Y)]=\sum _{y_{k}\in S}g(y_{k})f_{Y}(y_{k})}</annotation>
</semantics>
</math>
とくに
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>X</mi>
<mo stretchy="false">]</mo>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<mi>x</mi>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo stretchy="false">)</mo>
<mi>d</mi>
<mi>x</mi>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle E[X]=\int _{-\infty }^{+\infty }xf_{X}(x)dx}</annotation>
</semantics>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>Y</mi>
<mo stretchy="false">]</mo>
<mo>=</mo>
<munder>
<mo>∑</mo>
<mrow class="MJX-TeXAtom-ORD">
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo>∈</mo>
<mi>S</mi>
</mrow>
</munder>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo stretchy="false">)</mo>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle E[Y]=\sum _{y_{k}\in S}y_{k}f_{Y}(y_{k})}</annotation>
</semantics>
</math>
X と Y の分散は
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>V</mi>
<mo stretchy="false">[</mo>
<mi>X</mi>
<mo stretchy="false">]</mo>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo>−</mo>
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>X</mi>
<mo stretchy="false">]</mo>
<msup>
<mo stretchy="false">)</mo>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo stretchy="false">)</mo>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<msup>
<mi>x</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>x</mi>
<mo stretchy="false">)</mo>
<mi>d</mi>
<mi>x</mi>
<mo>−</mo>
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>X</mi>
<msup>
<mo stretchy="false">]</mo>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msup>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle V[X]=\int _{-\infty }^{+\infty }(x-E[X])^{2}f_{X}(x)dx=\int _{-\infty }^{+\infty }x^{2}f_{X}(x)dx-E[X]^{2}}</annotation>
</semantics>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<mi>V</mi>
<mo stretchy="false">[</mo>
<mi>Y</mi>
<mo stretchy="false">]</mo>
<mo>=</mo>
<munder>
<mo>∑</mo>
<mrow class="MJX-TeXAtom-ORD">
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo>∈</mo>
<mi>S</mi>
</mrow>
</munder>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo>−</mo>
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>Y</mi>
<mo stretchy="false">]</mo>
<msup>
<mo stretchy="false">)</mo>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo stretchy="false">)</mo>
<mo>=</mo>
<munder>
<mo>∑</mo>
<mrow class="MJX-TeXAtom-ORD">
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo>∈</mo>
<mi>S</mi>
</mrow>
</munder>
<msubsup>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<msub>
<mi>y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>k</mi>
</mrow>
</msub>
<mo stretchy="false">)</mo>
<mo>−</mo>
<mi>E</mi>
<mo stretchy="false">[</mo>
<mi>Y</mi>
<msup>
<mo stretchy="false">]</mo>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msup>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle V[Y]=\sum _{y_{k}\in S}(y_{k}-E[Y])^{2}f_{Y}(y_{k})=\sum _{y_{k}\in S}y_{k}^{2}f_{Y}(y_{k})-E[Y]^{2}}</annotation>
</semantics>
</math>
確率変数の変数変換による新しい変数の密度関数は、元の変数の密度関数で書くことができる。この公式は重積分における変数変換とほぼ同様である。
Rd から Rd への関数 T により、Rd 値確率変数 X と Y が
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <mi>X</mi> <mo>=</mo> <mi>T</mi> <mo stretchy="false">(</mo> <mi>Y</mi> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle X=T(Y)}</annotation> </semantics> </math>
と書けているとすると、Y の密度関数は X の密度関数を用いて
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>Y</mi> </mrow> </msub> <mo stretchy="false">(</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>d</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mo>=</mo> <mrow class="MJX-TeXAtom-ORD"> <mo stretchy="false">|</mo> </mrow> <mo stretchy="false">(</mo> <mo movablelimits="true" form="prefix">det</mo> <msub> <mi>J</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>T</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mo stretchy="false">(</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>d</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mrow class="MJX-TeXAtom-ORD"> <mo stretchy="false">|</mo> </mrow> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <mi>T</mi> <mo stretchy="false">(</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mo>…</mo> <mo>,</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>d</mi> </mrow> </msub> <mo stretchy="false">)</mo> <mo stretchy="false">)</mo> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle f_{Y}(y_{1},\ldots ,y_{d})=|(\det J_{T})(y_{1},\ldots ,y_{d})|f_{X}(T(y_{1},\ldots ,y_{d}))}</annotation> </semantics> </math>
となる。ただし J はヤコビアン とする。
例えばボックス-ミューラー変換は (0, 1]2 上の一様分布に従う確率変数 X = (X1, X2) を
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<msub>
<mi>Y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mrow class="MJX-TeXAtom-ORD">
<msqrt>
<mo>−</mo>
<mn>2</mn>
<mi>ln</mi>
<mo></mo>
<msub>
<mi>X</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>1</mn>
</mrow>
</msub>
</msqrt>
</mrow>
<mi>sin</mi>
<mo></mo>
<mo stretchy="false">(</mo>
<mn>2</mn>
<mi>π</mi>
<msub>
<mi>X</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msub>
<mo stretchy="false">)</mo>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle Y_{1}={\sqrt {-2\ln X_{1}}}\sin(2\pi X_{2})}</annotation>
</semantics>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<msub>
<mi>Y</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mrow class="MJX-TeXAtom-ORD">
<msqrt>
<mo>−</mo>
<mn>2</mn>
<mi>ln</mi>
<mo></mo>
<msub>
<mi>X</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>1</mn>
</mrow>
</msub>
</msqrt>
</mrow>
<mi>cos</mi>
<mo></mo>
<mo stretchy="false">(</mo>
<mn>2</mn>
<mi>π</mi>
<msub>
<mi>X</mi>
<mrow class="MJX-TeXAtom-ORD">
<mn>2</mn>
</mrow>
</msub>
<mo stretchy="false">)</mo>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle Y_{2}={\sqrt {-2\ln X_{1}}}\cos(2\pi X_{2})}</annotation>
</semantics>
</math>
によって変換する。X の密度関数は
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>X</mi> </mrow> </msub> <mo stretchy="false">(</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msub> <mo stretchy="false">)</mo> <mo>=</mo> <mrow class="MJX-TeXAtom-ORD"> <mrow> <mo>{</mo> <mtable columnalign="left left" rowspacing=".2em" columnspacing="1em" displaystyle="false"> <mtr> <mtd> <mn>1</mn> <mo>,</mo> </mtd> <mtd> <mo stretchy="false">(</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msub> <mo stretchy="false">)</mo> <mo>∈</mo> <mo stretchy="false">(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <msup> <mo stretchy="false">]</mo> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msup> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> <mo>,</mo> </mtd> <mtd> <mo stretchy="false">(</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>x</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msub> <mo stretchy="false">)</mo> <mo>∉</mo> <mo stretchy="false">(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <msup> <mo stretchy="false">]</mo> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msup> </mtd> </mtr> </mtable> <mo fence="true" stretchy="true"></mo> </mrow> </mrow> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle f_{X}(x_{1},x_{2})={\begin{cases}1,&(x_{1},x_{2})\in (0,1]^{2}\<add_contents_exp><m=6 date=20160913>,&(x_{1},x_{2})\notin (0,1]^{2}\end{cases}}}</annotation> </semantics> </math>
であり、上の公式を当てはめると Y の密度関数は
<math xmlns="http://www.w3.org/1998/Math/MathML" > <semantics> <mrow class="MJX-TeXAtom-ORD"> <mstyle displaystyle="true" scriptlevel="0"> <msub> <mi>f</mi> <mrow class="MJX-TeXAtom-ORD"> <mi>Y</mi> </mrow> </msub> <mo stretchy="false">(</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msub> <mo stretchy="false">)</mo> <mo>=</mo> <mrow class="MJX-TeXAtom-ORD"> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>π</mi> </mrow> </mfrac> </mrow> <mi>exp</mi> <mo></mo> <mrow> <mo>(</mo> <mo>−</mo> <mrow class="MJX-TeXAtom-ORD"> <mfrac> <mrow> <msubsup> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>1</mn> </mrow> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> <mrow class="MJX-TeXAtom-ORD"> <mn>2</mn> </mrow> </msubsup> </mrow> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> </mrow> </mstyle> </mrow> <annotation encoding="application/x-tex">{\displaystyle f_{Y}(y_{1},y_{2})={\frac {1}{2\pi }}\exp \left(-{\frac {y_{1}^{2}+y_{2}^{2}}{2}}\right)}</annotation> </semantics> </math>
となり、Y が二次元の標準正規分布に従うことが分かる。このように単純な分布を持つ変数を変換して、複雑な分布を作る操作は計算機による乱数の生成で重要となる。
2つの確率変数 X と Y の和 X + Y の確率分布や差 X - Y の確率分布は変数変換公式により計算できる。特に X と Y が独立で、確率密度関数がそれぞれ fX と fY だったとすると、和と差の確率密度関数は
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
<mo>+</mo>
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>t</mi>
<mo stretchy="false">)</mo>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>t</mi>
<mo>−</mo>
<mi>y</mi>
<mo stretchy="false">)</mo>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>y</mi>
<mo stretchy="false">)</mo>
<mspace width="thinmathspace" />
<mi>d</mi>
<mi>y</mi>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle f_{X+Y}(t)=\int _{-\infty }^{+\infty }f_{X}(t-y)f_{Y}(y)\,dy}</annotation>
</semantics>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" >
<semantics>
<mrow class="MJX-TeXAtom-ORD">
<mstyle displaystyle="true" scriptlevel="0">
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
<mo>−</mo>
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>t</mi>
<mo stretchy="false">)</mo>
<mo>=</mo>
<msubsup>
<mo>∫</mo>
<mrow class="MJX-TeXAtom-ORD">
<mo>−</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
<mrow class="MJX-TeXAtom-ORD">
<mo>+</mo>
<mi mathvariant="normal">∞</mi>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>X</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>y</mi>
<mo stretchy="false">)</mo>
<msub>
<mi>f</mi>
<mrow class="MJX-TeXAtom-ORD">
<mi>Y</mi>
</mrow>
</msub>
<mo stretchy="false">(</mo>
<mi>y</mi>
<mo stretchy="false">)</mo>
<mspace width="thinmathspace" />
<mi>d</mi>
<mi>y</mi>
</mstyle>
</mrow>
<annotation encoding="application/x-tex">{\displaystyle f_{X-Y}(t)=\int _{-\infty }^{+\infty }f_{X}(t+y)f_{Y}(y)\,dy}</annotation>
</semantics>
</math>
となる。
とくに和の確率密度関数は二つの分布の確率密度関数の畳み込みである。また、特性関数は確率密度関数のフーリエ変換であり、畳み込みのフーリエ変換は周波数領域における積であることから、和の特性関数は二つの分布の特性関数の積となる。
なお、確率変数の和の確率分布が元の分布族に従う場合、その分布は再生性があるという。
浸透(percolation) 確率に基づくモデル。具体的には森林火災の広がり、伝染病の伝搬、金属と絶縁体の混合物、強磁性元素と非磁性元素の混晶系、分子間の重合による巨大高分子のゲル化などがある[4]。
分岐過程(branching process)は、生命の数変化モデルです[5]。
酔歩(random walk)のモデル[6]。
無限粒子の遷移率の連続時間のモデル[7]。
拡散律速凝集(DLA: diffusion limited aggregation)と呼ぶ、ヴィッテンとサンダーによる粒子のクラスターが凝集によって成長するモデル。
バックたちによる砂山の斜面の崩壊を表すモデル。
交通流の渋滞モデル。
生命の時間的空間的モデル。セルオートマトンとも呼ぶ。生命競技(life game)は2次元セルオートマトンの一種である。
排他過程(exclusion porcess)は、連続時間で発展する確率モデル。上記生命モデルが離散時間の決定論的モデルであるのに対応している[8]。
|
|
</raw> </toggledisplay>
関連記事 | 「確率」「度」 |
.