Моделювання в RJAGS для регресії Пуассона
У попередньому відео ми побудували модель регресії Пуассона для обсягу \(Y\)i за статусом дня тижня \(X\)i та температурою \(Z\)i:
- правдоподібність: \(Y\)i \(\sim Pois(l\)i), де \(log(l\)i\() = a + b X\)i \(+ c Z\)i
- пріори: \(a \sim N(0, 200^2)\), \(b \sim N(0, 2^2)\) і \(c \sim N(0, 2^2)\)
Поєднавши ваші висновки з наявних даних RailTrail і зазначених пріорів, ви визначите, скомпілюєте та змоделюєте апостеріорну модель цієї залежності за допомогою RJAGS. Щоб додати виклику у цьому останньому моделюванні RJAGS у курсі, ви отримаєте менше допоміжного коду, ніж зазвичай!
Дані RailTrail є у вашому робочому середовищі.
Ця вправа є частиною курсу
Баєсівське моделювання з RJAGS
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# DEFINE the model
poisson_model <-