Определение, компиляция и симуляция регрессионной модели
Наблюдая зависимость между весом \(Y\)i и ростом \(X\)i для 507 испытуемых \(i\) из набора данных bdims, вы можете обновить апостериорную модель этой зависимости. Чтобы построить апостериорное распределение, необходимо объединить информацию из функции правдоподобия и априорных распределений:
- правдоподобие: \(Y\)i \(\sim N(m\)i, \(s^2)\), где \(m\)i \(= a + b X\)i
- априорные распределения: \(a \sim N(0, 200^2)\), \(b \sim N(1, 0.5^2)\) и \(s \sim Unif(0, 20)\)
В этой серии упражнений вы определите, скомпилируете и симулируете апостериорное распределение байесовской регрессии. Данные bdims уже загружены в вашем рабочем пространстве.
Это упражнение является частью курса
Байесовское моделирование с RJAGS
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# DEFINE the model
weight_model <- ___