entegre etmek istediğim bir fonksiyonum olduğunu varsayın Tabii ki 'nin uç noktalarda sıfıra gittiğini varsayarsak , patlama olmaz, güzel işlev. I ile ilgilenmek oldum bir yolu örnekleri bir listesini oluşturmak için Metropolis- Hastings yapısı kullanmaktır dağılımından orantılı için normalleştirme sabit eksik, N = \ int _ {- \ infty} ^ {\ infty} g (x) dx p (x) diyeceğim ve sonra bu x 'lerde bazı istatistik f (x) hesaplayacağım :
Yana , içinde yerini alabilir iptal etmek formunun bir ifade ile sonuçlanır entegrali gelen
Bunu örnek fonksiyonu için R'de test etmeyi denedim . Bu örnekte Metropolis-Hastings'i örnek üretmek için kullanmıyorum ama gerçek olasılıkları örnek üretmek için kullanıyorum (sadece test etmek için). Aradığım sonuçları tam olarak alamıyorum. Temel olarak hesaplayacağım şeyin tam ifadesi:
Bu benim . Yaklaşıyor ama kesinlikle beklenen şekilde birleşmiyor, yanlış bir şey mi yapıyorum? 1rnorm
1 / √
ys = rnorm(1000000, 0, 1/sqrt(2))
r = max(ys) - min(ys)
sum(sapply(ys, function(x) 1/( r * exp(-x^2))))/length(ys)
## evaluates to 0.6019741. 1/sqrt(pi) = 0.5641896
CliffAB için düzenle
Aralığı kullanmamın nedeni, noktalarımın bulunduğu bölgede sıfır olmayan ancak aralığında ile bütünleşen bir işlevi kolayca tanımlamaktır . Fonksiyonun tam özellikleri: Bu eşit yoğunluk olarak kullanmak zorunda kalmadım. bütünleştirilmiş başka bir yoğunluk kullanabilirdim , örneğin olasılık yoğunluğu Ancak bu, bireysel örneklerin toplanmasını önemsiz hale getirecektir. [ - ∞ , ∞ ] U ( x ) = { 1U(x)1P(x)=1
Bu tekniği entegre olan diğer dağıtımlar için deneyebilirim . Ancak yine de neden düzgün bir dağıtım için işe yaramadığını bilmek istiyorum.