「データ解析のための統計モデリング入門」第11章メモ - もうカツ丼はいいよなp.254の図11.7をstanを使って再現したかった。 stanコードはこんな感じ。 data { int<lower=1> N; // 欠測値を含むデータ数 int<lower=1> N_obs; // 欠測値以外のデータの数 int<lower=1, upper=N> Idx_obs[N_obs]; // 欠測値以外のデータのインデックス int<lower=0> Y_obs[N_obs]; // 欠測値を除いたデータベクトル } parameters { real r[N]; real<lower=0> s_r; } model { for (j in 2:N) r[j] ~ normal(r[j-1], s_r); // 1次のCAR model for (j in 1:N_obs) Y_obs[j] ~ poisson_log(r[Idx_obs[j]]); } generat
