確率分布
コーシー 分布は、 オーギュスタン=ルイ・コーシー にちなんで名付けられた 連続確率分布 です 。特に 物理学者 の間では、 ローレンツ分布 ( ヘンドリック・ローレンツ にちなんで)、 コーシー・ローレンツ分布 、 ローレンツ関数 、または ブライト・ウィグナー分布 としても知られています。コーシー分布は、一様分布の角度で 発せられる光線の x 切片の分布です。また、平均が 0 である 2 つの独立した 正規分布するランダム変数の 比 の分布でもあります 。
ふ
(
x
;
x
0
、
γ
)
{\displaystyle f(x;x_{0},\gamma )}
(
x
0
、
γ
)
{\displaystyle (x_{0},\gamma )}
コーシー分布は、その 期待値 と 分散 が定義されていないため、統計学では「 異常な 」分布の標準的な例としてよく使用されます(ただし、以下の § モーメントを参照してください)。コーシー分布には、1 以上の有限 モーメント はなく、分数絶対モーメントのみが存在します。 [1] コーシー分布には モーメント生成関数 がありません。
数学 においては、 上半平面 における ラプラス方程式 の 基本解 である ポアソン核 と密接な関係がある 。
これは、解析的に表現できる確率密度関数を持つ数少ない 安定した分布 の 1 つであり、他の 2 つは 正規分布 と レヴィ分布 です。
定義
ここに最も重要な構造があります。
回転対称性
線の前に立って、線に向かって一様にランダムに方向(より正確には角度)でボールを蹴った場合、ボールが線に当たる点の分布はコーシー分布になります。
より正式には、xy 平面の点 を考え 、その点を通る直線を選択します。直線の方向 ( - 軸との角度) は、一様 (-90° から +90° の間) にランダムに選択されます。直線と x 軸の交点は、位置 とスケールの Cauchy 分布です 。
(
x
0
、
γ
)
{\displaystyle (x_{0},\gamma )}
x
{\displaystyle x}
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma}
この定義は、標準コーシー分布からサンプルを採取する簡単な方法を提供します。 を からの一様分布からのサンプルとすると、 を 使って標準コーシー分布から
サンプル を生成できます。
あなた
{\displaystyle u}
[
0
、
1
]
{\displaystyle [0,1]}
x
{\displaystyle x}
x
=
日焼け
(
π
(
あなた
−
1
2
)
)
{\displaystyle x=\tan \left(\pi (u-{\frac {1}{2}})\right)}
とが、 期待値 0 および 分散 1の 2 つの独立した 正規分布の ランダム変数 である場合 、比は 標準コーシー分布に従います。
あなた
{\displaystyle U}
五
{\displaystyle V}
あなた
/
五
{\displaystyle U/V}
より一般的には、が 平面上の回転対称分布である場合、比は 標準コーシー分布に従います。
(
あなた
、
五
)
{\displaystyle (U,V)}
あなた
/
五
{\displaystyle U/V}
確率密度関数 (PDF)
コーシー分布は、次の確率密度関数 (PDF) [1] [2] を持つ確率分布である。
ふ
(
x
;
x
0
、
γ
)
=
1
π
γ
[
1
+
(
x
−
x
0
γ
)
2
]
=
1
π
[
γ
(
x
−
x
0
)
2
+
γ
2
]
、
{\displaystyle f(x;x_{0},\gamma )={\frac {1}{\pi \gamma \left[1+\left({\frac {x-x_{0}}{\gamma }}\right)^{2}\right]}}={1 \over \pi }\left[{\gamma \over (x-x_{0})^{2}+\gamma ^{2}}\right],}
ここで 、 は 分布のピークの位置を指定する 位置パラメータ 、は半値幅 (HWHM) を指定する スケールパラメータ、あるいは 半値全幅 (FWHM)です 。 は 四分位範囲の 半分に等しく、 確率誤差 と呼ばれることもあります。 この関数は ローレンツ関数 [3] としても知られ 、 新生デルタ関数 の例でもあり 、したがって の 極限で ディラックのデルタ関数 に近づきます。 オーギュスタン=ルイ・コーシーは1827 年にこのような密度関数を 無限 小スケールパラメータで利用し 、この ディラックのデルタ関数 を定義しました。
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma}
2
γ
{\displaystyle 2\gamma }
γ
{\displaystyle \gamma}
γ
→
0
{\displaystyle \gamma \to 0}
PDFのプロパティ
コーシー PDF の最大値または振幅は で 、 に位置します 。
1
π
γ
{\displaystyle {\frac {1}{\pi \gamma }}}
x
=
x
0
{\displaystyle x=x_{0}}
PDFを複素パラメータで表現すると便利な場合がある。
ψ
=
x
0
+
私
γ
{\displaystyle \psi =x_{0}+i\gamma }
ふ
(
x
;
ψ
)
=
1
π
私は
(
1
x
−
ψ
)
=
1
π
再
(
−
私
x
−
ψ
)
{\displaystyle f(x;\psi )={\frac {1}{\pi }}\,{\textrm {Im}}\left({\frac {1}{x-\psi }}\right)={\frac {1}{\pi }}\,{\textrm {Re}}\left({\frac {-i}{x-\psi }}\right)}
の特別な場合 、 確率密度関数 [4] [5]を持つ 標準コーシー分布 と呼ばれる。
x
0
=
0
{\displaystyle x_{0}=0}
γ
=
1
{\displaystyle \gamma =1}
ふ
(
x
;
0
、
1
)
=
1
π
(
1
+
x
2
)
。
{\displaystyle f(x;0,1)={\frac {1}{\pi (1+x^{2})}}.\!}
物理学では、3 パラメータのローレンツ関数がよく使用されます。
ふ
(
x
;
x
0
、
γ
、
私
)
=
私
[
1
+
(
x
−
x
0
γ
)
2
]
=
私
[
γ
2
(
x
−
x
0
)
2
+
γ
2
]
、
{\displaystyle f(x;x_{0},\gamma ,I)={\frac {I}{\left[1+\left({\frac {x-x_{0}}{\gamma }}\right)^{2}\right]}}=I\left[{\gamma ^{2} \over (x-x_{0})^{2}+\gamma ^{2}}\right],}
ここで、 ピークの高さです。示されている3パラメータのローレンツ関数は、一般には確率密度関数ではありません。なぜなら、次の特別な場合を除いて、1に積分されないからです。
私
{\displaystyle I}
私
=
1
π
γ
。
{\displaystyle I={\frac {1}{\pi \gamma }}.\!}
累積分布関数 (CDF)
コーシー分布は、次の 累積分布関数 (CDF) を持つ確率分布です。
ふ
(
x
;
x
0
、
γ
)
=
1
π
アークタンジェント
(
x
−
x
0
γ
)
+
1
2
{\displaystyle F(x;x_{0},\gamma )={\frac {1}{\pi }}\arctan \left({\frac {x-x_{0}}{\gamma }}\right)+{\frac {1}{2}}}
コーシー分布の
分位関数 ( 逆累積 分布関数)は
質問
(
p
;
x
0
、
γ
)
=
x
0
+
γ
日焼け
[
π
(
p
−
1
2
)
]
。
{\displaystyle Q(p;x_{0},\gamma )=x_{0}+\gamma \,\tan \left[\pi \left(p-{\tfrac {1}{2}}\right)\right].}
したがって、第 1 四分位数と第 3 四分位数は となり 、 四分位範囲 は となります 。
(
x
0
−
γ
、
x
0
+
γ
)
{\displaystyle (x_{0}-\gamma ,x_{0}+\gamma )}
2
γ
{\displaystyle 2\gamma }
標準分布の場合、累積分布関数は 逆正接関数 に簡略化されます。
アークタンジェント
(
x
)
{\displaystyle \arctan(x)}
ふ
(
x
;
0
、
1
)
=
1
π
アークタンジェント
(
x
)
+
1
2
{\displaystyle F(x;0,1)={\frac {1}{\pi }}\arctan \left(x\right)+{\frac {1}{2}}}
その他の構造
標準的なコーシー分布は自由度1の スチューデントの t 分布 であり、スチューデントのt分布を構築する任意の方法で構築できます。 [6]
が正の対角成分を持つ正半正定値共分散行列で ある 場合、 独立かつ同一分布する 、 および および から独立でかつ ( カテゴリ分布 を 定義する)任意のランダム-ベクトルに対して 、 次が成り立つ。
Σ
{\displaystyle \Sigma }
p
×
p
{\displaystyle p\times p}
バツ
、
はい
〜
いいえ
(
0
、
Σ
)
{\displaystyle X,Y\sim N(0,\Sigma )}
p
{\displaystyle p}
わ
{\displaystyle w}
バツ
{\displaystyle X}
はい
{\displaystyle Y}
わ
1
+
⋯
+
わ
p
=
1
{\displaystyle w_{1}+\cdots +w_{p}=1}
わ
私
≥
0
、
私
=
1
、
…
、
p
、
{\displaystyle w_{i}\geq 0,i=1,\ldots ,p,}
∑
じゅう
=
1
p
わ
じゅう
バツ
じゅう
はい
じゅう
〜
C
1つの
あなた
c
h
ええ
(
0
、
1
)
。
{\displaystyle \sum _{j=1}^{p}w_{j}{\frac {X_{j}}{Y_{j}}}\sim \mathrm {コーシー} (0,1)。}
[7]
プロパティ
コーシー分布は、平均 、 分散 、または高次の モーメント が定義されていない分布の例です 。その モード と 中央値は 明確に定義されており、どちらも に等しくなります 。
x
0
{\displaystyle x_{0}}
コーシー分布は 無限に分割可能な確率分布である。また、厳密に 安定な 分布でもある 。 [8]
すべての安定分布と同様に、 コーシー分布が属する 位置スケール族は、 実係数を持つ 線形変換 に対して閉じています。 さらに、コーシー分布のランダム変数の族は、実係数を持つ 線形分数変換 に対して閉じています。 [9] これに関連して、 McCullaghのコーシー分布のパラメータ化 も参照してください。
コーシー分布に従う確率変数の合計
標準コーシー分布からのIID サンプルの 場合 、その サンプル平均 も標準コーシー分布になります。特に、平均は平均に収束しないため、標準コーシー分布は大数の法則に従いません。
バツ
1
、
バツ
2
、
…
、
バツ
ん
{\displaystyle X_{1},X_{2},\ldots ,X_{n}}
バツ
¯
=
1
ん
∑
私
バツ
私
{\displaystyle {\bar {X}}={\frac {1}{n}}\sum _{i}X_{i}}
これは、PDF との繰り返し積分によって証明できますが、より便利な方法として、 標準コーシー分布の 特性関数 を使用することで証明できます (以下を参照)。これにより、 が得られ 、 標準コーシー分布も得られます。
φ
X
(
t
)
=
E
[
e
i
X
t
]
=
e
−
|
t
|
.
{\displaystyle \varphi _{X}(t)=\operatorname {E} \left[e^{iXt}\right]=e^{-|t|}.}
φ
∑
i
X
i
(
t
)
=
e
−
n
|
t
|
{\displaystyle \varphi _{\sum _{i}X_{i}}(t)=e^{-n|t|}}
X
¯
{\displaystyle {\bar {X}}}
より一般的には、 が独立で、位置パラメータ とスケール を持つコーシー分布であり 、が 実数である場合、 は 位置 とスケールを持つコーシー分布 です。独立したコーシー分布の加重和には大数の法則がないことがわかります。
X
1
,
X
2
,
…
,
X
n
{\displaystyle X_{1},X_{2},\ldots ,X_{n}}
x
1
,
…
,
x
n
{\displaystyle x_{1},\ldots ,x_{n}}
γ
1
,
…
,
γ
n
{\displaystyle \gamma _{1},\ldots ,\gamma _{n}}
a
1
,
…
,
a
n
{\displaystyle a_{1},\ldots ,a_{n}}
∑
i
a
i
X
i
{\displaystyle \sum _{i}a_{i}X_{i}}
∑
i
a
i
x
i
{\displaystyle \sum _{i}a_{i}x_{i}}
∑
i
|
a
i
|
γ
i
{\displaystyle \sum _{i}|a_{i}|\gamma _{i}}
これは、中心極限定理 における有限分散の条件は省略できないことを示しています。これは、すべての 安定分布 に特徴的な中心極限定理のより一般化されたバージョンの例でもあり 、コーシー分布はその特殊なケースです。
中心極限定理
が有限だがゼロでない PDFを持つIIDサンプルである 場合 、 分布 はスケール の Cauchy 分布に収束します 。 [10]
X
1
,
X
2
,
…
{\displaystyle X_{1},X_{2},\ldots }
ρ
{\displaystyle \rho }
lim
c
→
∞
1
c
∫
−
c
c
x
2
ρ
(
x
)
d
x
=
2
γ
π
{\displaystyle \lim _{c\to \infty }{\frac {1}{c}}\int _{-c}^{c}x^{2}\rho (x)\,dx={\frac {2\gamma }{\pi }}}
1
n
∑
i
=
1
n
X
i
{\displaystyle {\frac {1}{n}}\sum _{i=1}^{n}X_{i}}
γ
{\displaystyle \gamma }
特性関数
をコーシー分布の確率変数とします。 コーシー 分布の
特性関数は次のように与えられます。
X
{\displaystyle X}
φ
X
(
t
)
=
E
[
e
i
X
t
]
=
∫
−
∞
∞
f
(
x
;
x
0
,
γ
)
e
i
x
t
d
x
=
e
i
x
0
t
−
γ
|
t
|
.
{\displaystyle \varphi _{X}(t)=\operatorname {E} \left[e^{iXt}\right]=\int _{-\infty }^{\infty }f(x;x_{0},\gamma )e^{ixt}\,dx=e^{ix_{0}t-\gamma |t|}.}
これは確率密度のフーリエ変換 に過ぎません 。元の確率密度は、基本的に逆フーリエ変換を使用して特性関数で表現できます。
f
(
x
;
x
0
,
γ
)
=
1
2
π
∫
−
∞
∞
φ
X
(
t
;
x
0
,
γ
)
e
−
i
x
t
d
t
{\displaystyle f(x;x_{0},\gamma )={\frac {1}{2\pi }}\int _{-\infty }^{\infty }\varphi _{X}(t;x_{0},\gamma )e^{-ixt}\,dt\!}
分布の n 次モーメントは、特性関数の n 次導関数を で評価したものです 。 特性 関数 は 原点で 微分可能で はないことに注意してください。これは、コーシー分布がゼロ次モーメントよりも高い明確に定義されたモーメントを持たないという事実に対応しています。
t
=
0
{\displaystyle t=0}
カルバック・ライブラー距離
2つのコーシー分布間のカルバック ・ライブラー距離は、 次のような対称閉形式の式で表される: [11]
K
L
(
p
x
0
,
1
,
γ
1
:
p
x
0
,
2
,
γ
2
)
=
log
(
γ
1
+
γ
2
)
2
+
(
x
0
,
1
−
x
0
,
2
)
2
4
γ
1
γ
2
.
{\displaystyle \mathrm {KL} \left(p_{x_{0,1},\gamma _{1}}:p_{x_{0,2},\gamma _{2}}\right)=\log {\frac {\left(\gamma _{1}+\gamma _{2}\right)^{2}+\left(x_{0,1}-x_{0,2}\right)^{2}}{4\gamma _{1}\gamma _{2}}}.}
2つのコーシー分布間のfダイバージェンス は 対称であり、カイ2乗ダイバージェンスの関数として表すことができます。 [12] 全変動 、 ジェンセン・シャノン・ダイバージェンス 、 ヘリンガー距離 など
の閉じた形式の表現 が利用可能です。
エントロピ
コーシー分布のエントロピーは次のように表されます。
H
(
γ
)
=
−
∫
−
∞
∞
f
(
x
;
x
0
,
γ
)
log
(
f
(
x
;
x
0
,
γ
)
)
d
x
=
log
(
4
π
γ
)
{\displaystyle {\begin{aligned}H(\gamma )&=-\int _{-\infty }^{\infty }f(x;x_{0},\gamma )\log(f(x;x_{0},\gamma ))\,dx\\[6pt]&=\log(4\pi \gamma )\end{aligned}}}
コーシー分布の
分位関数( 分位密度関数)の導関数は次のようになります。
Q
′
(
p
;
γ
)
=
γ
π
sec
2
[
π
(
p
−
1
2
)
]
.
{\displaystyle Q'(p;\gamma )=\gamma \,\pi \,{\sec }^{2}\left[\pi \left(p-{\tfrac {1}{2}}\right)\right].\!}
分布の微分エントロピーは 、 その分位密度によって定義することができる。 [13] 具体的には、
H
(
γ
)
=
∫
0
1
log
(
Q
′
(
p
;
γ
)
)
d
p
=
log
(
4
π
γ
)
{\displaystyle H(\gamma )=\int _{0}^{1}\log \,(Q'(p;\gamma ))\,\mathrm {d} p=\log(4\pi \gamma )}
コーシー分布は、 ランダム変量 に対する 最大エントロピー確率分布であり、 [14]
X
{\displaystyle X}
E
[
log
(
1
+
(
X
−
x
0
)
2
/
γ
2
)
]
=
log
4
{\displaystyle \operatorname {E} [\log(1+(X-x_{0})^{2}/\gamma ^{2})]=\log 4}
瞬間
コーシー分布は、明確に定義された(または「不定」な)モーメントを持たない分布として、初等確率の授業で説明的な反例としてよく使用されます。
サンプルの瞬間
標準コーシー分布からIID サンプルを取ると 、そのサンプル平均のシーケンスは となり、これも標準コーシー分布に従います。したがって、項をいくつ取っても、サンプル平均は収束しません。
X
1
,
X
2
,
…
{\displaystyle X_{1},X_{2},\ldots }
S
n
=
1
n
∑
i
=
1
n
X
i
{\displaystyle S_{n}={\frac {1}{n}}\sum _{i=1}^{n}X_{i}}
同様に、標本分散 も収束しません。
V
n
=
1
n
∑
i
=
1
n
(
X
i
−
S
n
)
2
{\displaystyle V_{n}={\frac {1}{n}}\sum _{i=1}^{n}(X_{i}-S_{n})^{2}}
サンプル平均の典型的な軌跡は、ゼロへのゆっくりとした収束の長い期間のように見え、ゼロからの大きなジャンプで区切られますが、決して離れすぎることはありません。サンプル分散の典型的な軌跡も同様に見えますが、ジャンプは減衰よりも速く蓄積され、無限大に発散します。
の典型的な軌道は、 ゼロへのゆっくりとした収束の長い期間のように見えますが、ゼロから大きく離れるジャンプが中断されますが、離れすぎることはあり ません。 の典型的な軌道も同様に見えますが、ジャンプは減衰よりも速く蓄積され、無限大に発散します。これらの 2 種類の軌道が図にプロットされています。
S
1
,
S
2
,
.
.
.
{\displaystyle S_{1},S_{2},...}
V
1
,
V
2
,
.
.
.
{\displaystyle V_{1},V_{2},...}
1 次より低いサンプルのモーメントはゼロに収束します。2 次より高いサンプルのモーメントは、サンプル分散よりもさらに速く無限大に発散します。
平均
確率分布が 密度関数 を持つ 場合 、平均は、存在するならば、次のように与えられる。
f
(
x
)
{\displaystyle f(x)}
この両側 仮積分は、 2つの片側仮積分の和を計算することで評価できる。つまり、
任意の実数に対して 。
a
{\displaystyle a}
積分が存在するためには(たとえ無限大であっても)、この和の項のうち少なくとも1つは有限であるか、または両方の項が無限大で同じ符号を持っている必要があります。しかし、コーシー分布の場合、この和( 2 )の項は両方とも無限大で、反対の符号を持っています。したがって、( 1 )は未定義であり、平均も未定義です。 [15] 確率分布関数(PDF)の平均が未定義の場合、サンプルのサイズに関係なく、誰も実験データポイントの信頼できる平均を計算することはできません。
コーシー分布の平均の
コーシー主値 はゼロであることに注意する
。一方、積分を計算するとわかるように、関連する積分は
ゼロで
は ない。これもまた、平均( 1 )が存在できないことを示している 。
lim
a
→
∞
∫
−
a
a
x
f
(
x
)
d
x
{\displaystyle \lim _{a\to \infty }\int _{-a}^{a}xf(x)\,dx}
lim
a
→
∞
∫
−
2
a
a
x
f
(
x
)
d
x
{\displaystyle \lim _{a\to \infty }\int _{-2a}^{a}xf(x)\,dx}
確率論における期待値 に関する様々な結果 、例えば 大数の強い法則は コーシー分布には当てはまらない。 [15]
小さな瞬間
の絶対モーメントは定義されて い
ます。
p
∈
(
−
1
,
1
)
{\displaystyle p\in (-1,1)}
X
∼
C
a
u
c
h
y
(
0
,
γ
)
{\displaystyle X\sim \mathrm {Cauchy} (0,\gamma )}
E
[
|
X
|
p
]
=
γ
p
s
e
c
(
π
p
/
2
)
.
{\displaystyle \operatorname {E} [|X|^{p}]=\gamma ^{p}\mathrm {sec} (\pi p/2).}
より高い瞬間
コーシー分布には、いかなる次数の有限モーメントもありません。より高い生 のモーメント の中には、無限大の値を持つものもいくつか存在します。たとえば、生の 2 番目のモーメントです。
E
[
X
2
]
∝
∫
−
∞
∞
x
2
1
+
x
2
d
x
=
∫
−
∞
∞
1
−
1
1
+
x
2
d
x
=
∫
−
∞
∞
d
x
−
∫
−
∞
∞
1
1
+
x
2
d
x
=
∫
−
∞
∞
d
x
−
π
=
∞
.
{\displaystyle {\begin{aligned}\operatorname {E} [X^{2}]&\propto \int _{-\infty }^{\infty }{\frac {x^{2}}{1+x^{2}}}\,dx=\int _{-\infty }^{\infty }1-{\frac {1}{1+x^{2}}}\,dx\\[8pt]&=\int _{-\infty }^{\infty }dx-\int _{-\infty }^{\infty }{\frac {1}{1+x^{2}}}\,dx=\int _{-\infty }^{\infty }dx-\pi =\infty .\end{aligned}}}
式を整理すると、第 2 モーメントは本質的に定数 (ここでは 1) の無限積分であることがわかります。偶数乗の高い生のモーメントも無限大になります。ただし、奇数乗の生のモーメントは未定義であり、無限大の値で存在する場合とは明確に異なります。奇数乗の生のモーメントが未定義なのは、 積分の 2 つの半分が両方とも発散し、符号が反対であるため、その値が本質的に に等しいためです。第 1 の生のモーメントは平均であり、奇数であるため存在しません (これについては上記の説明も参照してください)。これは、 中心モーメント と 標準化モーメント がすべて平均に基づいているため、定義されていないことを意味します。分散 (第 2 の中心モーメント) も同様に存在しません (生の第 2 モーメントが値無限大で存在するにもかかわらず)。
∞
−
∞
{\displaystyle \infty -\infty }
高次のモーメントの結果は、 ヘルダーの不等式 から得られます。これは、低次のモーメントが発散すると、高次のモーメント(またはモーメントの半分)も発散することを意味します。
切断分布のモーメント
標準的なコーシー分布を区間 [−10 100 , 10 100 ]に制限して定義される 切断分布 を考えてみましょう 。このような切断分布はすべてのモーメントを持ちます(そして中心極限定理はそこからの iid 観測に適用されます)。しかし、ほとんどすべての実用目的では、それはコーシー分布のように動作します。 [16]
もし そうなら [17]
X
∼
Cauchy
(
x
0
,
γ
)
{\displaystyle X\sim \operatorname {Cauchy} (x_{0},\gamma )}
k
X
+
ℓ
∼
Cauchy
(
x
0
k
+
ℓ
,
γ
|
k
|
)
{\displaystyle kX+\ell \sim {\textrm {Cauchy}}(x_{0}k+\ell ,\gamma |k|)}
とが独立であれ ば 、 と
X
∼
Cauchy
(
x
0
,
γ
0
)
{\displaystyle X\sim \operatorname {Cauchy} (x_{0},\gamma _{0})}
Y
∼
Cauchy
(
x
1
,
γ
1
)
{\displaystyle Y\sim \operatorname {Cauchy} (x_{1},\gamma _{1})}
X
+
Y
∼
Cauchy
(
x
0
+
x
1
,
γ
0
+
γ
1
)
{\displaystyle X+Y\sim \operatorname {Cauchy} (x_{0}+x_{1},\gamma _{0}+\gamma _{1})}
X
−
Y
∼
Cauchy
(
x
0
−
x
1
,
γ
0
+
γ
1
)
{\displaystyle X-Y\sim \operatorname {Cauchy} (x_{0}-x_{1},\gamma _{0}+\gamma _{1})}
もし そうなら
X
∼
Cauchy
(
0
,
γ
)
{\displaystyle X\sim \operatorname {Cauchy} (0,\gamma )}
1
X
∼
Cauchy
(
0
,
1
γ
)
{\displaystyle {\tfrac {1}{X}}\sim \operatorname {Cauchy} (0,{\tfrac {1}{\gamma }})}
McCullaghのコーシー分布のパラメータ化 : [18] コーシー分布を 1つの複素パラメータ で表すと 、 は を意味すると定義します 。 のとき 、次のようになります。 ここで 、、、 は 実数です。
ψ
=
x
0
+
i
γ
{\displaystyle \psi =x_{0}+i\gamma }
X
∼
Cauchy
(
ψ
)
{\displaystyle X\sim \operatorname {Cauchy} (\psi )}
X
∼
Cauchy
(
x
0
,
|
γ
|
)
{\displaystyle X\sim \operatorname {Cauchy} (x_{0},|\gamma |)}
X
∼
Cauchy
(
ψ
)
{\displaystyle X\sim \operatorname {Cauchy} (\psi )}
a
X
+
b
c
X
+
d
∼
Cauchy
(
a
ψ
+
b
c
ψ
+
d
)
{\displaystyle {\frac {aX+b}{cX+d}}\sim \operatorname {Cauchy} \left({\frac {a\psi +b}{c\psi +d}}\right)}
a
{\displaystyle a}
b
{\displaystyle b}
c
{\displaystyle c}
d
{\displaystyle d}
上記と同じ規則を用いると、 次のようになる。 [18] ここで、 円コーシー分布 である 。
X
∼
Cauchy
(
ψ
)
{\displaystyle X\sim \operatorname {Cauchy} (\psi )}
X
−
i
X
+
i
∼
CCauchy
(
ψ
−
i
ψ
+
i
)
{\displaystyle {\frac {X-i}{X+i}}\sim \operatorname {CCauchy} \left({\frac {\psi -i}{\psi +i}}\right)}
CCauchy
{\displaystyle \operatorname {CCauchy} }
統計的推論
パラメータの推定
コーシー分布のパラメータは平均と分散に対応していないため、標本平均と標本分散を使用してコーシー分布のパラメータを推定しようとすると成功しません。 [19] たとえば、コーシー分布からサイズ n のiid標本を抽出した場合、標本平均は次のように計算できます。
x
¯
=
1
n
∑
i
=
1
n
x
i
{\displaystyle {\bar {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}}
標本値は 中心値 の周囲に集中しますが 、観測値が増えるにつれて、絶対値の大きい標本点に遭遇する確率が高くなるため、標本平均はますます変動しやすくなります。実際、標本平均の分布は、観測値自体の分布に等しくなります。つまり、大きな標本の標本平均は、 標本からの単一の観測値よりも の推定値として優れている (または劣っている) わけではありません。同様に、標本分散を計算すると、観測値が増えるにつれて値が大きくなります。
x
i
{\displaystyle x_{i}}
x
0
{\displaystyle x_{0}}
x
0
{\displaystyle x_{0}}
したがって、中心値 とスケーリングパラメータを推定するより堅牢な手段 が必要である。 1 つの簡単な方法は、サンプルの中央値を の推定値として取り 、サンプルの四分 位範囲 の半分をの推定値として取ることである 。 他にも、より正確で堅牢な方法が開発されている [20] [21]。 たとえば、サンプル 順序統計量 の中央 24% の 切り捨て平均は、 サンプルの中央値またはサンプル全体の平均を使用するよりも効率的な の 推定値を生成する。 [22] [23]ただし、コーシー分布の 太い裾 のため 、サンプルの 24% 以上を使用すると推定値の効率は低下する。 [22] [23]
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma }
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma }
x
0
{\displaystyle x_{0}}
最大尤度法は パラメータとを 推定するためにも使用できます 。しかし、これは高次多項式の根を見つける必要があり、局所的最大値を表す根が複数存在する可能性があるため、複雑になる傾向があります。 [24] また、最大尤度推定量は漸近的には効率的ですが、小さなサンプルでは比較的非効率的です。 [25] [26] サンプルサイズに対するコーシー分布の対数尤度関数は次 のとおりです。
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma }
n
{\displaystyle n}
ℓ
^
(
x
1
,
…
,
x
n
∣
x
0
,
γ
)
=
−
n
log
(
γ
π
)
−
∑
i
=
1
n
log
(
1
+
(
x
i
−
x
0
γ
)
2
)
{\displaystyle {\hat {\ell }}(x_{1},\dotsc ,x_{n}\mid \!x_{0},\gamma )=-n\log(\gamma \pi )-\sum _{i=1}^{n}\log \left(1+\left({\frac {x_{i}-x_{0}}{\gamma }}\right)^{2}\right)}
および に関して対数尤度関数を最大化するために一次導関数をとると 、 次の連立方程式が生成されます。
x
0
{\displaystyle x_{0}}
γ
{\displaystyle \gamma }
d
ℓ
d
x
0
=
∑
i
=
1
n
2
(
x
i
−
x
0
)
γ
2
+
(
x
i
−
x
0
)
2
=
0
{\displaystyle {\frac {d\ell }{dx_{0}}}=\sum _{i=1}^{n}{\frac {2(x_{i}-x_{0})}{\gamma ^{2}+\left(x_{i}-\!x_{0}\right)^{2}}}=0}
d
ℓ
d
γ
=
∑
i
=
1
n
2
(
x
i
−
x
0
)
2
γ
(
γ
2
+
(
x
i
−
x
0
)
2
)
−
n
γ
=
0
{\displaystyle {\frac {d\ell }{d\gamma }}=\sum _{i=1}^{n}{\frac {2\left(x_{i}-x_{0}\right)^{2}}{\gamma (\gamma ^{2}+\left(x_{i}-x_{0}\right)^{2})}}-{\frac {n}{\gamma }}=0}
ご了承ください
∑
i
=
1
n
(
x
i
−
x
0
)
2
γ
2
+
(
x
i
−
x
0
)
2
{\displaystyle \sum _{i=1}^{n}{\frac {\left(x_{i}-x_{0}\right)^{2}}{\gamma ^{2}+\left(x_{i}-x_{0}\right)^{2}}}}
は単調関数であり 、解 は
γ
{\displaystyle \gamma }
γ
{\displaystyle \gamma }
min
|
x
i
−
x
0
|
≤
γ
≤
max
|
x
i
−
x
0
|
.
{\displaystyle \min |x_{i}-x_{0}|\leq \gamma \leq \max |x_{i}-x_{0}|.}
についてのみ解くには 、 次数 の多項式を解く必要があります 。 [24] また、についてのみ解くには、 次数 の多項式を解く必要があります。したがって、1 つのパラメータについて解く場合でも、両方のパラメータを同時に解く場合でも、通常はコンピュータ上で 数値 解を求める必要があります。最尤推定の利点は漸近効率です。 標本中央値を使用して推定した場合の漸近効率は、 最尤推定値と比較して約 81% にすぎません。 [23] [27] 中間 24% の順序統計量を使用した切り捨て標本平均は、 最大尤度推定値と比較して約 88% の漸近効率の の推定値です。 [23] ニュートン法 を使用して最大尤度推定値の解を求める場合 、中間 24% の順序統計量を の初期解として使用できます 。
x
0
{\displaystyle x_{0}}
2
n
−
1
{\displaystyle 2n-1}
γ
{\displaystyle \,\!\gamma }
2
n
{\displaystyle 2n}
x
0
{\displaystyle x_{0}}
x
0
{\displaystyle x_{0}}
x
0
{\displaystyle x_{0}}
x
0
{\displaystyle x_{0}}
位置 0 の Cauchy 変数の場合、 形状パラメータは
絶対値の中央値を使用して推定できます。
X
∼
C
a
u
c
h
y
(
0
,
γ
)
{\displaystyle X\sim \mathrm {Cauchy} (0,\gamma )}
median
(
|
X
|
)
=
γ
{\displaystyle \operatorname {median} (|X|)=\gamma }
一般的な
Cauchy
(
0
,
1
)
∼
t
(
d
f
=
1
)
{\displaystyle \operatorname {Cauchy} (0,1)\sim {\textrm {t}}(\mathrm {df} =1)\,}
スチューデントの t 分布
Cauchy
(
μ
,
σ
)
∼
t
(
d
f
=
1
)
(
μ
,
σ
)
{\displaystyle \operatorname {Cauchy} (\mu ,\sigma )\sim {\textrm {t}}_{(\mathrm {df} =1)}(\mu ,\sigma )\,}
非標準化スチューデントの t 分布
独立している 場合は
X
,
Y
∼
N
(
0
,
1
)
X
,
Y
{\displaystyle X,Y\sim {\textrm {N}}(0,1)\,X,Y}
X
Y
∼
Cauchy
(
0
,
1
)
{\displaystyle {\tfrac {X}{Y}}\sim {\textrm {Cauchy}}(0,1)\,}
もし そうなら
X
∼
U
(
0
,
1
)
{\displaystyle X\sim {\textrm {U}}(0,1)\,}
tan
(
π
(
X
−
1
2
)
)
∼
Cauchy
(
0
,
1
)
{\displaystyle \tan \left(\pi \left(X-{\tfrac {1}{2}}\right)\right)\sim {\textrm {Cauchy}}(0,1)\,}
もし そうなら
X
∼
L
o
g
-
C
a
u
c
h
y
(
0
,
1
)
{\displaystyle X\sim \operatorname {Log-Cauchy} (0,1)}
ln
(
X
)
∼
Cauchy
(
0
,
1
)
{\displaystyle \ln(X)\sim {\textrm {Cauchy}}(0,1)}
もし そうなら
X
∼
Cauchy
(
x
0
,
γ
)
{\displaystyle X\sim \operatorname {Cauchy} (x_{0},\gamma )}
1
X
∼
Cauchy
(
x
0
x
0
2
+
γ
2
,
γ
x
0
2
+
γ
2
)
{\displaystyle {\tfrac {1}{X}}\sim \operatorname {Cauchy} \left({\tfrac {x_{0}}{x_{0}^{2}+\gamma ^{2}}},{\tfrac {\gamma }{x_{0}^{2}+\gamma ^{2}}}\right)}
コーシー分布は タイプ4の ピアソン分布の極限ケースである [ 要出典 ]
コーシー分布は タイプ7の ピアソン分布の特殊なケースである 。[1]
コーシー分布は 安定分布 です。つまり、 の場合 、 となります 。
X
∼
Stable
(
1
,
0
,
γ
,
μ
)
{\displaystyle X\sim {\textrm {Stable}}(1,0,\gamma ,\mu )}
X
∼
Cauchy
(
μ
,
γ
)
{\displaystyle X\sim \operatorname {Cauchy} (\mu ,\gamma )}
コーシー分布は双曲分布 の特異極限である [ 要出典 ]
円上の値を取るラップされたコーシー分布は 、 コーシー分布を円の周りでラップすることによって派生します。
、 の場合 、 となります 。半コーシー分布の場合、 と設定することで関係式が成立します 。
X
∼
N
(
0
,
1
)
{\displaystyle X\sim {\textrm {N}}(0,1)}
Z
∼
I
n
v
e
r
s
e
-
G
a
m
m
a
(
1
/
2
,
s
2
/
2
)
{\displaystyle Z\sim \operatorname {Inverse-Gamma} (1/2,s^{2}/2)}
Y
=
μ
+
X
Z
∼
Cauchy
(
μ
,
s
)
{\displaystyle Y=\mu +X{\sqrt {Z}}\sim \operatorname {Cauchy} (\mu ,s)}
X
∼
N
(
0
,
1
)
I
{
X
≥
0
}
{\displaystyle X\sim {\textrm {N}}(0,1)I\{X\geq 0\}}
レヴィ尺度
コーシー分布は 指数1の 安定分布 です。このようなパラメータの安定分布の レヴィ・ヒンチン表現は 次のように与えられます 。
γ
{\displaystyle \gamma }
X
∼
Stable
(
γ
,
0
,
0
)
{\displaystyle X\sim \operatorname {Stable} (\gamma ,0,0)\,}
E
(
e
i
x
X
)
=
exp
(
∫
R
(
e
i
x
y
−
1
)
Π
γ
(
d
y
)
)
{\displaystyle \operatorname {E} \left(e^{ixX}\right)=\exp \left(\int _{\mathbb {R} }(e^{ixy}-1)\Pi _{\gamma }(dy)\right)}
どこ
Π
γ
(
d
y
)
=
(
c
1
,
γ
1
y
1
+
γ
1
{
y
>
0
}
+
c
2
,
γ
1
|
y
|
1
+
γ
1
{
y
<
0
}
)
d
y
{\displaystyle \Pi _{\gamma }(dy)=\left(c_{1,\gamma }{\frac {1}{y^{1+\gamma }}}1_{\left\{y>0\right\}}+c_{2,\gamma }{\frac {1}{|y|^{1+\gamma }}}1_{\left\{y<0\right\}}\right)\,dy}
と 明示的に表現することができる。 [28] コーシー分布の 場合、となる 。
c
1
,
γ
,
c
2
,
γ
{\displaystyle c_{1,\gamma },c_{2,\gamma }}
γ
=
1
{\displaystyle \gamma =1}
c
1
,
γ
=
c
2
,
γ
{\displaystyle c_{1,\gamma }=c_{2,\gamma }}
この最後の表現は、次の式の結果である。
π
|
x
|
=
PV
∫
R
∖
{
0
}
(
1
−
e
i
x
y
)
d
y
y
2
{\displaystyle \pi |x|=\operatorname {PV} \int _{\mathbb {R} \smallsetminus \lbrace 0\rbrace }(1-e^{ixy})\,{\frac {dy}{y^{2}}}}
多変量コーシー分布
ランダム ベクトルは 、その成分のすべての線形結合がコーシー分布に従う場合、多変量コーシー分布に従うと言われます 。つまり、任意の定数ベクトルに対して 、ランダム変数は 単変量コーシー分布に従う必要があります。 [29] 多変量コーシー分布の特性関数は次のように与えられます。
X
=
(
X
1
,
…
,
X
k
)
T
{\displaystyle X=(X_{1},\ldots ,X_{k})^{T}}
Y
=
a
1
X
1
+
⋯
+
a
k
X
k
{\displaystyle Y=a_{1}X_{1}+\cdots +a_{k}X_{k}}
a
∈
R
k
{\displaystyle a\in \mathbb {R} ^{k}}
Y
=
a
T
X
{\displaystyle Y=a^{T}X}
φ
X
(
t
)
=
e
i
x
0
(
t
)
−
γ
(
t
)
,
{\displaystyle \varphi _{X}(t)=e^{ix_{0}(t)-\gamma (t)},\!}
ここで 、およびは 1次の 同次関数 と 1次の正同次関数 を持つ実関数である。 [29] より正式には: [29]
x
0
(
t
)
{\displaystyle x_{0}(t)}
γ
(
t
)
{\displaystyle \gamma (t)}
x
0
(
t
)
{\displaystyle x_{0}(t)}
γ
(
t
)
{\displaystyle \gamma (t)}
x
0
(
a
t
)
=
a
x
0
(
t
)
,
{\displaystyle x_{0}(at)=ax_{0}(t),}
γ
(
a
t
)
=
|
a
|
γ
(
t
)
,
{\displaystyle \gamma (at)=|a|\gamma (t),}
すべてに対して 。
t
{\displaystyle t}
二変量コーシー分布の例は次のようになる。 [30]
f
(
x
,
y
;
x
0
,
y
0
,
γ
)
=
1
2
π
[
γ
(
(
x
−
x
0
)
2
+
(
y
−
y
0
)
2
+
γ
2
)
3
/
2
]
.
{\displaystyle f(x,y;x_{0},y_{0},\gamma )={1 \over 2\pi }\left[{\gamma \over ((x-x_{0})^{2}+(y-y_{0})^{2}+\gamma ^{2})^{3/2}}\right].}
この例では、と間の共分散 が0である にもかかわらず 、 とが 統計的に独立してい ないことに注意してください 。 [30]
x
{\displaystyle x}
y
{\displaystyle y}
x
{\displaystyle x}
y
{\displaystyle y}
この式は複素変数に対しても書くことができます。複素コーシーの確率密度関数は次のようになります。
f
(
z
;
z
0
,
γ
)
=
1
2
π
[
γ
(
|
z
−
z
0
|
2
+
γ
2
)
3
/
2
]
.
{\displaystyle f(z;z_{0},\gamma )={1 \over 2\pi }\left[{\gamma \over (|z-z_{0}|^{2}+\gamma ^{2})^{3/2}}\right].}
標準コーシー分布が自由度 1 のスチューデント t 分布であるのと同様に、多次元コーシー密度は自由度 1 の 多変量スチューデント分布 です。自由度 1 の次元スチューデント分布
の密度は次のとおりです。
k
{\displaystyle k}
f
(
x
;
μ
,
Σ
,
k
)
=
Γ
(
1
+
k
2
)
Γ
(
1
2
)
π
k
2
|
Σ
|
1
2
[
1
+
(
x
−
μ
)
T
Σ
−
1
(
x
−
μ
)
]
1
+
k
2
.
{\displaystyle f({\mathbf {x} };{\mathbf {\mu } },{\mathbf {\Sigma } },k)={\frac {\Gamma \left({\frac {1+k}{2}}\right)}{\Gamma ({\frac {1}{2}})\pi ^{\frac {k}{2}}\left|{\mathbf {\Sigma } }\right|^{\frac {1}{2}}\left[1+({\mathbf {x} }-{\mathbf {\mu } })^{T}{\mathbf {\Sigma } }^{-1}({\mathbf {x} }-{\mathbf {\mu } })\right]^{\frac {1+k}{2}}}}.}
多次元コーシー分布の特性は、多変量スチューデント分布の特殊なケースになります。
発生と応用
一般的に
分光学 において、コーシー分布は、 スペクトル線 の形状を表します。このスペクトル 線は、すべての原子がスペクトル線の形状に含まれる周波数範囲で同じように相互作用し、均一な広がりが生じます 。 均一な広がりは多くのメカニズムによって発生しますが、最も顕著なのは 衝突による広がりです 。 [31] 寿命や自然広がり によっても、コーシー分布で表されるスペクトル線の形状が生じます。
コーシー分布またはその変換の応用は、 指数関数的増加を 扱う分野で見出すことができます。1958年のホワイトの論文 [32] は、方程式の 推定値の検定統計量を導き出し 、通常の最小二乗法を使用して最大尤度推定値を求めると、統計量の標本分布がコーシー分布になることを示しました。
β
^
{\displaystyle {\hat {\beta }}}
x
t
+
1
=
β
x
t
+
ε
t
+
1
,
β
>
1
{\displaystyle x_{t+1}=\beta {x}_{t}+\varepsilon _{t+1},\beta >1}
CumFreq を用いて最大1日降雨量に累積コーシー分布を当てはめた 。 分布の当てはめも参照 [33]
コーシー分布は、回転する物体の観測値の分布であることが多い。これに関する古典的な例は、カモメの灯台問題 [34] と呼ばれ、上記のセクションでは素粒子物理学におけるブライト・ウィグナー分布と呼ばれている。
水文学 では 、コーシー分布は年間最大日雨量や河川流量などの極端な事象に適用されます。青い図は、順位付けされた月間最大日雨量にコーシー分布を当てはめた例を示しており、 二項分布 に基づく90% 信頼区間も表示しています。降雨量データは、 累積頻度分析 の一部として 位置をプロットする ことによって表されます 。
ローレンツモデルによれば、複素 誘電率 の虚数部の表現はコーシー分布です。
計算ファイナンス における ファットテール をモデル化するための追加の分布として、コーシー分布は ガウス分布 よりもはるかに大きな極端なリスクの確率を生み出すVAR( リスク値 )をモデル化するために使用できます 。 [35]
相対論的ブライト・ウィグナー分布
原子核 物理学および 素粒子物理学 では、 共鳴 のエネルギープロファイルは 相対論的なブライト・ウィグナー分布 によって記述さ れ、コーシー分布は(非相対論的な)ブライト・ウィグナー分布である。 [ 要出典 ]
歴史
コーシー分布 (下) のサンプルから平均と標準偏差を推定すると、 正規分布 (上) の場合のようにサンプルのサイズが大きくなっても収束しません。下のグラフに見られるように、推定値には任意の大きなジャンプが生じる可能性があります。(クリックして拡大)
コーシー分布の密度関数の形を持つ関数は、 1659年に フェルマー によって幾何学的に研究され、後に マリア・ガエターナ・アニェージが 1748年に微積分の教科書に例として含めたことから、 アニェージの魔女 として知られるようになった。その名前にもかかわらず、コーシー分布の特性の最初の明示的な分析は1824年にフランスの数学者 ポアソン によって発表され、コーシーがそれに関連付けられたのは1853年の学術的論争のときだけであった。 [36] ポアソンは、そのような分布に従う観測の平均をとると、 標準偏差は どの有限の数にも収束しないことに注目した。したがって、 ラプラスがそのような分布で 中心極限定理 を使用したのは、有限の平均と分散を仮定していたため不適切であった。それにもかかわらず、ポアソンはこの問題を重要視せず、 ビエネメ とは対照的であった。ビエネメ はこの問題をめぐってコーシーと長い論争を繰り広げた。
参照
参考文献
^ abc NL Johnson; S. Kotz; N. Balakrishnan (1994). 連続一変量分布、第1巻 。ニューヨーク:Wiley。 、第16章。
^ フェラー、ウィリアム(1971)。 確率論とその応用入門、第2巻 (第2版)。ニューヨーク:ジョン・ワイリー・アンド・サンズ社。pp. 704。ISBN 978-0-471-25709-7 。
^ 「ローレンツ関数」 。MathWorld。Wolfram Research 。 2024年 10月27日 閲覧 。
^ ライリー、ケン F.; ホブソン、マイケル P.; ベンス、スティーブン J. (2006). 物理学と工学のための数学的手法 (第 3 版). ケンブリッジ、イギリス: ケンブリッジ大学出版局. pp. 1333. ISBN 978-0-511-16842-0 。
^ Balakrishnan, N.; Nevrozov, VB (2003). 統計分布入門(第1版). ホーボーケン、ニュージャージー:John Wiley & Sons Inc. pp. 305. ISBN 0-471-42798-5 。
^ Li, Rui; Nadarajah, Saralees (2020-03-01). 「スチューデントのt分布とその一般化のレビュー」. 経験的経済学 . 58 (3): 1461–1490. doi :10.1007/s00181-018-1570-0. ISSN 1435-8921.
^ Pillai N.; Meng, XL (2016). 「コーシーとレヴィとの予期せぬ出会い」. 統計年報 . 44 (5): 2089–2097. arXiv : 1505.01957 . doi :10.1214/15-AOS1407. S2CID 31582370.
^ キャンベル B. 読む; N.バラクリシュナン。ブラニ・ヴィダコビッチ。サミュエル・コッツ (2006)。 統計科学百科事典 (第 2 版)。 ジョン・ワイリー&サンズ 。 p. 778.ISBN 978-0-471-15044-2 。
^ナイト、フランク B. ( 1976 ) 。 「コーシー型の特徴付け」。 アメリカ数学会紀要 。55 (1): 130–135。doi : 10.2307/2041858。JSTOR 2041858 。
^ 「コーシー中心極限のアップデート」。 量子計算 。2022年11月13日。 2023年 6月21日 閲覧 。
^ Frederic, Chyzak; Nielsen, Frank (2019). 「コーシー分布間のカルバック・ライブラー距離の閉形式の公式」 arXiv : 1905.10965 [cs.IT].
^ Nielsen, Frank; Okamura, Kazuki (2023). 「コーシー分布間のfダイバージェンスについて」. IEEE Transactions on Information Theory . 69 (5): 3150–3171. arXiv : 2101.12459 . doi :10.1109/TIT.2022.3231645. S2CID 231728407.
^ Vasicek, Oldrich (1976). 「サンプルエントロピーに基づく正規性検定」. 英国王立統計学会誌、シリーズ B. 38 ( 1): 54–59. doi :10.1111/j.2517-6161.1976.tb01566.x.
^ Park, Sung Y.; Bera, Anil K. (2009). 「最大エントロピー自己回帰条件付き異分散モデル」 (PDF) . Journal of Econometrics . 150 (2). Elsevier: 219–230. doi :10.1016/j.jeconom.2008.12.014. 2011-09-30 に オリジナル (PDF)からアーカイブ 。2011-06-02 に取得 。
^ ab Kyle Siegrist. 「Cauchy Distribution」. ランダム . 2021年7月9日時点のオリジナルよりアーカイブ 。 2021年 7月5日 閲覧。
^ Hampel, Frank (1998)、「統計は難しすぎるか?」 (PDF) 、 Canadian Journal of Statistics 、 26 (3): 497–513、 doi :10.2307/3315772、 hdl : 20.500.11850/145503 、 JSTOR 3315772、 S2CID 53117661、2022-01-25にオリジナルからアーカイブ、 2019-09-25に 取得 。
^ レモンズ、ドン・S. (2002)、「物理学における確率過程入門」、 アメリカ物理学ジャーナル 、 71 (2)、ジョンズ・ホプキンス大学出版局: 35、 Bibcode :2003AmJPh..71..191L、 doi :10.1119/1.1526134、 ISBN 0-8018-6866-1
^ ab McCullagh, P. 、「Conditional inference and Cauchy models」、 Biometrika 、第79巻(1992年)、247~259ページ。PDFは2010年6月10日にWayback Machine でMcCullaghのホームページからアーカイブされています 。
^ 「標本平均の不安定性の図解」。2017年3月24日時点のオリジナルよりアーカイブ 。 2014年11月22日 閲覧。
^ Cane, Gwenda J. (1974). 「 サンプル分位数に基づくコーシー分布のパラメータの線形推定」。 アメリカ統計学会誌 。69 (345): 243–245. doi :10.1080/01621459.1974.10480163. JSTOR 2285535.
^ Zhang, Jin (2010). 「コーシー分布の位置パラメータの非常に効率的なL推定量」. 計算統計 . 25 (1): 97–105. doi :10.1007/s00180-009-0163-y. S2CID 123586208.
^ ab Rothenberg, Thomas J.; Fisher, Franklin, M.; Tilanus, CB (1964). 「コーシー標本からの推定に関する注記」 アメリカ統計学会誌 . 59 (306): 460–463. doi :10.1080/01621459.1964.10482170. {{cite journal}}: CS1 maint: multiple names: authors list (link)
^ abcd Bloch, Daniel (1966). 「コーシー分布の位置パラメータの推定に関する注記」 アメリカ統計学会誌 . 61 (316): 852–855. doi :10.1080/01621459.1966.10480912. JSTOR 2282794.
^ ab Ferguson, Thomas S. (1978). 「サンプルサイズ 3 および 4 の Cauchy 分布のパラメータの最大尤度推定値」。 アメリカ 統計学会誌 。73 (361): 211–213。doi : 10.1080/01621459.1978.10480031。JSTOR 2286549。
^ Cohen Freue, Gabriella V. (2007). 「コーシー位置パラメータのピットマン推定量」 (PDF) . Journal of Statistical Planning and Inference . 137 (6): 1901. doi :10.1016/j.jspi.2006.05.002. 2011-08-16にオリジナル (PDF) からアーカイブ。
^ Wilcox, Rand (2012). ロバスト推定と仮説検定入門 . Elsevier.
^ Barnett, VD (1966). 「コーシー分布の位置の順序統計推定量」 アメリカ統計学会誌 . 61 (316): 1205–1218. doi :10.1080/01621459.1966.10482205. JSTOR 2283210.
^ Kyprianou, Andreas (2009). Lévy processes and continuously-state branching processes:part I (PDF) . p. 11. 2016-03-03 のオリジナルからアーカイブ (PDF) 。2016-05-04 に取得 。
^ abc Ferguson, Thomas S. (1962). 「対称二変量コーシー分布の表現」. 数理統計年報 . 33 (4): 1256–1266. doi : 10.1214/aoms/1177704357 . JSTOR 2237984. 2017年1月7日 閲覧 。
^ ab Molenberghs, Geert; Lesaffre, Emmanuel (1997). 「与えられた周辺値と依存関数で二変量密度を近似する非線形積分方程式」 (PDF) 。 Statistica Sinica 。 7 : 713–738。 2009-09-14 にオリジナル (PDF) からアーカイブ。
^ E. Hecht (1987). Optics (第2版). Addison-Wesley . p. 603.
^ White, JS (1958年12月). 「爆発的ケースにおける連続相関係数の極限分布」. 数理統計年報 . 29 (4): 1188–1197. doi : 10.1214/aoms/1177706450 .
^ 「CumFreq、累積頻度分析と確率分布フィッティングのための無料ソフトウェア」。2018年2月21日時点のオリジナルよりアーカイブ。
^ Gull, SF (1988) ベイジアン帰納的推論と最大エントロピー。Kluwer Academic Publishers、ベルリン。https://doi.org/10.1007/978-94-009-3049-0_4 2022-01-25 に Wayback Machineでアーカイブ
^ Tong Liu (2012)、ガウス分布とコーシー分布の中間分布。https://arxiv.org/pdf/1208.5109.pdf 2020-06-24に Wayback Machineでアーカイブ
^ コーシーとアグネージの魔女、 統計学の表 、S.M.スティグラー、ハーバード大学、1999年、第18章
外部リンク