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
![R: hSDMパッケージをためす: Taglibro de H](https://cdn-ak-scissors.b.st-hatena.com/image/square/bebce60b964ff3b840869ad631117eec4c37395f/height=288;version=1;width=512/https%3A%2F%2Fito-hi.c.blog.ss-blog.jp%2F_images%2Fblog%2F_5c3%2Fito-hi%2FRplot01-88fa6.png%3Fc%3Da2)