コーシー分布 – Wikipedia

Posted on by

コーシー分布（コーシーぶんぷ、英語: Cauchy distribution）は、連続確率分布の一種である。分布の名称は、フランスの数学者オーギュスタン＝ルイ・コーシーに因む。確率密度関数は以下の式で与えられる。

f(x;x0,γ)=1πγ[1+(x−x0γ)2]=1πγ(x−x0)2+γ2{displaystyle f(x;x_{0},gamma )={frac {1}{pi gamma left[1+left({dfrac {x-x_{0}}{gamma }}right)^{2}right]}}={frac {1}{pi }}{frac {gamma }{(x-x_{0})^{2}+gamma ^{2}}}}

ここで x0 は分布の最頻値を与える位置母数、γ は半値半幅を与える尺度母数である。

この分布は、ヘンドリック・ローレンツの名を取ってローレンツ分布と呼ばれることもあり、またこれら2人の名前を合わせてコーシー-ローレンツ分布とも呼ばれる。また物理学の分野では、ブライト・ウィグナー分布という名前で知られている。この分布は強制共鳴を記述する微分方程式の解となることから、物理学では重要な存在となっている。また分光学では共鳴広がりを含む多くのメカニズムによって広げられたスペクトル線の形状を記述するために用いられる。以下では、統計学における名称であるコーシー分布を用いて説明する。

x0 = 0, γ = 1 である場合、この分布は標準コーシー分布と呼ばれ、以下の確率密度関数で与えられる。

f(x;0,1)=1π(1+x2){displaystyle f(x;0,1)={frac {1}{pi (1+x^{2})}}}

累積分布関数は以下のようになる。

F(x;x0,γ)=1πarctan⁡(x−x0γ)+12{displaystyle F(x;x_{0},gamma )={frac {1}{pi }}arctan left({frac {x-x_{0}}{gamma }}right)+{frac {1}{2}}}

また、逆累積分布関数は次の通りである。

F−1(p;x0,γ)=x0+γtan⁡(π(p−1/2)).{displaystyle F^{-1}(p;x_{0},gamma )=x_{0}+gamma tan(pi (p-1/2)).}

コーシー分布は、期待値や分散（およびより高次のモーメント）が定義されない分布の例として知られる。最頻値と中央値は常に定義され、それらはいずれも x0 で与えられる。

X をコーシー分布に従う確率変数とする。コーシー分布の特性関数は以下のように与えられる。

ϕx(t;x0,γ)=E(eiXt)=exp⁡(ix0t−γ|t|).{displaystyle phi _{x}(t;x_{0},gamma )=mathrm {E} (e^{iXt})=exp(ix_{0}t-gamma |t|).}

U と V を標準正規分布（期待値0、分散1の正規分布）に従う互いに独立な確率変数であるとすると、それらの比 U/V は標準コーシー分布に従う。

X1, X2, …, Xn をあるコーシー分布に従う独立な確率変数列とすると、それらの算術平均

X¯=X1+⋯+Xnn{displaystyle {overline {X}}={frac {X_{1}+dotsb +X_{n}}{n}}}

は再び同じ位置母数、尺度母数を持つコーシー分布に従う（再生性）。この性質は、算術平均の特性関数が

ϕX¯(t)=E⁡(eiX¯,t)=E⁡(eiX1t/n⋯eiXnt/n)=E⁡(eiX1t/n)⋯E⁡(eiXnt/n)=E⁡(eiXt/n)n=exp⁡(ix0t/n−γ|t/n|)n=ϕX(t){displaystyle {begin{aligned}phi _{overline {X}}(t)&=operatorname {E} left(e^{i{overline {X}},t}right)=operatorname {E} left(e^{iX_{1}t/n}dotsb e^{iX_{n}t/n}right)\&=operatorname {E} left(e^{iX_{1}t/n}right)dotsb operatorname {E} left(e^{iX_{n}t/n}right)=operatorname {E} left(e^{iXt/n}right)^{n}\&=exp(ix_{0}t/n-gamma |t/n|)^{n}=phi _{X}(t)end{aligned}}}

となることから導かれる。このように、コーシー分布に従う確率変数を増やしても、その算術平均の分布は正規分布に近づかないため、中心極限定理における有限分散の仮定は必須であることが分かる。また、これは安定分布族（コーシー分布は安定分布族に含まれる）における一般化中心極限定理の例でもある。

コーシー分布は無限分解可能な分布である。

自由度1のT分布は、標準コーシー分布と一致する。

コーシー分布が属している位置尺度母数分布族は、実係数メビウス変換に関して閉じている。

特性関数の求め方[編集]

コーシー分布の特性関数の求め方は、標準コーシー分布の確率密度関数

1π(1+x2){displaystyle {frac {1}{pi (1+x^{2})}}}

が複素平面上で x = ±i のみに1位の極を持つことを利用し、留数定理を用いて算出する。

期待値が定義されない理由[編集]

コーシー分布の期待値は、

z=1γ(x−x0){displaystyle z={frac {1}{gamma }}(x-x_{0})}

と置換すると

E⁡[x]=∫−∞∞xf(x)dx=∫−∞∞x0f(x)dx+∫−∞∞(x−x0)f(x)dx=x0+γπ∫−∞∞x−x0(x−x0)2+γ2dx=x0+limR1,R2→∞γπ∫−R1R2z1+z2dz=x0+limR1,R2→∞γ2π[log⁡(1+z2)]−R1R2=x0+limR1,R2→∞γ2πlog⁡(1+R221+R12){displaystyle {begin{aligned}operatorname {E} [x]&=int _{-infty }^{infty }xf(x),dx=int _{-infty }^{infty }x_{0}f(x),dx+int _{-infty }^{infty }(x-x_{0})f(x),dx\&=x_{0}+{frac {gamma }{pi }}int _{-infty }^{infty }{frac {x-x_{0}}{(x-x_{0})^{2}+gamma ^{2}}},dx\&=x_{0}+lim _{R_{1},R_{2}to infty }{frac {gamma }{pi }}int _{-R_{1}}^{R_{2}}{frac {z}{1+z^{2}}},dz=x_{0}+lim _{R_{1},R_{2}to infty }{frac {gamma }{2pi }}left[log(1+z^{2})right]_{-R_{1}}^{R_{2}}\&=x_{0}+lim _{R_{1},R_{2}to infty }{frac {gamma }{2pi }}log left({frac {1+{R_{2}}^{2}}{1+{R_{1}}^{2}}}right)end{aligned}}}

となるが、この広義積分の値は存在せず（十分大きな R1, R2 について、

log⁡((1+R22)/(1+R12)){displaystyle log left((1+{R_{2}}^{2})/(1+{R_{1}}^{2})right)}

は

(−∞,∞){displaystyle (-infty ,infty )}

のどのような値でも取りうるので、二重極限としての収束値は存在しない）、このため期待値は存在しない。なお、コーシーの主値

limR→∞∫−RRxf(x)dx{displaystyle lim _{Rto infty }int _{-R}^{R}xf(x),dx}

は x0 である。

大数の強法則など、期待値に関する確率論のさまざまな結果は、このようなケースでは成立しない。

また、コーシー分布に従う母集団から無作為抽出された標本に関する算術平均は、ただ一つの抽出による結果からは一切改善されない。これは、標本に極端に大きな（あるいは小さな）値が含まれる可能性がかなり高いからである。しかし、標本中央値（これは極端な値には影響を受けない）は中心（最頻値）を知るための一つの尺度となりうる。

2次モーメントが無限大になる理由[編集]

期待値が定義されない限り、分散や標準偏差を考えることは不可能である。しかし、原点を中心とした2次モーメントを考えることは可能である。しかし、これもまた無限大となる。

E⁡(X2)∝∫−∞∞x21+x2dx=∫−∞∞dx−∫−∞∞11+x2dx=∞−π=∞.{displaystyle operatorname {E} (X^{2})propto int _{-infty }^{infty }{frac {x^{2}}{1+x^{2}}},dx=int _{-infty }^{infty }dx-int _{-infty }^{infty }{frac {1}{1+x^{2}}},dx=infty -pi =infty .}

相対論的ブライト・ウィグナー分布[編集]

原子核物理学および素粒子物理学において、共鳴のエネルギー特性は相対論的ブライト・ウィグナー分布によって記述される。

関連項目[編集]