Визначте, скомпілюйте й змоделюйте Normal-Normal

Спостерігаючи зміну часу реакції \(Y\)_i для кожного з 18 учасників \(i\) у дослідженні сну, ви можете оновити свою апостеріорну модель впливу нестачі сну на час реакції. Для цього потрібно поєднати інформацію з правдоподібності та апріорних моделей:

правдоподібність: \(Y\)_i \(\sim N(m, s^2)\)
апріорі: \(m \sim N(50, 25^2)\) та \(s \sim Unif(0, 200)\)

У цій серії вправ ви будете визначати, компілювати та імітувати свою баєсівську апостеріорну модель. Спостережувані дані sleep_study є у вашому робочому середовищі.

Ця вправа є частиною курсу

Баєсівське моделювання з RJAGS

Переглянути курс

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

# DEFINE the model    
___ <- "model{
    # Likelihood model for Y[i]
    for(i in 1:___) {
        Y[i] ~ ___
    }

    # Prior models for m and s
    m ~ ___
    s ~ ___
}"

Редагувати та запускати код