Kesalahan dalam perkiraan normal untuk distribusi jumlah yang seragam

20

Salah satu metode naif untuk mendekati distribusi normal adalah dengan menambahkan bersama-sama mungkin variabel acak IID yang terdistribusi secara seragam pada , kemudian masuk kembali dan skala, bergantung pada Teorema Batas Pusat. ( Catatan : Ada metode yang lebih akurat seperti transformasi Box-Muller .) Jumlah variabel acak IID dikenal sebagai distribusi jumlah seragam atau distribusi Irwin-Hall .[ 0 , 1 ]100[0,1]U(0,1)

Berapa besar kesalahan dalam memperkirakan distribusi jumlah seragam dengan distribusi normal?

Kapan pun jenis pertanyaan ini muncul untuk mendekati jumlah variabel acak IID, orang-orang (termasuk saya) memunculkan Teorema Berry – Esseen , yang merupakan versi efektif dari Teorema Batas Pusat mengingat bahwa momen ketiga ada:

|Fn(x)-Φ(x)|Cρσ3n

di mana adalah fungsi distribusi kumulatif untuk jumlah yang dihitung kembali dari variabel acak IID, adalah momen sentral ketiga mutlak, adalah standar deviasi, dan adalah konstanta absolut yang dapat dianggap 1 atau bahkan 1/2 . n ρ E | ( X - E X ) 3 | σ C 1 1 / 2FnnρE|(X-EX)3|σC11/2

Ini tidak memuaskan. Tampak bagi saya bahwa perkiraan Berry-Esseen paling dekat dengan tajam pada distribusi binomial yang terpisah, dengan kesalahan terbesar pada 0 untuk distribusi binomial simetris. Kesalahan terbesar terjadi pada lompatan terbesar. Namun, distribusi jumlah seragam tidak memiliki lompatan.

Tes numerik menunjukkan bahwa kesalahan menyusut lebih cepat daripada c/n .

Menggunakan C=1/2 , estimasi Berry – Esseen adalah

|Fn(x)Φ(x)|121321123n0.650n

yang untuk n=10,20,40 sekitar 0.205 , 0,145 , dan 0,103 , masing-masing. Perbedaan maksimum sebenarnya untuk n=10,20,40 tampaknya masing-masing sekitar 0,00281 , 0,00139 , dan 0,000692 , yang jauh lebih kecil dan tampaknya jatuh sebagai c/n bukannya c/n .

Douglas Zare
sumber
7
Jika Anda memperluas distribusi jumlah dalam ekspansi Edgeworth , Anda menemukan bahwa seragam dalam sebagai (karena distribusi seragam simetris), jadi terdengar benar. Karena istilah , itu tidak memberi Anda batasan ...x n c / n o ( n - 1 )Fn(x)=Φ(x)+n-1g(x)+Hai(n-1)xnc/nHai(n-1)
MånsT
1
Terima kasih, sepertinya ini menjelaskan pola untuk banyak distribusi lainnya juga. c/n
Douglas Zare

Jawaban:

17

Misalkan menjadi iid variabel acak dan pertimbangkan jumlah dinormalisasi dan norm mana adalah distribusi .U ( - b , b ) S n = U1,U2,...U(-b,b)sup δ n = sup x R | F n ( x ) - Φ ( x ) |

Sn=3saya=1nUsayabn,
supF n S n
δn=supxR|Fn(x)-Φ(x)|,
FnSn

Lemma 1 ( Uspensky ): berikut pada memegang . δ n < 1δn

δn<17.5πn+1π(2π)n+12π3nexp(-π2n/24).

Bukti . Lihat JV Uspensky (1937), Pengantar probabilitas matematika , New York: McGraw-Hill, hal. 305.

Ini kemudian diperbaiki oleh R. Sherman sebagai berikut.

Lemma 2 ( Sherman ): Peningkatan berikut ini berlaku untuk Uspensky.

δn<17.5πn-(π180+17.5πn)e-π2n/24+1(n+1)π(2π)n+12π3ne-π2n/24.

Bukti : Lihat R. Sherman, Kesalahan perkiraan normal dengan jumlah variabel acak N , Biometrika , vol. 58, tidak. 2, 396–398.

Buktinya adalah aplikasi yang cukup mudah dari ketidaksetaraan segitiga dan batas klasik pada ekor distribusi normal dan pada diterapkan pada fungsi karakteristik masing-masing dari dua distribusi.(dosax)/x

kardinal
sumber
2
N=n
@Prastrastator: Tangkapan yang bagus.
kardinal
1
2