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
![「データ解析のための統計モデリング入門」第11章メモ - もうカツ丼はいいよな](https://cdn-ak-scissors.b.st-hatena.com/image/square/e1e1ac83e8f9552b84a4b298f4e7043faf7e10d8/height=288;version=1;width=512/https%3A%2F%2Fcdn-ak.f.st-hatena.com%2Fimages%2Ffotolife%2FR%2FRion778%2F20161008%2F20161008220158.png)