JAGSとStanによる負の二項分布のパラメータ推定: Taglibro de H
熊谷さんが負の二項分布のパラメータ推定を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
Models
Bayesian Analysis for Epidemiologists III: Regression Modeling Acknowledgements Much of this material comes from a workshop run by Drs. David Spiegelhalter and Nicky Best in Cambridge, UK, in 2005. Other major contributions come (again) from Lyle Gurrin and Shane Jensen. Jim Albert David Spiegelhalter Nicky Best Andrew Gelman Bendix Carstensen Lyle Guerrin Shane Jensen and Statistical Horizons Sec
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
Negative Binomial Reparameterization