Haftaya göre hacim için çıkarım
10.000 yinelemelik RJAGS benzetim çıktısı rail_sim_1 çalışma alanında; ayrıca Markov zinciri çıktısının bir veri çerçevesi de mevcut:
> head(rail_chains_1, 2)
a b.1. b.2. s
1 420.6966 0 -54.30783 118.2328
2 399.5823 0 -52.02570 119.9499
Bu zincirler, a için (hafta sonu günlerindeki tipik patika hacmi) ve b.2. için (tipik hafta içi hacmi ile hafta sonu hacmi arasındaki karşıtlık) 10.000 benzersiz değer kümesi sağlar. Örneğin, ilk parametre kümesi, hafta sonu günlerinde tipik olarak 420.6966 yolcu ve hafta içi günlerde ise 54.30783 yolcu daha az olduğunu gösterir. Dolayısıyla hafta içi günlerinde tipik olarak 420.6966 - 54.30783 = 366.3888 yolcu vardır. Bu benzetim verilerini, hafta içi patika hacmi hakkında çıkarımlar yapmak için kullanacaksın.
Bu egzersiz, kursun bir parçasıdır
RJAGS ile Bayesçi Modelleme
Egzersiz talimatları
- Tipik hafta içi patika hacmi için 10.000 değerlik bir zincir oluşturmak üzere
aveb.2.zincir değerlerini birleştir. Bunurail_chains_1içindeweekday_meanolarak sakla. ggplot()kullanarakweekday_meanzincir değerlerinin bir yoğunluk grafiğini oluştur.- Tipik hafta içi patika hacmi için %95 güvenilir aralık oluştur.
Uygulamalı etkileşimli egzersiz
Bu egzersizi bu örnek kodu tamamlayarak deneyin.
# Construct a chain of values for the typical weekday volume
rail_chains_1 <- rail_chains_1 %>%
mutate(weekday_mean = ___)
# Construct a density plot of the weekday chain
ggplot(___, aes(x = ___)) +
geom_density()
# 95% credible interval for typical weekday volume