Art değer olasılıkları
RJAGS çıktısını kullanarak art trendini ve belirsizliği \(b\) keşfettin ve nicelleştirdin. RJAGS çıktısını belirli hipotezleri değerlendirmek için de kullanabilirsin. Örneğin: Ortalama olarak, boydaki her 1 cm artış için ağırlığın 1.1 kg'dan fazla arttığına dair art olasılık nedir? Yani, \(b > 1.1\) olmasının art olasılığı nedir?
Bu olasılığı, 1.1'i aşan \(b\) Markov zinciri değerlerinin oranıyla yaklaşık olarak hesaplayacaksın. 100.000 iterasyonluk Markov zinciri çıktısını içeren weight_chains veri çerçevesi çalışma alanında.
Bu egzersiz, kursun bir parçasıdır
RJAGS ile Bayesçi Modelleme
Egzersiz talimatları
- \(b\) Markov zinciri değerlerinin bir yoğunluk grafiğini oluştur ve 1.1 konumuna dikey bir çizgi bindirmek için
geom_vline()kullan. - 1.1'i aşan \(b\) Markov zinciri değerlerinin sayısını özetlemek için
table()kullan. - 1.1'i aşan \(b\) Markov zinciri değerlerinin oranını hesaplamak için
mean()kullan.
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# Mark 1.1 on a posterior density plot for b
ggplot(___, aes(x = ___)) +
geom_density() +
geom_vline(xintercept = ___, color = "red")
# Summarize the number of b chain values that exceed 1.1
# Calculate the proportion of b chain values that exceed 1.1