Симуляция RJAGS для многомерной регрессии
Рассмотрите следующую байесовскую модель объёма \(Y\)i в зависимости от статуса рабочего дня \(X\)i и температуры \(Z\)i:
- правдоподобие: \(Y\)i \(\sim N(m\)i, \(s^2)\), где \(m\)i \(= a + b X\)i \(+ c Z\)i.
- априорные распределения: \(a \sim N(0, 200^2)\), \(b \sim N(0, 200^2)\), \(c \sim N(0, 20^2)\), \(s \sim Unif(0, 200)\)
Предыдущий анализ взаимосвязи между volume, weekday и hightemp в наборе данных RailTrail дал вам определённое представление об этой зависимости. Теперь вы объедините его с информацией из априорных распределений, чтобы построить апостериорную модель этой взаимосвязи с помощью RJAGS. Набор данных RailTrail доступен в вашем рабочем пространстве.
Это упражнение является частью курса
Байесовское моделирование с RJAGS
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# DEFINE the model
rail_model_2 <-