カイ2乗変量の無限コレクションの統計(たとえば、最小)を注文しますか?


11

今回が初めてなので、フォーマットやタグなど、何らかの形で質問を明確にできるかどうか教えてください。(うまくいけば、後で編集できます!)参照を見つけて、誘導を使用して自分自身を解決しようとしましたが、両方で失敗しました。

私は、自由度の異なる独立した確率変数の無数に無限の集合の次数統計に減少するように見える分布を単純化しようとしています。具体的には、独立した中で番目に小さい値の分布は何ですか? m個χ 2 2χ 2 4χ 2 6χ 2 8...χ2mχ22,χ42,χ62,χ82,

特別なケース興味があります:(独立した)の最小値の分布はどうですか?χ 2 2χ 2 4χ 2 6...m=1χ22,χ42,χ62,

最小値の場合、累積分布関数(CDF)を無限積として書くことができましたが、それをさらに単純化することはできません。のCDF が (場合、これにより、期待値2の指数分布との等価性に関する以下の2番目のコメントが確認されます。)最小値のCDFは、として記述できます。 製品の最初の項はであり、「最後の」項は F 2 MX = γ M X / 2 / Γ M = γ M X / 2 /M - 1 = 1 e x / 2 m 1 k = 0 x k /2 kχ2m2M = 1 Fは、mは、I NX = 1 - 1 - F 2X 1 - F 4X ... = 1 - Π M = 11 - F 2 MX = 1 - Π M =

F2m(x)=γ(m,x/2)/Γ(m)=γ(m,x/2)/(m1)!=1ex/2k=0m1xk/(2kk!).
m=1
Fmin(x)=1(1F2(x))(1F4(x))=1m=1(1F2m(x))
=1m=1(ex/2k=0m1xk2kk!).
ex/2ex/2k=0xk/(2kk!)=1です。しかし、そこから単純化する方法(可能な場合)がわかりません。あるいは、まったく異なるアプローチの方が良いかもしれません。

もう1つの潜在的に役立つ注意:は期待値2の指数分布と同じであり、はそのような2つの指数の合計などです。χ22χ42

誰かが興味を持っている場合は、定数の回帰(すべてのに対して)の場合について、このペーパーの定理1を簡略化しようとしています。(を掛けたので、分布の代わりにます。)xi=1iχ2Γ2κ


これは、あなたの質問に答えますか?
mpiktas 2011

@mpiktas:提案をありがとう。それは似ていますが、レートパラメーターが異なる指数の代わりに、自由度が異なるカイ2乗があります(有限ではなく、無限の数です)。また、は指数関数ですが、は指数関数ではありません。それらは指数の合計ですが、指数の合計自体は指数ではありません。(そして、理想的には、私は一般的な注文統計を期待していますが、最小値は素晴らしいスタートです。)χ22χ42,χ62,
デビッドMカプラン

1
これには閉じたフォームがあるとは思えません。ただし、奇妙な特徴がありますが Poisson()変量である場合、であり、はすべての。Xkλ/2k=1,2,1Fmin(λ)Xkk
whuber

1
@whuber:ポアソンプロセスの観点から考えると、それほど興味深いことではないかもしれません。LET BE IID対応ポアソン過程とランダム変数、レート。ましょう、、、等そして、独立したポアソン過程の定常独立インクリメントプロパティであり、我々そのます。T1,T2,Exp(1/2)N(t):=sup{n:i=1nTit}1/2U1=T1U2=T2+T3U3=T4+T5+T6Uiχ2i2P(Uit)=P(N(t)i)
枢機卿

@Cardinalもちろん:それはそれを見る良い方法です。好奇心はポワソンとガンマの関係ではありません。イベント自体の説明にあります!
whuber

回答:


8

無限積のゼロは、項のゼロの和集合になります。20項まで計算すると、一般的なパターンがわかります。

複素数ゼロのプロット

この複素平面のゼロのプロットは、さまざまな記号によって積の個々の項の寄与を区別します。各ステップで、見かけの曲線がさらに拡張され、新しい曲線がさらに左に開始されます。

この図の複雑さは、Whittakerのような古典的なテキストで調査されるように、高次分析のよく知られた関数(ガンマ、シータ、超幾何関数など)と基本関数の観点から閉じた形式のソリューションが存在しないことを示しています&ワトソン)。

したがって、問題はもう少し実り多い形でもたらされる可能性があります。注文統計の分布について何を知る必要がありますか?それらの特徴的な機能の推定値?低次モーメント?分位数への近似?他に何か?


製品のゼロが重要なのはなぜですか?些細なことを見逃しているように感じます。
mpiktas 2011

2
@mp零点と極は、関数の複雑さについて何かを示しています。有理関数にはそれらの有限数があります。初等関数は通常、に対して、積分などのゼロの行を持ちます。典型的な「超越」関数は、すべての非正の整数(ガンマ関数の逆数)や点の格子(シータ関数と楕円関数)など、やや複雑なゼロのパターンを持っています。ここに示されている複雑なパターンは、CDFをこれらの使い慣れた関数の観点から表現することが困難または不可能であることを示唆しています。2iπnnexp()
whuber

2
@whuber(1/2)、ありがとう!複素平面にゼロのそれらの異なるパターンを持つ関数の異なるクラスについては知りませんでした。それは非常に便利に聞こえ、あなたのグラフは私の質問に答えているようです(提起されたとおり)。
デビッドMカプラン

@whuber(2/2)、これは別の論文で与えられた推定量の(複雑な)分布の特別なケースをチェックしていました。彼らは、ブートストラップの使用を正当化するためにディストリビューションの存在を使用しました。私の顧問は、分布を概算しようと提案しました。この特別なケース(私はそれが何であるかを知っている)では、彼らの配布がオフになっているようですので、助成金の締め切り後にアドバイザーと一緒にチェックします。しかし、潜在的に、より複雑な設定で番目の次数統計(除算)の高次展開をとして取得しようとしています。その場合は再度投稿します。再度、感謝します!mmm
デビッドMカプラン

4

(独立した)の最小値の分布は何ですか?χ22,χ42,χ62,

約6年遅れで到着したことをお詫びします。OPはおそらく他の問題に移行しているようですが、問題はまだ新鮮で、別のアプローチを提案するかもしれません。


where where with pdf'sが与えられ:(X1,X2,X3,)XiChisquared(vi)vi=2ifi(xi)

ここに画像の説明を入力してください

以下は、対応するpdfののプロットです。サンプルサイズが増加するにつれて、:fi(xi)i=1 to 8

ここに画像の説明を入力してください

の分布に興味があり。min(X1,X2,X3,)

追加の用語を追加するたびに、追加された最後の最後の用語のpdfはさらに右にシフトします。そのため、用語を追加する効果は、ますます関連性が低くなるだけでなく、ほんの数用語後になります。 、サンプルの最小値では、ほとんど無視できます。これは、実際には、非常に少数の用語のみが実際に問題になる可能性が高いことを意味します。用語の追加(または無数の用語の存在)は、サンプルの最小問題にはほとんど関係ありません。

テスト

これをテストするために、のpdf を1項、2項、3項、4項、5用語、6用語、7用語、8用語に計算しました。 9項、および10項。これを行うために、私はmathStaticaの関数を使用して、ここでサイズのサンプルの最小値(オーダー統計)のpdfを計算するように指示し、パラメーター(代わりに修正される)は:1 st j i v imin(X1,X2,X3,)OrderStatNonIdentical1stjivi

ここに画像の説明を入力してください ここに画像の説明を入力してください

項の数が増えると少し複雑になりますが、上記の1項(1行目)、2項(2行目)、3項(3行目)、4項の出力を示しました。

次の図は、サンプル最小値のpdfを1項(青)、2項(オレンジ)、3項、および10項(赤)と比較しています。結果が3項vs 10項でどれほど似ているかに注意してください。 ここに画像の説明を入力してください

次の図は、5つの用語(青)と10の用語(オレンジ)を比較しています。プロットは非常に類似しており、互いに消去され、違いを確認することもできません。

ここに画像の説明を入力してください

つまり、項の数を5から10に増やしても、サンプルの最小値の分布に視覚的な影響はほとんどありません。

半ロジスティック近似

最後に、サンプルminのpdfの優れた単純な近似は、pdfを使用した半ロジスティック分布です。

g(x)=2ex(ex+1)2 for x>0

次の図は、正確な解と10項(5項または20項との区別がつかない)と半ロジスティック近似(破線)を比較しています。

ここに画像の説明を入力してください

20項に増やしても、識別できる違いはありません。

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