T分布確率変数の二乗和の分布


11

T分布確率変数の二乗和の分布と、テール指数調べています。Xがrvである場合、フーリエ変換であるは、畳み込み前の正方形の解を与えます。 αX2F(t)F(t)n

F(t)=0exp(itx2)((αα+x2)α+12α B(α2,12))dx

で、溶液をフーリエ逆変換を行うと逆することができるが、扱いにくいと不可能である。したがって、問題は、標本分散の分布またはT分布確率変数の標準偏差で作業が行われたかどうかです。(Gaussianのカイ二乗とは、StudentTにとってのものです)。ありがとうございました。Ft nα=3F(t)n

(考えられる解決策)はフィッシャー分布であるため、フィッシャー分布変数の合計を確認します。 F 1 α X2F(1,α)

(可能な解決策)特性関数から、合計された平均 は、分布の最初の2つのモーメントが存在する場合、それらのモーメントは同じです。したがって、uを平方根とし、確率分布内で変数を変更すると、n個のサンプルのT変数の標準偏差の密度は、次のように近似できます X 2 F N α G U = 2 α α / 2 N N / 2 UはN - 1 α + N U 2 - αnX2F(n,α)

g(u)=2αα/2nn/2un1(α+nu2)α2n2B(n2,α2)

1
F F 1 α T2は分布です。独立した分布変数の合計の平均と分散は容易に導き出されますが、分布は閉じた形式では利用できません。詳細については、この質問を参照してください。リンクされた論文が役に立つかもしれません。特性関数は、FのWikipediaページにも記載されています。[t分布変数の標本分散は、かなり異なる質問です。]FF(1,α)
Glen_b -Reinstate Monica

回答:


7

あなたの質問の明確化(私には2つの関連しているが異なる部分があるようです):あなたは(1)独立した二乗確率変数の合計の分布、および(2)サンプリングを探しています分布から抽出されたランダムサンプルの分散(または関連する標準偏差)の分布(おそらく(1)について尋ねる理由)。T αトンαn tαtα

独立二乗変数の合計の分布tα

もし(独立した)ランダムでの変数が DF、それが偽である(これは2番目の「可能な解決策」で主張しているように見えます)。これは、それぞれの最初の瞬間を考慮することで簡単に検証できます(後者の最初の瞬間は最初の瞬間の倍です)。 Titαtαi=1nTi2F(n,α)n

最初の「可能な解決策」の主張は正しいです:。特性関数に頼るのではなく、比率分布としての分布の特徴付けを検討するとき、この結果はより透明であると思いますここで、は標準正規変数です。そして有するカイ自乗変数であるの独立した自由度、。この比率の2乗は、2つの独立したカイ2乗変数の比率であり、それぞれの自由度でスケーリングされます。つまり、Ti2F(1,α)tZU/αZUαZV/1U/αV=Z2これは、分布の標準的な特性です(分子dfは1に等しく、分母dfは等しい)。F(1,α)α

上記の最初の段落の最初の瞬間に書いたメモを考えると、より良い主張は [I haveここでは、分布に同じ式を使用し、その分布を持つ確率変数を使用して、表記を少し乱用しています。]。最初のモーメントは一致しますが、2番目の中心モーメントは一致しません(の場合、最初の式の分散は、後者の式の分散よりも小さくなります)。したがって、この主張も誤りです。[とはいえ、であることを観察するのは興味深いことです。これは、2乗(標準)を合計したときに期待される結果です。通常の変量。]i=1nTi2nF(n,α)α>4limαnF(n,α)=χn2

分布からサンプリングするときの分散のサンプリング分布tα

上で書いたことを考えると、「n個のサンプルのT変数の標準偏差の密度」に対して得られる式は正しくありません。ただし、が正しい分布であったとしても、標準偏差は二乗和の平方根ではありません(密度に到達するために使用していたようです)。代わりに、の(スケーリングされた)サンプリング分布を探し。通常の場合、この式のLHSは、2乗された正規変数の合計として書き直すことができます(2乗内の項は、正規分布の正規結合の線形結合として書き直すことができます)。おなじみF(n,α)g(u)i=1n(TiT¯)2=i=1nTi2nT¯2トントンχ2分布。残念ながら、変数の線形結合(自由度が同じであっても)はとして分布されないため、同様のアプローチは利用できません。tt

おそらく、何を実証したいのかを再考する必要がありますか?たとえば、いくつかのシミュレーションを使用して目的を達成できる可能性があります。ただし、した例を示しています。これは、の最初のモーメントのみが有限であるため、シミュレーションはそのようなモーメントの計算に役立ちません。 F 1 α α=3F(1,α)


マークありがとうございます。実際、最初の2つの瞬間は保持されますが、畳み込みは失敗します。カイ二乗を試みて元に戻します。
Nero

質問を言い換えました。または、ページの他の場所に変更を投稿する必要がありますか?
Nero

Nero-質問への変更は質問に表示されます。それが役立つ場合は、いつでも質問の中で質問がどのように変更されたかを通知できます(ただし、必要に応じて質問と回答の編集履歴全体が利用可能であることを覚えておいてください)。
Glen_b-2015

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.