前回の記事についてbaibai先生から、dlmと尤度がガウスの1次元CARで結果を比較してみるとよいかもといわれましたので、やってみました。 前回のJAGSのモデルのうち、観測モデルの部分を正規分布にしたものと、dlmの結果とを比較してみます。dlmのシステムモデルはローカルレベルモデルです。 Rコード Y <- c( 0, 3, 2, 5, 6, 16, 8, 14, 11, 10, 17, 19, 14, 19, 19, 18, 15, 13, 13, 9, 11, 15, 18, 12, 11, 17, 14, 16, 15, 9, 6, 15, 10, 11, 14, 7, 14, 14, 13, 17, 8, 7, 10, 4, 5, 5, 7, 4, 3, 1) N.site <- length(Y) ## dlm library(dlm) build1 <- function