Question

Difficulty: HardBayes Tahmin Yöntemi ve Sonsal (Posterior) Dağılım

Bir tarım il müdürlüğünde görevli ziraat mühendisi, belirli bir tarım arazisinde metrekareye düşen zararlı böcek sayısını (XX) modellemektedir. Zararlı böcek sayısının θ\theta parametreli Poisson dağılımına uyduğu varsayılmaktadır.

Geçmiş yıllara ait veriler incelendiğinde, θ\theta parametresinin önsel (prior) dağılımının α=4\alpha=4 ve β=2\beta=2 parametreli Gamma dağılımı olduğu belirlenmiştir. (Gamma dağılımının olasılık yoğunluk fonksiyonu f(θ)=βαΓ(α)θα1eβθ,θ>0f(\theta) = \frac{\beta^\alpha}{\Gamma(\alpha)} \theta^{\alpha-1} e^{-\beta\theta}, \theta > 0 olarak verilmektedir.)

Mühendis, araziden rastgele seçilen n=5n=5 farklı metrekarelik alanda yaptığı sayımlarda sırasıyla 2,4,3,22, 4, 3, 2 ve 55 adet zararlı böcek tespit etmiştir.

Tahmin işleminde kullanılacak kayıp fonksiyonu L(θ,δ)=θ(θδ)2L(\theta, \delta) = \theta(\theta - \delta)^2 olarak belirlendiğine göre, θ\theta parametresi için Bayes tahmin edicisinin aldığı değer aşağıdakilerden hangisidir?

  1. 33Answer
  2. B
    207\frac{20}{7}
  3. C
    197\frac{19}{7}
  4. D
    165\frac{16}{5}
  5. E
    52\frac{5}{2}

Answer

Bayes tahmin edicisi 33 olarak bulunur.
Bayes tahmin edicisi, verilen kayıp fonksiyonunun sonsal (posterior) beklenen değerini minimize eden tahmin edicidir. Sonsal dağılım, önsel (prior) Gamma(4, 2) ve olabilirlik Poisson(16, 5) fonksiyonlarının çarpımı ile Gamma(20, 7) olarak bulunur. Soruda verilen ağırlıklı karesel kayıp fonksiyonu L(θ,δ)=θ(θδ)2L(\theta, \delta) = \theta(\theta - \delta)^2 için beklenen risk δ\delta'ya göre minimize edildiğinde Bayes tahmin edicisi genel formülü δ=E[θ2x]E[θx]\delta = \frac{E[\theta^2|x]}{E[\theta|x]} olarak elde edilir. Gamma(20, 7) dağılımının özellikleri (E[θ]=αβE[\theta]=\frac{\alpha}{\beta} ve E[θ2]=α(α+1)β2E[\theta^2]=\frac{\alpha(\alpha+1)}{\beta^2}) kullanıldığında sonuç 217=3\frac{21}{7} = 3 olarak bulunur.

Step-by-Step Solution

1
Olabilirlik fonksiyonunun ve önsel dağılımın belirlenmesi.
Olabilirlik: L(θx)θxienθ=θ16e5θL(\theta|x) \propto \theta^{\sum x_i} e^{-n\theta} = \theta^{16} e^{-5\theta}. Önsel: π(θ)θα1eβθ=θ3e2θ\pi(\theta) \propto \theta^{\alpha-1} e^{-\beta\theta} = \theta^{3} e^{-2\theta}.
Sonsal dağılımı bulmak için Bayes teoremi gereği olabilirlik fonksiyonu ile önsel dağılım çarpılmalıdır.
2
Sonsal (posterior) dağılımın hesaplanması.
π(θx)L(θx)π(θ)=θ19e7θ\pi(\theta|x) \propto L(\theta|x) \pi(\theta) = \theta^{19} e^{-7\theta}. Bu ifade, αpost=20\alpha_{post}=20 ve βpost=7\beta_{post}=7 parametreli Gamma dağılımının çekirdeğidir.
Sonsal dağılım, Bayes parametre tahmininde kullanılacak temel dağılımdır.
3
Verilen kayıp fonksiyonu altında Bayes tahmin edicisi formülünün türetilmesi.
L(θ,δ)=θ(θδ)2L(\theta, \delta) = \theta(\theta - \delta)^2 için beklenen kayıp E[Lx]=E[θ32δθ2+δ2θx]E[L|x] = E[\theta^3 - 2\delta\theta^2 + \delta^2\theta | x] olur. Bu ifadenin δ\delta'ya göre türevi alınıp sıfıra eşitlendiğinde δ=E[θ2x]E[θx]\delta = \frac{E[\theta^2|x]}{E[\theta|x]} elde edilir.
Bayes tahmin edicisi, sonsal beklenen kaybı (riski) en küçük yapan δ\delta değeridir.
4
Gamma dağılımının momentlerini kullanarak tahmin edicinin sayısal olarak hesaplanması.
Gamma(20,7)Gamma(20, 7) dağılımı için E[θx]=207E[\theta|x] = \frac{20}{7} ve E[θ2x]=20×2172E[\theta^2|x] = \frac{20 \times 21}{7^2}'dir. Formülde yerine konulduğunda δ=20×21/4920/7=217=3\delta = \frac{20 \times 21 / 49}{20 / 7} = \frac{21}{7} = 3 bulunur.
Sonsal dağılımın momentleri türetilen formülde kullanılarak nihai nokta tahmini elde edilir.

Key Concept

Bayes Tahmin Yöntemi ve Ağırlıklı Kayıp Fonksiyonları
Estimated Time:2m 30s
Rate this question