N normal iid ürününün yaklaşık dağılımı? Özel durum μ≈0


12

Verilen iid ve , arayan:N30XnN(μX,σX2)μX0

  1. doğru kapalı form dağıtım yaklaşımı YN=1NXn
  2. aynı ürünün asimtotik ( üstel ?) yaklaşımı

Bu daha genel bir sorunun özel durumudur .μX0


1. ve mı? ( tüm olması iyi olur .) (2) Asimptotik normal bir yaklaşım korkunç olacaktır , çünkü asimptotik uzaktan normal görünmeyecektir. μXσXμX/σX0Y
whuber

Bununla hızlıca oynadım. İlgileniyorsanız, olan rasgele değişkenin ürünü için kesin bir kapalı form çözümü elde etmek mümkündür . Sıfırdan farklı vaka işler çok daha karmaşık hale getirir. nN(0,σ2)μ
wolfies

@whuber (1) bazı farklı ve ile bazı monte carlo yaptıktan sonra , dağılımının ve için oldukça iyi davrandığını ; Şimdi ben güzel bulmak istiyoruz ifadesini için ve nasıl benzer birkaç güzel yaklaşımları vardır. Taylor genişlemesi ile birkaç yaklaşım oluşturdum, ama kötü davranıyorlar. (2) iyi, kesinlikle chi kare ile normalin bir toplamı "görünüyor", bu yüzden eğer "bunu" kanıtlarsa, normale düşürülebilir. μσFN>30|μX|10σXμFσFχ2FF
Andrei Pozolotin

3
Tüm , güzel bir şekilde bir lognormal dağılım ile yaklaşık olacaktır (bir uygulama olarak Barry-Esseen teoremi için gösterir). μX10σXYlog(X)
whuber

@whuber Barry-Esseen doğrudan uygulaması verir güzel gerçekten de, ancak bazı yapı kalktığında,: negatif olmalıdır, bağlı olmalıdır , etc., belki, uygulamanın daha iyi yolları var mı? FN0+1NZμFσFα
Andrei Pozolotin

Yanıtlar:


10

Sıfır ortalama durumunda kesin bir çözüm elde etmek mümkündür (bölüm B).

Sorun

Let göstermektedirler IID değişkenler, ortak pdf her :(X1,,Xn)nN(0,σ2)f(x)

resim açıklamasını buraya girin

, için pdfi=1nXin=2,3,

Çözüm

Bu tür iki Normalin ürününün pdf'si basitçe:

resim açıklamasını buraya girin

... burada MathematicaTransformProduct için mathStatica paketindeki işlevi kullanıyorum . Destek alanı:

resim açıklamasını buraya girin

3, 4, 5 ve 6 Normalin ürünü, aynı fonksiyonun tekrar tekrar uygulanmasıyla (burada dört kez) elde edilir:

resim açıklamasını buraya girin

... burada MeijerGbelirtmektedir Meijer G fonksiyonunu

İndüksiyonla, iid rasgele değişkenlerin ürününün pdf'si :nN(0,σ2)


1(2π)n2σnMeijerG[{{},{}},{{01,,0n},{}},x22nσ2n] for xR

Hızlı Monte Carlo kontrolü

İşte karşılaştırmayı hızlı bir şekilde kontrol edin:

  • yeni elde edilen teorik pdf ( ve ): KIRMIZI DASHED eğrisin=6σ=3
  • ampirik Monte Carlo pdf'ye: dalgalı MAVİ eğri

resim açıklamasını buraya girin

İyi görünüyor! [Mavi dalgalı Monte eğrisi tam kırmızı kesikli eğriyi gizliyor]


Olağanüstü, teşekkürler, Colin. Şimdi kitabınızı neden satın almam gerektiğini görüyorum :-) Ayrıca daha basit görünüyor mu merak ediyor . Wolfram becerilerimin tozunu alma zamanı. log(...MeijerG(...))
Andrei Pozolotin
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.