切断正規分布

確率密度関数
異なるパラメータセットにおける切断正規分布の確率密度関数。すべてのケースにおいて、a = −10、b = 10です。黒:μ = −8、σ = 2、青:μ = 0、σ = 2、赤:μ = 9、σ = 10、オレンジ:μ = 0、σ = 10。
累積分布関数
さまざまなパラメーターのセットに対する切り捨て正規分布の累積分布関数。すべての場合において、a = −10 およびb = 10 です。黒の場合: μ = −8、σ = 2。青: μ = 0、σ = 2。赤: μ = 9、σ = 10。オレンジ色: μ = 0、σ = 10。
表記
パラメータ
(ただし定義を参照)— の最小値— の最大値

サポート
PDF[1]
CDF
平均
中央値
モード
分散
エントロピ
MGF

確率論と統計学において、切断正規分布とは、正規分布に従う確率変数の確率分布において、その確率変数を下方または上方(あるいはその両方)から制限することによって導かれる分布である。切断正規分布は、統計学や計量経済学において広く応用されている

定義

が平均と分散 の正規分布に従い、区間 内にあると仮定します。この場合、 の条件付きで、は切断正規分布に従います。

確率密度関数はに対して、次のように与えられる。

およびそれ以外の場合。

ここで、 は標準正規分布の確率密度関数でありはその累積分布関数です 。定義により、 の場合、 となり、同様に、の場合、 となります

上記の式は、切断された正規分布の尺度パラメータが負の値をとることが許される場合を示しています。この場合、パラメータは虚数ですが、関数は実数で正であり、正規化可能です。切断されていない正規分布の尺度パラメータは正でなければなりません。そうでなければ、分布は正規化できないからです。一方、二重に切断された正規分布は、そのような積分可能性の問題は有界領域では発生しないため、原理的には負の尺度パラメータ(分散とは異なります。まとめの式を参照)を持つことができます。この場合、分布は を条件とする切断されていない正規分布として解釈することはできませんが、第 1 および第 2 モーメントを制約条件とする最大エントロピー分布として解釈することはでき、さらにおよびに位置する 1 つではなく2 つの極大値を示すという特異な特徴があります

プロパティ

切断正規分布は、平均と分散が固定され、区間[a,b]に制約された2つの最大エントロピー確率分布のうちの1つであり、もう1つは切断U分布である。[2]固定サポートを持つ切断正規分布は指数族を形成する。ニールセン[3] は、最初の分布のサポートが2番目の分布のサポートにネストされた2つの切断正規分布間のカルバック・ライブラー距離とバタチャリヤ距離を計算するための閉形式の公式を報告した。

瞬間

確率変数が下側からのみ切り捨てられた場合、確率質量の一部がより高い値にシフトし、一次確率的に支配的な分布となり、平均は元の正規分布の平均よりも高い値に増加します。同様に、確率変数が上側からのみ切り捨てられた場合、切り捨てられた分布の平均は以下になります。

ランダム変数が上、下、またはその両方に制限されているかどうかに関係なく、切り捨ては平均を保存する縮小と平均を変更する固定シフトを組み合わせたものなので、切り捨てられた分布の分散は元の正規分布の分散よりも小さくなります。

両側切り捨て[4]

とします。すると、

これらの式の数値評価には注意が必要です。区間にが含まれていない場合、破滅的な相殺が生じる可能性があります。この問題を回避するために、式を書き直すより良い方法があります。[5]

片側切断(下尾部)[6][7]

この場合

そして

どこ

片側切断(上尾部)

この場合

Barr & Sherrill (1999) は、片側切断の分散をより簡略化した式を提示しています。彼らの式は、標準的なソフトウェアライブラリに実装されているカイ二乗累積分布関数に基づいています。Bebu & Mathew (2009) は、切断モーメントの(一般化)信頼区間の式を提示しています。

再帰式

切り捨てられていない場合には、切り捨てられたモーメントに対する再帰式が存在する。[8]

特に、 については、

証拠

変数の変更により、次式が得られます。部分積分 を使用すると、証明する方程式が得られます。

多変量

多変量切断正規分布のモーメントを計算するのは困難です。

切断正規分布から値を生成する

として定義される確率変量は、正規分布の累積分布関数(つまり、正しい平均と分散を持つ)からサンプリングされ、その逆関数である 上の一様乱数は、範囲 に切り捨てられた分布に従います。これは、確率変数をシミュレートするための逆変換法です。この方法は最も単純な方法の1つですが、正規分布の裾野でサンプリングすると失敗する可能性があり、[9]、非常に遅くなることもあります。[10]そのため、実際には、代替のシミュレーション方法を見つける必要があります。

そのような切り捨て正規分布生成器の1つ(MatlabおよびR(プログラミング言語)でtrandn.Rとして実装)は、Marsagliaによる受け入れ拒否の考えに基づいています。[11] Marsaglia(1964)の受け入れ率はRobert(1995)と比較してわずかに最適ではありませんが、Marsagliaの方法は通常より高速です。[10]これは、指数関数のコストのかかる数値評価を必要としないためです。

切断正規分布からの抽出のシミュレーションの詳細については、Robert (1995)、Lynch (2007、セクション8.1.3 (200~206ページ))、Devroye (1986) を参照してください。RのMSMパッケージには、切断正規分布からの抽出を計算する関数rtnormが含まれています。Rのtruncnormパッケージにも、切断正規分布からの抽出を行う関数が含まれています。

Chopin (2011) は、Marsaglia と Tsang (1984, 2000) のZiggurat アルゴリズムに着想を得たアルゴリズムを提案しました (arXiv) 。このアルゴリズムは、通常、最速のガウスサンプラーと考えられており、Ahrens のアルゴリズム (1995) にも非常に近いものです。実装は C、C++、Matlab、Python で見つけることができます。

多変量切断正規分布からのサンプリングは、はるかに困難です。[12]正確な、あるいは完全なシミュレーションは、正規分布を多面体領域に切断した場合にのみ可能です。[12] [13]より一般的なケースとして、Damien & Walker (2001) は、ギブスサンプリングの枠組みの中で切断密度をサンプリングするための一般的な手法を導入しています。彼らのアルゴリズムは1つの潜在変数を導入し、ギブスサンプリングの枠組みの中で、Robert (1995) のアルゴリズムよりも計算効率に優れています。

参照

注記

  1. ^ 「講義 4: 選択」(PDF) . web.ist.utl.pt高技術研究所。 2002 年 11 月 11 日。p. 1 . 2015 年7 月 14 日に取得
  2. ^ Dowson, D.; Wragg, A. (1973年9月). 「規定された第一モーメントと第二モーメントを持つ最大エントロピー分布(対応)」. IEEE Transactions on Information Theory . 19 (5): 689– 693. doi :10.1109/TIT.1973.1055060. ISSN  1557-9654.
  3. ^ Frank Nielsen (2022). 「ネストされたサポートを持つ切断指数族の密度間の統計的乖離:デュオ・ブレグマン・ダイバージェンスとデュオ・ジェンセン・ダイバージェンス」.エントロピー. 24 (3). MDPI: 421. Bibcode :2022Entrp..24..421N. doi : 10.3390 /e24030421 . PMC 8947456. PMID  35327931. 
  4. ^ ジョンソン、ノーマン・ロイド、コッツ、サミュエル、バラクリシュナン、N. (1994).連続一変量分布. 第1巻(第2版). ニューヨーク: ワイリー. 第10.1節. ISBN 0-471-58495-9. OCLC  29428092。
  5. ^ Fernandez-de-Cossio-Diaz, Jorge (2017-12-06), TruncatedNormal.jl: 単変量切断正規分布の平均と分散を計算する(ピークから離れたところで動作する)2017-12-06取得
  6. ^ グリーン、ウィリアム・H. (2003).計量経済分析(第5版). プレンティス・ホール. ISBN 978-0-13-066189-0
  7. ^ デル・カスティージョ、ジョアン (1994年3月). 「単一切断正規分布:非急峻な指数分布族」(PDF) .統計数学研究所紀要. 46 (1): 57– 66. doi :10.1007/BF00773592.
  8. ^ エリック・オルジェビンによる文書、「https://people.smp.uq.edu.au/YoniNazarathy/teaching_projects/studentWork/EricOrjebin_TruncatedNormalMoments.pdf」
  9. ^ Kroese, DP ; Taimre, T.; Botev, ZI (2011).モンテカルロ法ハンドブック. John Wiley & Sons.
  10. ^ ab Botev, ZI; L'Ecuyer, P. (2017). 「正規分布の裾の区間を切り捨てたシミュレーション」.第10回EAI国際パフォーマンス評価方法論およびツール会議. 2016年10月25日~28日 タオルミーナ(イタリア): ACM. pp.  23– 29. doi :10.4108/eai.25-10-2016.2266879. ISBN 978-1-63190-141-6{{cite conference}}: CS1 maint: location (link)
  11. ^ Marsaglia, George (1964). 「正規分布の裾からの変数生成」. Technometrics . 6 (1): 101– 102. doi :10.2307/1266749. JSTOR  1266749.
  12. ^ ab Botev, ZI (2016). 「線形制約下における正規則:ミニマックス傾斜法によるシミュレーションと推定」. Journal of the Royal Statistical Society, Series B. 79 : 125–148 . arXiv : 1603.04166 . doi : 10.1111/rssb.12162. S2CID  88515228.
  13. ^ Botev, Zdravko & L'Ecuyer, Pierre (2018). 「第8章:単変量および多変量正規分布の裾からのシミュレーション」. Puliafito, Antonio (編).システムモデリング:方法論とツール. EAI/Springer Innovations in Communication and Computing . Springer, Cham. pp.  115– 132. doi :10.1007/978-3-319-92378-9_8. ISBN 978-3-319-92377-2. S2CID  125554530。
  14. ^ Sun, Jingchao; Kong, Maiying; Pal, Subhadip (2021年6月22日). 「修正半正規分布:特性と効率的なサンプリング手法」. Communications in Statistics - Theory and Methods . 52 (5): 1591– 1613. doi :10.1080/03610926.2021.1934700. ISSN  0361-0926. S2CID  237919587.

参考文献

  • Botev, Zdravko & L'Ecuyer, Pierre (2018). 「第8章:単変量および多変量正規分布の裾からのシミュレーション」. Puliafito, Antonio (編).システムモデリング:方法論とツール. EAI/Springer Innovations in Communication and Computing. Springer, Cham. pp.  115– 132. doi :10.1007/978-3-319-92378-9_8. ISBN 978-3-319-92377-2. S2CID  125554530。
  • Devroye, Luc (1986). 非一様乱数変量生成(PDF) . ニューヨーク: Springer-Verlag. オリジナル(PDF)から2014年8月18日にアーカイブ。 2012年4月12日閲覧
  • グリーン、ウィリアム・H. (2003).計量経済分析(第5版) . プレンティス・ホール. ISBN 978-0-13-066189-0
  • ノーマン・L・ジョンソン、サミュエル・コッツ (1970).連続一変量分布-1、第13章. John Wiley & Sons.
  • リンチ、スコット(2007年)『社会科学者のための応用ベイズ統計と推定入門』ニューヨーク:シュプリンガー、ISBN 978-1-4419-2434-6
  • Robert, Christian P. (1995). 「切断正規分布変数のシミュレーション」.統計と計算. 5 (2): 121– 125. arXiv : 0907.4010 . doi :10.1007/BF00143942. S2CID  15943491.
  • Barr, Donald R.; Sherrill, E.Todd (1999). 「切断正規分布の平均と分散」. The American Statistician . 53 (4): 357– 361. doi :10.1080/00031305.1999.10474490.
  • Bebu, Ionut; Mathew, Thomas (2009). 「正規分布と対数正規分布における有限モーメントと切断モーメントの信頼区間」.統計と確率論文集. 79 (3): 375– 380. doi :10.1016/j.spl.2008.09.006.
  • ダミアン, ポール; ウォーカー, スティーブン G. (2001). 「サンプリングによる切断正規分布、ベータ分布、ガンマ分布の密度」.計算およびグラフィカル統計ジャーナル. 10 (2): 206– 215. doi :10.1198/10618600152627906. S2CID  123156320.
  • Chopin, Nicolas (2011-04-01). 「切断ガウス分布の高速シミュレーション」. Statistics and Computing . 21 (2): 275– 288. arXiv : 1201.6140 . doi :10.1007/s11222-009-9168-1. ISSN  1573-1375.
  • Burkardt, John. 「切断正規分布」(PDF) .フロリダ州立大学科学計算学部ウェブサイト. 2018年2月15日閲覧.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Truncated_normal_distribution&oldid=1315201905"