ベルヌーイを証明することはベータの限界です


7

それは私たちが解決した場合、その検査によって私には明らかだ(それによって、平均固定)としましょう(、ベータ分布はベルヌーイに近づく)配布。β=1μμαα0μ

例えば:

par(mfrow = c(1, 2),
    oma = c(0, 0, 1.5, 0))
xx = seq(0, 1, length.out = 1000)
mus = c(.2, .7)
for (ii in 1:2) {
  mu = mus[ii]
  matplot(xx, sapply(10^(-1:-5), function(al) 
    pbeta(xx, al, (1-mu)/mu * al)),
    type = 'l', lty = ii,
    main = paste('Mean:', mu),
    ylab = 'Cumulative Probability', xlab = 'x')
}
title('Beta Approaches Bernoulli', outer = TRUE)

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

つまり、場合、そのCDF、XB(α,1μμα)FX(x;α,μ)

limα0FX(x;α,μ)=1μ=FY(x)x(0,1)

ここで、ベルヌーイ()、および収束が不均一です。私はこれをより正式に証明しようとしましたが、単純なケース(ような)でさえ、何も進歩することができませんでした。ベータ版ウィキペディアのページでは、ケースが証明されずに何度か言及されています。おそらく、積分を行うことについて簡単なこと、またはベータ/不完全ベータ関数についてのいくつかの事実が欠けています。何か助けていただければ幸いです!Yμμ=12μ=12


1
stats.stackexchange.com/a/237849/919での私の回答は、限界確率がと集中することを示しているため、これを暗示しています。図については、2番目のグラフィックスのペアを参照してください。01
whuber

回答:


7

分析は、特に完全で基本的な厳密さで実行される場合、少し面倒になる可能性がありますが、アイデアはシンプルで簡単に把握できます。とに非常に近い小さな領域に焦点を当てます。及びアプローチ、ベータのほぼすべての確率分布は、これらの領域内に位置するようになります。領域のサイズを縮小すると、限界分布が存在する場合、それはベルヌーイ分布のみになる可能性があることがわかります。質問で説明されているように、比率定数に近づけることによってのみ、制限分布を作成できます。01αβ0(α,β)α:β

この分析の良いところは、相対領域を見ると、正規化定数であるベータ関数の動作を考慮する必要がなくなることです。これはかなり単純化されています。(ベータ関数の回避は、ベータ分布の2つの分位点でパラメーターを決定するか?でのベータ分布分位点の私の分析と精神的に似ていますB(α,β)

この分析のさらなる特徴は、定数形式の単純な積分によって不完全ベータ関数を近似することです。これにより、すべてが微積分の最も基本的な演算と代数的不等式に削減されます。tcdtc>1


ベータ版PDFは、比例し検討小と下領域への寄与を調べる 3つの間隔以内、、及びとして及びは小さくなります(ただし、プラスのままです)。

f(x)=xα1(1x)β1.
ϵ>0f(0,ϵ](ϵ,1ϵ)[1ϵ,1)αβ

最終的には両方のと以下になります両方の:したがって、両方の極を持つことになりますとのように見て、:αβ1f01

図1

のグラフは上の青い線です。それと比較して、(赤い曲線、極のみ)と(金の曲線、のみの極)のグラフです。。fxα10(1x)β11

極限では、下の3つの領域は、互いに相対的にどうなりますか?f

表記の問題として、と記述しますは、と間ののグラフの下の領域です。、、およびの相対サイズについて質問しています。

F(x)=0xf(t)dt=0xtα1(1t)β1dt
f0xF(ϵ)F(1ϵ)F(ϵ)F(1)F(1ϵ)

レッツ推定これらの領域常に仮定一つずつ、と及び。これらの仮定の下で0<α<10<β<1, 0<x<1,0<ϵ<1/2

xα1>1;(1x)β1>1,

xxα1(赤)は減少関数であり、(金)は増加関数です。 x,x(1x)β1

  1. 左側では、青と赤の曲線が近づいているように見えます。実際、場合、前述の不等式は境界間の各積分と簡単であり、圧搾二近い境界間0<x<ϵ

    xα1<xα1(1x)β1<xα1(1ϵ)β1.
    0ϵF(ϵ)
    (1)ϵαα<F(ϵ)<(1ϵ)β1ϵαα.
  2. 同じ分析が右側にも適用され、同様の結果が得られます。

  3. は凹形であるため、中央の区間では、端点で極値が得られます。その結果、面積はこれらの点がまたがる台形の面積より小さくなります:f[ϵ,1ϵ]

    (2)F1εFε<12fε+f1ε1εε=12ε2εα11εβ1+1εα1εβ1

これは厄介を取得する恐れが、のは、一時的に修正してみましょう比率に何が起こるかを検討としてとアプローチ。式とでは、と両方が近づきます。したがって、制限内で重要な唯一の用語はεF1εFεFεαβ0121εα11εβ11ε0=1

(3)F1εFεFεεα1+εβ1/2εα/α=α2ε+α2εαβαε

なぜなら。その結果、、最終的に中央の領域は左側の領域と比較して重要ではなくなります。αβ0α0

同じ議論は、最終的に中央の領域がに正しい領域を近くなり、これも重要ではなくなることを示しています。これはβ/ε

が何であろうと、と両方を十分に小さくすると、基本的に下のすべての領域は左の間隔と正しい間隔。0<ε<1/2 αβf0ε1ε1

残りは簡単です:平均は右極の近くの領域に非常に近くなります(証明:左と中央の区間の積分でをに置き換えて右の間隔で、次いで置換することにより、それを過大評価することによって左側に、中央に、およびでどちらの式も、非常に近似しています。)ただし、、相対面積はおよそバツfバツ0fバツ1εfバツバツfバツεfバツ1εfバツfバツF1F1ε

F1F1εFεε/βε/α=αβ

平均を一定に保つことにより、この比率は一定に保たれ、もう1つの観測を追加できます。

私たちは聞かせている場合してような方法でその制限定数近づい、の領域に右の領域の後、最終的に比率を左も任意に近くなります。α0β0α/βλλ

次に、ゼロに縮小する考えます。その結果、限界分布が存在し、すべての確率が値と集中している必要がありますこれはベルヌーイ分布のクラスです。 どちらをピン留めするか:平均があるベルヌーイ分布は確率を、確率をに割り当てるため、比率は制限比率でなければなりませんε01ppp11p0p/1pλ

質問の用語では、

λ=α/1μμα=μ1μ=p1p

主張どおり。

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