確率不等式


37

無制限のランダム変数の合計の確率不等式を探しています。誰かが私にいくつかの考えを提供できるなら、本当に感謝しています。

私の問題は、実際には2つのiidガウスの乗算である無制限のiid確率変数の合計が特定の値、つまりを超える確率に関する指数の上限を見つけることです。、、とからIIDが生成される。Pr[X≥ϵσ2N]≤exp⁡(?)X=∑i=1NwiviwiviN(0,σ)

モーメント生成関数(MGF)を使用してChernoff境界を使用しようとしましたが、派生境界は次のようになります。

Pr[X≥ϵσ2N]≤minsexp⁡(−sϵσ2N)gX(s)=exp⁡(−N2(1+4ϵ2−1+log⁡(1+4ϵ2−1)−log⁡(2ϵ2)))

ここで、gX(s)=(11−σ4s2)N2はXのMGFですX。しかし、限界はそれほど厳しくありません。私の問題の主な問題は、ランダム変数が制限されていないことであり、残念ながら、ヘフディング不等式の境界を使用することはできません。

あなたが私にいくつかのきつい指数関数的境界を見つけるのを手伝ってくれれば幸いです。


3
圧縮センシング関連の問題のように聞こえます。R. Vershyninの非漸近ランダム行列理論、特に彼がサブ指数ランダム変数と呼ぶものの限界に関するノートを調べてください。それで始められます。さらにポインタが必要な場合はお知らせください。さらに情報を投稿しようと思います。
— 枢機

1
math.SEのこのトピックには、少なくとも2つの関連する質問と回答があります(免責事項:私が参加したものを含む)。
— 枢機

1
製品wiviは、「正規の製品」分布があります。この積の平均はゼロで、分散はσ4で、σ2はwiとv_iの分散ですvi。以下のためにN大きめの、あなたはのおおよそのnorality取得するために中心極限定理を使用することができますX。正規の製品分布のスキューを計算できる場合、Berry-Esseenの定理を適用してCDFの収束速度を制限できると思います。
— みすぼらしいシェフ

1
@ shabbychef、Berry-Esseenは収束がかなり遅いです。これは、すべての分布関数クラスに対する一様な境界だからです。F
— 枢機卿

4
@DilipSarwate:しばらく前にコメントが表示されてすみません。math.SEでも何度かリンクした次の小さな論文に興味があるかもしれないと思います:TK Phillips and R. Nelson(1995)、モーメントバインドはChernoffのポジティブテールのバインドよりもタイトです確率、アメリカ統計学、第42巻、No。2.、175-178。
— 枢機

回答:


1

後で指定されるいくつかに対して提案したChernoff境界を使用して、 おかげで2番目の不等式が成り立つ任意の。ここで、および取る と、右側はこれは 、任意のに対して。s≤1/(2σ2)− log (1 − x )

P[X>t]≤exp⁡(−st)exp⁡(−(N/2)log⁡(1−σ4s2))≤exp⁡(−st+σ4s2N)
−log⁡(1−x)≤2xx∈(0,1/2)t=ϵσ2Ns=t/(2σ4N)exp⁡(−t2/(4σ4N)=exp⁡(−ϵ2N/4)ε ∈ (0 、1 )
P[X>ϵσ2N]≤exp⁡(−ϵ2N/4).
ϵ∈(0,1)

別の手段は、ハンソン・ライトの不等式などの濃度不等式、または関心のあるランダム変数を含む次数2のガウスカオスの濃度不等式を直接適用することです。

モーメント生成機能を使用しないシンプルなアプローチ

取る(そうでなければ、いずれかで割ることによってスケール変更することができる簡単にするために)。σ=1σ2

ライトと。上限を求めています。v=(v1,...,vn)Tw=(w1,...,wn)TP(vTw>ϵN)

してみましょう。その場合、 および独立性によるは、自由度の分布の独立しています。Z=wTv/‖v‖Z∼N(0,1)v,w‖v‖2Zχ2n

標準正規および確率変数の標準境界、 ユニオン境界と組み合わせると、 にの形式の上限が得られます。χ2

P(|Z|>ϵn/2)≤2exp⁡(−ϵ2n/4),P(‖v‖>2n)≤exp⁡(−n(2−1)2/2).
P(vTw>ϵN)2exp⁡(−ϵ2n/4)+exp⁡(−n(2−1)2/2)


0

取得する境界は、順序としてです。一般的な方がずっといいとは思わない。製品変数のウィキペディアのページから、の分布はで、は修正ベッセル関数です。DLMF関数リストの(10.25.3)から、なので、十分に大きいは、サブガウス境界を与えません。e−ϵϵ→∞ϵwはiは、V 、I K 0(Z )/ π K 0 K 0(T )〜E - T / √wiviK0(z)/πK0K0(t)∼e−t/txP(wivi>x)∼∫x∞e−t/tdt

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