確率論では、重裾分布とは、裾が指数的に制限されていない確率分布のことである。 [1]つまり、指数分布よりも裾が重い分布である。多くの応用では、分布の右裾が重要になるが、分布の左裾が重い場合や、両方の裾が重い場合もある。
重裾分布には、太裾分布、長裾分布、準指数分布という3つの重要なサブクラスがあります。実際には、一般的に使用される重裾分布はすべて、 Jozef Teugelsによって導入された準指数クラスに属します。[2]
ヘビーテールという用語の使用については、まだ意見の相違があります。使用されている定義は他に 2 つあります。一部の著者は、この用語を、すべてのべき乗モーメントが有限ではない分布を指すために使用し、他の著者は、有限の分散を持たない分布を指すために使用します。この記事で示す定義は、最も一般的に使用されている定義であり、代替定義に含まれるすべての分布と、すべてのべき乗モーメントを持ちながらも一般にヘビーテールであると見なされる対数正規分布などの分布を含みます。(ヘビーテールは、正規分布よりも重い裾を持つ分布に対して使用されることがあります。)
定義
ヘビーテール分布の定義
分布関数Fを持つランダム変数 Xの分布は、Xのモーメント生成関数MX ( t )がすべてのt>0に対して無限大である場合、 重い(右の)裾を持つと言われる。[3]
つまり
- [4]
これは裾分布関数で表すこともできる。
として
ロングテール分布の定義
分布関数Fを持つ確率変数 Xの分布は、すべてのt > 0 に対して、
または同等
これは、右裾が長く裾が広がる分布量に対する直感的な解釈であり、長く裾が広がる量がある高い水準を超えると、他のより高い水準を超える確率が 1 に近づくというものです。
すべてのロングテール分布はヘビーテールですが、その逆は偽であり、ロングテールではないヘビーテール分布を構築することは可能です。
準指数分布
指数関数未満性は、確率分布の畳み込みによって定義されます。共通の分布関数を持つ2 つの独立した同一分布のランダム変数 について、とそれ自身との畳み込み(畳み込み平方と表記) は、ルベーグ・スティルチェス積分を使用して次のように定義されます。
n倍畳み込みは次の規則によって帰納的に定義されます。
裾分布関数は次のように定義されます。
正の半直線上の分布は、指数関数的分布である[1] [5] [2]。
これは[6]、任意のに対して、
これの確率論的解釈[6]は、共通の分布を持つ独立した確率変数の和に対して、
これは、単一の大きなジャンプの原理[7]またはカタストロフィー原理[8]としてよく知られています。
実数直線上の 分布が指数関数的であるのは、分布が指数関数的である場合である。[9]ここで、正の半直線の指示関数である 。あるいは、実数直線上でサポートされるランダム変数が指数関数的である場合、かつ、が指数関数的である場合に限ります。
すべての準指数分布は長い裾を持ちますが、準指数分布ではない長い裾を持つ分布の例を作成することもできます。
一般的な重尾分布
一般的に使用される重裾分布はすべて準指数分布である。[6]
片側検定には次のものがあります:
- パレート分布;
- 対数正規分布;
- レヴィ分布;
- 形状パラメータが 0 より大きく 1 より小さいワイブル分布。
- バー分布;
- 対数ロジスティック分布;
- 対数ガンマ分布;
- フレシェ分布;
- q-ガウス分布;
- 対数コーシー分布は対数減衰を示し、パレート分布よりも重い裾を持つため、「超重い裾」を持つと言われることもある。 [10] [11]
両側型のものには次のものがあります:
- コーシー分布は、安定分布と t 分布の両方の特殊なケースです。
- 安定分布の族 [12] 。ただし、その族内の正規分布の特殊なケースは除く。安定分布の中には片側(または半直線でサポートされている)のものもあります。例えば、レヴィ分布を参照してください。また、ロングテール分布を持つ金融モデルやボラティリティクラスタリングも参照してください。
- t分布。
- 歪んだ対数正規カスケード分布。[13]
太い裾を持つ分布との関係
裾が厚い分布とは、大きな x に対して確率密度関数がべき乗としてゼロになる分布です。このようなべき乗は常に指数分布の確率密度関数によって下方に制限されるため、裾が厚い分布は常に裾が重い分布です。ただし、分布によっては、ゼロになるのが指数関数よりも遅い (つまり裾が重い) ものの、べき乗よりも速い (つまり裾が厚くない) ものもあります。例としては、対数正規分布[矛盾]があります。ただし、対数ロジスティック分布やパレート分布など、他の多くの裾が重い分布も裾が厚い分布です。
末尾インデックスの推定
テールインデックス推定の問題にはパラメトリック[6]とノンパラメトリック[14]のアプローチがある。 [どのように定義されるか? ]
パラメトリックアプローチを使用してテールインデックスを推定するために、一部の著者は GEV 分布またはパレート分布を使用します。また、最大尤度推定量 (MLE) を適用する場合もあります。
Pickand の末尾インデックス推定量
独立かつ同一の密度関数のランダムシーケンスを用いると、 一般化極値密度の最大吸引領域[15]は、 となる。 かつ の場合には、 ピックアンドのテールインデックス推定は[6] [15]となる。
ここで、この推定量は確率的に に収束します。
ヒルのテールインデックス推定量
を分布関数(一般化極値分布の最大吸引領域)を持つ独立かつ同一に分布するランダム変数のシーケンスとします。ここで、サンプルパスはであり、はサンプルサイズです。 が中間次数シーケンス、つまり、 の場合、ヒルテールインデックス推定量は[16]です。
ここで、 はの- 番目の順序統計量です。この推定量は確率的に に収束し、 が高次の正規変動特性[17]に基づいて制限されている場合、漸近的に正規です 。[18]一貫性と漸近正規性は、 が観測されるか、または誤って指定されたモデルや従属する誤差を持つモデルを含む、モデルと推定量の大規模なクラスから計算された残差またはフィルタリングされたデータであるかに関係なく、従属および異質なシーケンスの大規模なクラスに拡張されます。[21] [22] [23] Pickand と Hill のテールインデックス推定量はどちらも、順序統計量の対数を一般的に使用することに注意してください。[24]
テールインデックスの比率推定値
テールインデックスの比率推定量(RE推定量)は、ゴールドとスミスによって導入されました。[25] これはヒルの推定量と同様に構築されていますが、非ランダムな「調整パラメータ」を使用します。
ヒル型推定量とRE型推定量の比較についてはNovakの論文を参照のこと。[14]
ソフトウェア
- aest Archived 2020-11-25 at the Wayback Machine、ヘビーテール指数を推定するCツール。 [26]
重尾密度の推定
ヘビーテールおよびスーパーヘビーテールの確率密度関数を推定するノンパラメトリックアプローチは、マルコビッチによって提案されました。[27]これらは、可変帯域幅とロングテールカーネル推定量に基づくアプローチ、有限または無限間隔で予備データを新しいランダム変数に変換するアプローチ(推定と得られた密度推定値の逆変換に便利です)、および密度のテールに特定のパラメトリックモデルを提供し、密度のモードを近似するノンパラメトリックモデルを提供する「つなぎ合わせアプローチ」です。ノンパラメトリック推定量では、カーネル推定量の帯域幅やヒストグラムのビン幅などの調整(平滑化)パラメータを適切に選択する必要があります。このような選択を行うデータ駆動型のよく知られた方法は、クロスバリデーションとその修正、平均二乗誤差(MSE)の最小化とその漸近線、およびそれらの上限に基づく方法です。[28]コルモゴロフ・スミルノフ、フォン・ミーゼス、アンダーソン・ダーリングなどのよく知られたノンパラメトリック統計を分布関数(dfs)の空間における測定基準として、そして後者の統計量の分位数を既知の不確実性または不一致値として使用する不一致法は、 [27] で見つけることができます。ブートストラップは、異なる再サンプル選択スキームによる未知のMSEの近似値を使用して平滑化パラメータを見つけるための別のツールです。たとえば、[29]を参照してください。
参照
参考文献
- ^ abc Asmussen, SR (2003). 「GI/G/1 の定常特性」.応用確率とキュー. 確率モデルと応用確率. 第 51 巻. pp. 266–301. doi :10.1007/0-387-21525-5_10. ISBN 978-0-387-00211-8。
- ^ ab Teugels, Jozef L. (1975). 「The Class of Subexponential Distributions」Annals of Probability . 3 (6). University of Louvain . doi : 10.1214/aop/1176996225 . 2019年4月7日閲覧。
- ^ Rolski、Schmidli、Scmidt、Teugels、保険と金融のための確率過程、1999年
- ^ S. Foss、D. Korshunov、S. Zachary、「Heavy-Tailed and Subexponential Distributions 入門」、Springer Science & Business Media、2013 年 5 月 21 日
- ^ Chistyakov, VP (1964). 「独立した正のランダム変数の合計に関する定理と分岐ランダムプロセスへのその応用」。ResearchGate 。 2019年4月7日閲覧。
- ^ abcde Embrechts P.; Klueppelberg C.; Mikosch T. (1997).保険と金融のための極限イベントのモデリング. 確率モデルと応用確率. 第33巻. ベルリン: Springer. doi :10.1007/978-3-642-33483-2. ISBN 978-3-642-08242-9。
- ^ Foss, S.; Konstantopoulos, T.; Zachary, S. (2007). 「Heavy-Tailed Increments による離散および連続時間変調ランダムウォーク」(PDF) . Journal of Theoretical Probability . 20 (3): 581. arXiv : math/0509605 . CiteSeerX 10.1.1.210.1699 . doi :10.1007/s10959-007-0081-2. S2CID 3047753.
- ^ Wierman, Adam (2014 年 1 月 9 日)。「大惨事、陰謀、および準指数分布 (パート III)」。Rigor + Relevance ブログ。RSRG、Caltech。2014年1 月 9 日閲覧。
- ^ Willekens、E. (1986)。 「実数直線上の準指数関数性」。技術レポート。 KU ルーヴェン。
- ^ Falk, M., Hüsler, J. & Reiss, R. (2010).少数の法則: 極値と稀な事象. Springer. p. 80. ISBN 978-3-0348-0008-2。
{{cite book}}: CS1 maint: multiple names: authors list (link) - ^ Alves, MIF, de Haan, L. & Neves, C. (2006 年 3 月 10 日). 「重い裾と超重い裾の分布の統計的推論」(PDF) 。2007年 6 月 23 日時点のオリジナル(PDF)からアーカイブ。2011 年11 月 1 日に閲覧。
{{cite web}}: CS1 maint: multiple names: authors list (link) - ^ John P. Nolan (2009). 「安定した分布: 裾の厚いデータのためのモデル」(PDF)。2011-07-17 にオリジナル(PDF)からアーカイブ。2009-02-21に閲覧。
- ^ Stephen Lihn (2009). 「Skew Lognormal Cascade Distribution」。2014年4月7日時点のオリジナルよりアーカイブ。2009年6月12日閲覧。
- ^ ab Novak SY (2011).極値法と金融への応用. ロンドン: CRC. ISBN 978-1-43983-574-6。
- ^ ab Pickands III, James (1975年1月). 「極度順序統計を用いた統計的推論」. 統計年報. 3 (1): 119–131. doi : 10.1214/aos/1176343003 . JSTOR 2958083.
- ^ Hill BM (1975) 分布の裾についての推論に対する単純で一般的なアプローチ。Ann. Stat., v. 3, 1163–1174。
- ^ Hall, P.(1982) 正規分布の指数のいくつかの推定について JR Stat. Soc. Ser. B., v. 44, 37–42.
- ^ Haeusler, E. および JL Teugels (1985) 正規分布の指数に対する Hill 推定量の漸近正規性について。Ann. Stat., v. 13, 743–756。
- ^ Hsing, T. (1991) 従属データを使用したテールインデックス推定について。Ann. Stat., v. 19, 1547–1569。
- ^ Hill, J. (2010) 従属的かつ異質なデータのテールインデックス推定について Econometric Th., v. 26, 1398–1436.
- ^ Resnick, S. および Starica, C. (1997)。自己回帰データに対する Hill 推定量の漸近挙動。Comm. Statist. Stochastic Models 13、703–721。
- ^ Ling, S. および Peng, L. (2004). ARMA モデルのテール インデックスに対する Hill 推定量。J. Statist. Plann. Inference 123, 279–293。
- ^ Hill, JB (2015). フィルタリングされた従属時系列のテールインデックス推定。Stat. Sin. 25, 609–630。
- ^ Lee, Seyoon; Kim, Joseph HT (2019). 「指数化一般化パレート分布: 極値理論に向けた特性と応用」. Communications in Statistics - Theory and Methods . 48 (8): 2014–2038. arXiv : 1708.01686 . doi :10.1080/03610926.2018.1441418. S2CID 88514574.
- ^ Goldie CM、Smith RL (1987) 剰余を伴う緩やかな変化:理論と応用。Quart. J. Math. Oxford、v. 38、45–71。
- ^ Crovella, ME; Taqqu, MS (1999). 「スケーリング特性からヘビーテール指数を推定する」. Methodology and Computing in Applied Probability . 1 : 55–79. doi :10.1023/A:1010012224103. S2CID 8917289. 2007-02-06 にオリジナルからアーカイブ。 2015-09-03に取得。
- ^ ab Markovich NM (2007).単変量ヘビーテールデータのノンパラメトリック分析: 研究と実践. Chitester: Wiley. ISBN 978-0-470-72359-3。
- ^ Wand MP、Jones MC (1995)。カーネルスムージング。ニューヨーク:チャップマンアンドホール。ISBN 978-0412552700。
- ^ Hall P. (1992).ブートストラップとエッジワース展開. Springer. ISBN 9780387945088。
