R Package for Hierarchical Bayesian species distribution models now available | Adam M. Wilsonという記事で紹介されていたhSDMパッケージをためしてみた。 例題として、久保本11章の空間構造のある階層ベイズモデルをやってみる。 library(hSDM) load(url("http://hosho.ees.hokudai.ac.jp/~kubo/stat/iwanamibook/fig/spatial/Y.RData")) n.site <- length(Y) adj <- c(2, sapply(2:(n.site - 1), function(i) c(i - 1, i + 1)), n.site - 1) n.adj <- c(1, rep(2, n.site - 2), 1) data