学生のt
確率密度関数
累積分布関数
パラメータ 自由度実数、ほとんどの場合正の整数
サポート
PDF
CDF

超幾何関数は どこにあるか
平均それ以外の場合定義されていません
中央値
モード
分散それ以外の場合未定義
歪度それ以外の場合定義されていません
過剰尖度それ以外の場合未定義
エントロピ


ここで、はディガンマ関数、はベータ関数である。
MGF未定義
CF

のために

ここで第二種修正ベッセル関数[ 1 ]
予想される不足額

ここでは逆標準化スチューデント t CDF、 は標準化スチューデントt PDFである。[ 2 ]

確率論および統計学においてスチューデントt 分布(または単にt 分布)は、標準正規分布を一般化した連続確率分布です。後者と同様に、ゼロを中心に対称で、ベル型の形状をしています。

しかし、はより重い裾を持ち、裾の確率質量の量はパラメータ によって制御されます。 については、スチューデントt分布は標準的なコーシー分布となり、非常に「太い」裾を持ちます。一方、 については、非常に「細い」裾を持つ標準正規分布となります

「Student」という名前は、ウィリアム・シーリー・ゴセットがアイルランドのダブリンのギネス醸造所で働いていたときに科学論文の発表で使用した仮名です

スチューデントt分布は、2 つの標本平均値の差の統計的有意性を評価するスチューデントt検定、2 つの母集団平均値の差の信頼区間の構築、線形回帰 分析など、広く使用されている多くの統計分析で役割を果たします

位置スケールt 分布 の形式では、正規分布が一般化され、分散パラメータを周辺化するときに複合分布として正規ファミリーのデータのベイズ分析でも発生します。

定義

[編集]

確率密度関数

[編集]

スチューデントt 分布の確率密度関数(PDF)は次式で表されます。ここでは自由度ガンマ関数です。これは と表記されることもあります 。ここではベータ関数です。特に、自由度が整数値の場合、以下の式が成り立ちます。

そしてさらに、

奇数の場合、

確率密度関数は対称であり、その全体的な形状は平均0、分散1の正規分布変数のベル型に似ていますが、やや低く幅が広いという点が異なります。自由度が大きくなるにつれて、t 分布は平均0、分散1の正規分布に近づきます。このため、 t分布は正規性パラメータとも呼ばれます。[ 3 ]

以下の画像は、 t分布の密度 を の値の増加ごとに示しています。正規分布は比較のために青い線で示されています。t分布(赤い線)は が増加するにつれて正規分布に近づくことに注意し ください

 自由度1、2、3、5、10、30におけるt分布(赤)の密度を標準正規分布(青)と比較したグラフ。
以前のグラフは緑で表示されています。
1自由度
2自由度
3自由度
5つの自由度
10自由度
30自由度

累積分布関数

[編集]

累積分布関数(CDF)、正規化不完全ベータ関数Iで表すことができます。t > 0の場合、

どこ

他の値は対称性によって得られる。

ここで、は超幾何関数の特定のインスタンスです

逆累積分布関数の詳細については、分位関数 § スチューデント t 分布を参照してください。

特殊なケース

[編集]

特定の値は、スチューデントの t 分布の単純な形を与えます。

PDFCDFメモ
1コーシー分布を参照
2
3
4
5
正規分布誤差関数を参照


プロパティ

[編集]

瞬間

[編集]

の場合t分布の生のモーメント は

秩序以上の瞬間は存在しない。[ 4 ]

の項(kは偶数)は、ガンマ関数の性質を使って次のように 簡略化できる。

自由度t 分布の場合、期待値は次のようになります。分散は、次のようになります歪度0になります。過剰尖度は、次のようになります。

t分布の発生原理 (特徴づけ)

[編集]

検定統計量の分布として

[編集]

自由度付きスチューデントのt分布は、確率変数Tの分布として定義される。[ 5 ] [ 6 ]

どこ

異なる分布は、与えられた定数 μに対して、 によって 定義される確率変数の分布として定義されます。この確率変数は、非心度パラメータμを持つ非心t分布に従います。この分布は、スチューデントのt検定の検出の研究において重要です

導出
[編集]

X 1 , ..., X n が正規分布に従う確率変数Xの独立した実現であり、期待値がμ分散がσ 2であるとする

標本平均であり、

標本からの分散の不偏推定値となる。確率変数は

は自由度を持つカイ二乗分布に従う(コクランの定理による)。[ 7 ]

標本平均は平均μ、分散σ 2 / nの正規分布に従うため、平均0、分散1の正規分布に従う。さらに、これら2つの確率変数(正規分布に従うZとカイ2乗分布に従うV )は独立であることを示すことも可能である。したがって[説明が必要] 、重要な量

これはZとは異なり、正確な標準偏差σが標本標準誤差sに置き換えられており、上で定義したスチューデントのt分布に従う。未知の母分散σ 2は分子と分母の両方に存在するためTには現れず、打ち消されている点に注意されたい。ゴセットは上記の確率密度関数をn  − 1と等しく直感的に導き出し、フィッシャーは1925年にそれを証明した。 [ 8 ]

検定統計量Tの分布はに依存しますが、μσには依存しません。 μσに依存しないことが、 t分布を理論と実践の両方で重要なものにしているのです

t統計量の標本分布

[編集]

t分布はt統計量 の標本分布 として生じます 。以下では1標本t 統計量について論じますが、対応する2標本t 統計量についてはスチューデントのt検定を参照してください。

不偏分散推定
[編集]

平均と分散を持つ正規分布から独立かつ同一に分布するサンプルを とします。サンプル平均と不偏サンプル分散は次のように与えられます。

結果として得られる(1標本の)t 統計量は次のように与えられる。

自由度 を持つスチューデントのt 分布に従って分布します。

したがって、推論の目的において、平均と分散が未知の母数パラメータである場合、 t統計量は、 t統計量がいずれ にも依存しない確率分布を持つという意味で、 有用な「重要な量」である。

ML分散推定
[編集]

不偏推定値の代わりに、最大尤度推定値を使用して統計量を得る こともできる。これは位置スケールt 分布に従って分布する。

正規分布と逆ガンマ分布の複合分布

[編集]

位置スケールt 分布は、平均と未知の分散を持つガウス分布(正規分布)と、その分散上にパラメータ と を持つ逆ガンマ分布を合成して生成されます。言い換えると、ランダム変数Xは、未知の分散が逆ガンマとして分布するガウス分布を持つと仮定され、その後、分散が周辺化(統合) されます。

同様に、この分布は、パラメータを持つスケーリングされた逆カイ二乗分布とガウス分布を組み合わせた結果であるスケーリングされた逆カイ二乗分布は、逆ガンマ分布とまったく同じ分布であるが、パラメータ化が異なる。すなわち、

この特徴付けが有用な理由は、ベイズ統計学において逆ガンマ分布はガウス分布の分散の共役事前分布であるためである。その結果、位置尺度t 分布は多くのベイズ推論問題において自然に現れる。[ 9 ]

最大エントロピー分布

[編集]

スチューデントのt 分布は、特定の値を持つランダム変数Xの最大エントロピー確率分布です。[ 10 ] [説明が必要] [より良い情報源が必要]これは、pdfが十分な統計量として指数族形式で表すことができるという観察から直接導かれます

スチューデントの確率密度関数の積分とp

[編集]

関数A ( t | ν )は、 t ≥ 0の場合の、スチューデントの確率密度関数f ( t )の-ttの間の  積分です。したがって、観測データから計算された値よりも小さいtの値が偶然に発生する確率を示します。そのため、関数A ( t | ν )は、2 つのデータ セットの平均の差が統計的に有意かどうかをテストするときに使用できます。これは、2 つのデータ セットが同じ母集団から抽出された場合のtの対応する値とその発生確率を計算することにより行います。これはさまざまな状況で使用されますが、特にt 検定で使用されます。自由度がνの統計量tの場合、A ( t | ν )は、2 つの平均が同じである場合にt が観測値よりも小さくなる確率です(小さい方の平均から大きい方の平均を引いてt ≥ 0 とする場合)。これは、 t分布の累積分布関数F ν ( t )から簡単に計算できます 。

ここでI x ( a , b )は正規化された不完全ベータ関数である

統計的仮説検定では、この関数を使用してpを構築します。

[編集]

一般的に

[編集]
  • 非心t 分布は、t分布を一般化し 、非心パラメータを含むようにしたものです。非標準化t 分布とは異なり、非心t分布は対称ではありません(中央値と最頻値は異なります)。
  • 離散スチューデントt 分布は、 rにおける確率質量関数が次式に比例する形で定義される: [ 11 ]ここでabkはパラメータである。この分布は、連続分布に対するピアソン分布と同様の離散分布系の構築から生じる[ 12 ]
  • 正規分布とχ²分布の平方根から変数の比を取ることで、Student A ( t | ν ) サンプルを生成できます正規分布代わりに、例えばアーウィン・ホール分布を使用すると、正規分布、一様分布、三角分布、スチューデント t分布、コーシー分布を含む、全体として対称的な4パラメータ分布が得られます。これは、正規分布の他の対称的な一般化よりも柔軟性に優れています。
  • t分布は比率分布 の一例です
  • t nに分布するランダム変数の二乗は 、スネデコールの F 分布 F 1, nに従って分布します

位置スケールt分布

[編集]

場所規模の変革

[編集]

スチューデントt分布は、位置パラメータ尺度パラメータを導入することで 、3パラメータの位置尺度t 分布に一般化されます。 位置尺度ファミリー変換により次の式が得られます。

結果として得られる分布は、非標準化スチューデントのt 分布とも呼ばれます。

密度と最初の2つのモーメント

[編集]

位置尺度t分布の密度は次のように定義される: [ 13 ]

同様に、密度は次のように表すことができます

このバージョンのディストリビューションの他の特性は次のとおりです。[ 13 ]

特殊なケース

[編集]
  • が位置尺度t 分布に従う場合に対しては平均と分散で正規分布する。
  • 自由度のある位置尺度t 分布はコーシー分布と同等である。
  • との位置スケールt 分布はスチューデントのt 分布に帰着する。

発生と応用

[編集]

頻度主義的統計的推論では

[編集]

スチューデントt 分布は、データが加法誤差を伴って観測される状況において、平均値などの未知のパラメータを推定することを目的とする様々な統計推定問題で用いられます。(ほぼすべての実用的な統計作業と同様に)これらの誤差の母標準偏差が未知であり、データから推定する必要がある場合、この推定から生じる追加の不確実性を考慮するためにt分布がよく用いられます。このような問題のほとんどにおいて、誤差の標準偏差が既知であれば、 t分布 ではなく正規分布が用いられます 。

信頼区間仮説検定は、特定の統計量(例えば標準得点)の標本分布の分位点が必要となる2つの統計手法です。この統計量がデータ線形関数であり、それを標準偏差の通常の推定値で割った場合、得られた値はスチューデントのt 分布に従うように再尺度化および中心化されます。平均値、加重平均値、回帰係数を含む統計分析はすべて、この形式の統計量となります。

教科書の問題では、母集団の標準偏差が既知であるかのように扱われ、スチューデントのt分布を用いる必要がない場合がよくあります。こうした問題は一般的に2種類あります。(1) 標本サイズが非常に大きいため、データに基づく分散 の推定値が確実であるかのように扱える問題、(2) 数学的推論を例示する問題で、著者や講師が解説する点が標準偏差の推定ではないため、標準偏差の推定の問題が一時的に無視される問題です。

仮説検定

[編集]

多くの統計量は、関心のある帰無仮説 の下では中程度の標本数に対してt分布を示すことが示されており、 t分布は有意性検定の基礎となる。例えば、帰無仮説(相関ゼロ)の場合のスピアマンの順位相関係数ρ の分布は、標本数が約20を超える場合のt分布によってよく近似される。 [要出典]

信頼区間

[編集]

Aが次のように選ばれると 仮定する。

Tが自由度n − 1のt 分布に従うとき、対称性により、これはA

つまりAこの確率分布の「95パーセンタイル」であり、

ここでS n は観測値の標本標準偏差である。これは次式と等価である。

したがって、端点が

はμの90%信頼区間です。したがって、正規分布に従うと合理的に期待できる観測値集合の平均を求める場合、t分布を用いて、その平均の信頼限界に理論的に予測される値(例えば帰無仮説 で予測される値)が含まれているかどうかを調べることができます

この結果はスチューデントのt 検定で使用されます。2 つの正規分布からのサンプルの平均の差自体が正規分布するため、t 分布を使用して、その差が合理的にゼロであると想定できるかどうかを調べることができます。

データが正規分布している場合、平均値の片側(1 − α上側信頼限界(UCL)は次の式を使用して計算できます。

結果として得られるUCLは、与えられた信頼区間と母集団サイズにおいて発生する最大の平均値となります。言い換えれば、UCLは観測値集合の平均であるため、分布の平均がUCL 1 − αより小さい確率は、信頼水準1 − αに等しくなります。

予測区間

[編集]

t分布は、平均値と分散が不明な正規分布から、観測されていないサンプルの予測区間 を構築するために使用できます。

ベイズ統計学では

[編集]

スチューデントt 分布、特に3パラメータ(位置スケール)版は、正規分布との関連からベイズ統計学で頻繁に用いられる。正規分布に従う確率変数分散が未知であり、その上に逆ガンマ分布に従う共役事前分布が配置される場合、結果として得られる変数の周辺分布はスチューデントt分布に従う 同じ結果をもたらす同等の構成としては、分散上の共役スケール逆カイ二乗分布、または精度上の共役ガンマ分布が挙げられる。⁠1/σ²​ を分散の上に置くと、t 分布も生じます。これは、正規分布する変数の平均が既知であるか、共役正規分布の事前分布に従って未知であるか、あるいは不適正な定数事前分布に従って未知であるかに関わらず当てはまります。

t 分布が生成される関連する状況は次のとおりです。

堅牢なパラメトリックモデリング

[編集]

t分布 は、正規分布よりも裾が厚いデータに対するモデルとして、正規分布の代替としてよく用いられます。例えば、Lange et al. [ 14 ]を参照してください。従来のアプローチは、外れ値を特定し(例えば、Grubbsの検定を用いて)、何らかの方法でそれらを除外または重み付けを下げることでした。しかし、外れ値を特定することは必ずしも容易ではなく(特に高次元の場合)、t分布はそのようなデータに対する自然なモデル選択であり、ロバスト統計 へのパラメトリックなアプローチを提供します

ベイズ的な説明はGelmanら[ 15 ]に示されています。自由度パラメータは分布の尖度を制御し、尺度パラメータと相関しています。尤度は複数の極大値を持つ可能性があるため、自由度をかなり低い値に固定し、これを既知として他のパラメータを推定することがしばしば必要になります。一部の著者[引用要]は、3から9の間の値がしばしば良い選択であると報告しています。VenablesとRipley [引用要]は、5がしばしば良い選択であると示唆しています。

スチューデントのt 過程

[編集]

実用的な回帰予測のニーズのために、関数に対するスチューデントt分布 の一般化であるスチューデントt過程が導入されました。ガウス過程がガウス分布から構築されるの と同じように、 スチューデントt 過程はスチューデントt分布から構築されますガウス過程の場合、すべての値のセットは多次元ガウス分布を持ちます。同様に、過程)の対応する値が結合多変量スチューデントt 分布を持つ場合、は 区間上のスチューデントt過程です[ 16 ]これらの過程は、回帰、予測、ベイズ最適化、および関連問題に使用されます。多変量回帰と多出力予測のために、多変量スチューデントt 過程が導入され、使用されています。[ 17 ]

選択された値の表

[編集]

以下の表は、片側または両側の臨界領域の範囲における自由度νのt 分布の値を示しています。最初の列はν、上部のパーセンテージは信頼水準、表本体の数値は信頼区間のセクションで説明されている係数です

無限のνを持つ最後の行は、正規分布の臨界点を示します。これは、 無限の自由度を持つt分布が正規分布であるためです。(上記の関連分布を参照)。

一方的な75%80%85%90%95%97.5%99%99.5%99.75%99.9%99.95%
両面50%60%70%80%90%95%98%99%99.5%99.8%99.9%
11.0001.3761.9633.0786.31412.70631.82163.657127.321318.309636.619
20.8161.0611.3861.8862.9204.3036.9659.92514.08922.32731.599
30.7650.9781.2501.6382.3533.1824.5415.8417.45310.21512.924
40.7410.9411.1901.5332.1322.7763.7474.6045.5987.1738.610
50.7270.9201.1561.4762.0152.5713.3654.0324.7735.8936.869
60.7180.9061.1341.4401.9432.4473.1433.7074.3175.2085.959
70.7110.8961.1191.4151.8952.3652.9983.4994.0294.7855.408
80.7060.8891.1081.3971.8602.3062.8963.3553.8334.5015.041
90.7030.8831.1001.3831.8332.2622.8213.2503.6904.2974.781
100.7000.8791.0931.3721.8122.2282.7643.1693.5814.1444.587
110.6970.8761.0881.3631.7962.2012.7183.1063.4974.0254.437
120.6950.8731.0831.3561.7822.1792.6813.0553.4283.9304.318
130.6940.8701.0791.3501.7712.1602.6503.0123.3723.8524.221
140.6920.8681.0761.3451.7612.1452.6242.9773.3263.7874.140
150.6910.8661.0741.3411.7532.1312.6022.9473.2863.7334.073
160.6900.8651.0711.3371.7462.1202.5832.9213.2523.6864.015
170.6890.8631.0691.3331.7402.1102.5672.8983.2223.6463.965
180.6880.8621.0671.3301.7342.1012.5522.8783.1973.6103.922
190.6880.8611.0661.3281.7292.0932.5392.8613.1743.5793.883
200.6870.8601.0641.3251.7252.0862.5282.8453.1533.5523.850
210.6860.8591.0631.3231.7212.0802.5182.8313.1353.5273.819
220.6860.8581.0611.3211.7172.0742.5082.8193.1193.5053.792
230.6850.8581.0601.3191.7142.0692.5002.8073.1043.4853.767
240.6850.8571.0591.3181.7112.0642.4922.7973.0913.4673.745
250.6840.8561.0581.3161.7082.0602.4852.7873.0783.4503.725
260.6840.8561.0581.3151.7062.0562.4792.7793.0673.4353.707
270.6840.8551.0571.3141.7032.0522.4732.7713.0573.4213.690
280.6830.8551.0561.3131.7012.0482.4672.7633.0473.4083.674
290.6830.8541.0551.3111.6992.0452.4622.7563.0383.3963.659
300.6830.8541.0551.3101.6972.0422.4572.7503.0303.3853.646
400.6810.8511.0501.3031.6842.0212.4232.7042.9713.3073.551
500.6790.8491.0471.2991.6762.0092.4032.6782.9373.2613.496
600.6790.8481.0451.2961.6712.0002.3902.6602.9153.2323.460
800.6780.8461.0431.2921.6641.9902.3742.6392.8873.1953.416
1000.6770.8451.0421.2901.6601.9842.3642.6262.8713.1743.390
1200.6770.8451.0411.2891.6581.9802.3582.6172.8603.1603.373
0.6740.8421.0361.2821.6451.9602.3262.5762.8073.0903.291
一方的な75%80%85%90%95%97.5%99%99.5%99.75%99.9%99.95%
両面50%60%70%80%90%95%98%99%99.5%99.8%99.9%
信頼区間の計算

 標本サイズ11、標本平均10、標本分散2の標本があるとします。自由度10で90%の信頼度を得るには、表の片側t値は1.372です。そして、信頼区間は次のように計算されます。

我々は90%の信頼度で真の平均値が以下の値にあると判断した。

つまり、特定のサンプルからこの方法で上限しきい値を計算すると、その上限しきい値が真の平均を超える回数は 90% になります。

そして90%の信頼度で真の平均値は上にあります

言い換えれば、この方法で特定のサンプルから下限しきい値が計算された場合、その下限しきい値は 90% の確率で真の平均を下回ります。

したがって、80%の信頼度(100% − 2 × (1 − 90%) = 80%から計算)では、真の平均は区間内に存在する。

この方法で特定のサンプルから上限と下限の閾値を計算すると、80% の確率で真の平均値が上限閾値を下回り、下限閾値を上回るということは、この方法で計算された特定の上限と下限の閾値の間に真の平均値が位置する確率が 80% であるということとは異なります。「信頼区間」「検察官の誤謬」を参照してください。

現在では、 R プログラミング言語などの統計ソフトウェアや、多くのスプレッドシート プログラムで使用できる関数により、表を使用せずにt分布とその逆分布の値が計算されます 。

計算方法

[編集]

モンテカルロサンプリング

[編集]

スチューデントt 分布からランダム標本を構築するには様々なアプローチがある。問題は、標本が単独で必要なのか、それとも一様標本に分位関数を適用して構築するのか(例えば、コピュラ依存性の多次元応用における)によって異なる[要出典]単独標本抽出の場合、ボックス・ミュラー法とその極形式の拡張は容易に展開できる。[ 18 ]この方法の利点は、自由度νが0に近い場合に他の多くの候補手法が失敗するのに対し、この方法はすべての正の実数に対して同様に適用できる点である。[ 18 ]

歴史

[編集]
統計学者ウィリアム・シーリー・ゴセット(通称「学生」)

統計学において、t分布は1876年にヘルメルト[ 19 ] [ 20 ] [ 21]とリュロス[ 22 ]によって事後分布 として初めて導出されました[ 23 ] [ 24 ]このように、スチューデントのt分布はスティグラーのエポニミー法則の一例です。t分布は、カール・ピアソンの1895年の論文において より一般的な形でピアソンIV型分布として登場しました[ 25 ]

英語文献では、この分布は、ウィリアム・シーリー・ゴセットがアイルランドのダブリンにあるギネス醸造所で働いていた当時、「Student」というペンネームで1908年にバイオメトリカ誌に発表した論文に由来するとされています。[ 26 ]ペンネームの由来に関する一説は、ゴセットの雇用主が、科学論文を発表する際には実名ではなくペンネームを使うことをスタッフに望んでいたため、彼は身元を隠すために「Student」という名前を使ったというものです。また、ギネスが競合他社にt 検定を用いて原材料の品質を判断していることを知られたくなかったという説もあります。[ 27 ] [ 28 ]

ゴセットはギネスで働き、少量サンプルの問題、例えば大麦の化学的性質など、サンプル数が3個程度しかない場合の問題に関心を持っていました。ゴセットの論文では、この分布を「正規母集団から抽出されたサンプルの標準偏差の頻度分布」と呼んでいます。この分布はロナルド・フィッシャーの研究によって広く知られるようになりました。フィッシャーはこの分布を「スチューデント分布」と呼び、検定値を文字tで表しました。[ 8 ] [ 29 ]

参照

[編集]

注記

[編集]
  1. ^ Hurst, Simon. 「スチューデントt分布の特性関数」金融数学研究報告書。統計研究報告書番号SRR044-95。 2010年2月18日時点のオリジナルよりアーカイブ。
  2. ^ Norton, Matthew; Khokhlov, Valentyn; Uryasev, Stan (2019). 「共通確率分布のCVaRとbPOEの計算とポートフォリオ最適化および密度推定への応用」(PDF) . Annals of Operations Research . 299 ( 1– 2). Springer: 1281– 1315. arXiv : 1811.11301 . doi : 10.1007/s10479-019-03373-1 . S2CID 254231768. 2023年2月27日閲覧.  
  3. ^ Kruschke, JK (2015). 『ベイズ統計データ分析の実践』(第2版). Academic Press. ISBN 9780124058880. OCLC  959632184 .
  4. ^ Casella G, Berger RL (1990).統計的推論. ダックスベリー・リソース・センター. p. 56. ISBN 9780534119584
  5. ^ Johnson NL, Kotz S, Balakrishnan N (1995). 「第28章」.連続一変量分布. 第2巻(第2版). Wiley. ISBN 9780471584940
  6. ^ Hogg RV、Craig AT (1978). 『数理統計入門(第4版)』ニューヨーク: Macmillan. ASIN B010WFO0SA . セクション4.4および4.8 {{cite book}}: CS1 maint: postscript (link)
  7. ^ Cochran, WG (1934). 「正規分布系における二次形式の分布と共分散分析への応用」.ケンブリッジ哲学協会数学紀要. 30 (2): 178– 191. Bibcode : 1934PCPS...30..178C . doi : 10.1017/S0305004100016595 . S2CID 122547084 . 
  8. ^ a b Fisher RA (1925). 「『スチューデント分布』の応用」(PDF) . Metron . 5 : 90–104 . 2016年3月5日時点のオリジナル(PDF)からのアーカイブ。
  9. ^ Gelman AB, Carlin JS, Rubin DB, Stern HS (1997).ベイズデータ分析(第2版). ボカラトン, フロリダ州: Chapman & Hal lp 68. ISBN 9780412039911
  10. ^ Park SY, Bera AK (2009). 「最大エントロピー自己回帰条件付き異分散モデル」. J​​ournal of Econometrics . 150 (2): 219– 230. doi : 10.1016/j.jeconom.2008.12.014 .
  11. ^ Ord JK (1972).頻度分布のファミリー. ロンドン, イギリス: Griffin. 表5.1. ISBN 9780852641378
  12. ^ Ord JK (1972).頻度分布のファミリー. ロンドン, イギリス: Griffin. 第5章. ISBN 9780852641378
  13. ^ a b Jackman, S. (2009).社会科学のためのベイズ分析. Wileyシリーズ確率統計. Wiley. p.  507. doi : 10.1002 /9780470686621 . ISBN 9780470011546
  14. ^ Lange KL, Little RJ, Taylor JM (1989). 「 t分布を用いたロバスト統計モデリング(PDF) .アメリカ統計学会誌. 84 (408): 881– 896. doi : 10.1080/01621459.1989.10478852 . JSTOR 2290063 .  
  15. ^ Gelman AB, Carlin JB, Stern HS, et al. (2014). 「計算効率の高いマルコフ連鎖シミュレーション」.ベイジアンデータ分析. フロリダ州ボカラトン: CRC Press. p. 293. ISBN 9781439898208
  16. ^ Shah, Amar; Wilson, Andrew Gordon; Ghahramani, Zoubin (2014). 「ガウス過程の代替としてのスチューデント t過程」(PDF) . JMLR . 33 (第17回国際人工知能統計会議 (AISTATS) 2014, レイキャビク, アイスランド): 877–885 . arXiv : 1402.4306 .
  17. ^ Chen, Zexun; Wang, Bo; Gorban, Alexander N. (2019). 「多変量ガウス分布とスチューデント t過程回帰による多出力予測」 .ニューラルコンピューティングとアプリケーション. 32 (8): 3005– 3028. arXiv : 1703.04455 . doi : 10.1007/s00521-019-04687-8 .
  18. ^ a b Bailey RW (1994). 「t分布によるランダム変数の極生成 」.計算数学. 62 (206): 779– 781. Bibcode : 1994MaCom..62..779B . doi : 10.2307/2153537 . JSTOR 2153537. S2CID 120459654 .  
  19. ^ ヘルマート FR (1875)。 「Über die Berechnung des wahrscheinlichen Fehlers aus einer endlichen Anzahl wahrer Beobachtungsfehler」。Zeitschrift für Angewandte Mathematik und Physik (ドイツ語)。20300~ 303
  20. ^ ヘルマート FR (1876)。 「Uber die Wahrscheinlichkeit der Potenzsummen der Beobachtungsfehler und uber einige damit in Zusammenhang stehende Fragen」。Zeitschrift für Angewandte Mathematik und Physik (ドイツ語)。21 : 192–218 .
  21. ^ ヘルマート FR (1876)。"Die Genauigkeit der Formel von Peters zur Berechnung des wahrscheinlichen Beobachtungsfehlers directer Beobachtungen gleicher Genauigkeit" [同じ精度の直接観測の確率観測誤差を計算するためのピータースの公式の精度]。Astronomische Nachrichten (ドイツ語)。88 ( 8–9 ): 113–132ビブコード: 1876AN....88....113H土井10.1002/asna.18760880802
  22. ^ ルーロス J (1876)。「Vergleichung von zwei Werten des wahrscheinlichen Fehlers」Astronomische Nachrichten (ドイツ語)。87 (14): 209–220ビブコード: 1876AN....87....209L土井10.1002/asna.18760871402
  23. ^ Pfanzagl J, Sheynin O (1996). 「確率統計史研究 XLIV. t 分布の先駆者」. Biometrika . 83 (4): 891– 898. doi : 10.1093/biomet/83.4.891 . MR 1766040 . 
  24. ^ Sheynin O (1995). 「ヘルメルトの誤差理論における研究」.正確科学史アーカイブ. 49 (1): 73– 104. doi : 10.1007/BF00374700 . S2CID 121241599 . 
  25. ^ ピアソン, K. (1895). 「進化の数学的理論への貢献. II. 均質物質における歪変化」 .王立協会哲学論文集A: 数学・物理学・工学. 186 (374): 343– 414. Bibcode : 1895RSPTA.186..343P . doi : 10.1098/rsta.1895.0010 . ISSN 1364-503X . 
  26. ^ 「学生」[ウィリアム・シーリー・ゴセットの偽名] (1908). 「平均値の推定誤差」(PDF) . Biometrika . 6 (1): 1– 25. doi : 10.1093/biomet/6.1.1 . hdl : 10338.dmlcz/143545 . JSTOR 2331554 .  {{cite journal}}: CS1 maint: numeric names: authors list (link)
  27. ^ Wendl MC (2016). 「偽名の名声」. Science . 351 (6280): 1406. Bibcode : 2016Sci...351.1406W . doi : 10.1126/science.351.6280.1406 . PMID 27013722 . 
  28. ^ Mortimer RG (2005).物理化学のための数学第3版). バーリントン, マサチューセッツ州: エルゼビア. pp.  326. ISBN 9780080492889. OCLC  156200058 .
  29. ^ Walpole RE, Myers R, Myers S, Ye K (2006).エンジニアと科学者のための確率と統計(第7版). ニューデリー, インド: ピアソン. p. 237. ISBN 9788177584042. OCLC  818811849 .
  30. ^ 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 .  

参考文献

[編集]
[編集]