コーシー分布

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:確率分布 コーシー分布(コーシーぶんぷ、テンプレート:Lang-en)は、連続確率分布の一種である。分布の名称は、フランスの数学者オーギュスタン=ルイ・コーシーに因む。確率密度関数は以下の式で与えられる。

f(x;x0,γ)=1πγ[1+(xx0γ)2]=1πγ(xx0)2+γ2

ここで テンプレート:Math は分布の最頻値を与えるテンプレート:Illテンプレート:Mvar半値半幅を与えるテンプレート:Illである。

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

テンプレート:Math2, テンプレート:Math2 である場合、この分布は標準コーシー分布と呼ばれ、以下の確率密度関数で与えられる。

f(x;0,1)=1π(1+x2)

性質

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

F(x;x0,γ)=1πarctan(xx0γ)+12

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

F1(p;x0,γ)=x0+γtan(π(p1/2)).

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

テンプレート:Mvar をコーシー分布に従う確率変数とする。コーシー分布の特性関数は以下のように与えられる。

ϕx(t;x0,γ)=E(eiXt)=exp(ix0tγ|t|).

テンプレート:Mvarテンプレート:Mvar を標準正規分布(期待値0、分散1の正規分布)に従う互いに独立な確率変数であるとすると、それらの比 テンプレート:Math は標準コーシー分布に従う。

テンプレート:Math2 をあるコーシー分布に従う独立な確率変数列とすると、それらの算術平均X=X1++Xnn は再び同じ位置母数、尺度母数を持つコーシー分布に従う(再生性)。この性質は、算術平均の特性関数が

ϕX(t)=E(eiX,t)=E(eiX1t/neiXnt/n)=E(eiX1t/n)E(eiXnt/n)=E(eiXt/n)n=exp(ix0t/nγ|t/n|)n=ϕX(t)

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

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

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

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

特性関数の求め方

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

1π(1+x2)

複素平面上で テンプレート:Math2 のみに1位の極を持つことを利用し、留数定理を用いて算出する。

期待値が定義されない理由

コーシー分布の期待値は、z=(xx0)/γと置換すると

E[x]=xf(x)dx=x0f(x)dx+(xx0)f(x)dx=x0+γπxx0(xx0)2+γ2dx=x0+γπlimR1,R2R1R2z1+z2dz=x0+γ2πlimR1,R2[log(1+z2)]R1R2=x0+γ2πlimR1,R2log(1+R221+R12)

となるが、この広義積分の値は存在せず(十分大きな テンプレート:Math2 について、log((1+R22)/(1+R12))(,) のどのような値でも取りうるので、二重極限としての収束値は存在しない)、このため期待値は存在しない。なお、コーシーの主値 limRRRxf(x)dxテンプレート:Math である。

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

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

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

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

E(X2)x21+x2dx=dx11+x2dx=π=.

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

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

関連項目

テンプレート:確率分布の一覧