熊谷さんが負の二項分布のパラメータ推定をWinBUGSでやっておられたので、JAGSとStanでもやってみました。 Rのコード ## ## Negative Binomial ## ## JAGS library(rjags) set.seed(1234) n <- 100 mu <- exp(2) r <- 1.5 Y <- rnbinom(n, size = r, mu = mu) params <- c("mu", "r") inits <- lapply(1:3, function(i) list(p = runif(1, 0, 1), r = rpois(1, 20) + 1)) model <- jags.model("negbin.bug", data = list(N = n, Y = Y), inits = inits, n.chains = 3, n.adapt = 5