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]; // 欠測値を除いたデー</lower=0></lower=1,></lower=1></lower=1>…