負の二項分布のベイズ推定 [統計]
負の二項分布のパラメータをJAGSで推定するテスト。
データは
BUGSコード
データは
n <- 100 p <- 0.8 r <- 10 num <- rnbinom(n, r, p)
BUGSコード
var N, p, r, num[N] model { for (i in 1:N) { num[i] ~ dnegbin(p, r); } # priors p ~ dunif(0, 1); r ~ dgamma(1.0E-3, 1.0E-3); }
rは整数でなくてもいいみたい。
コメント 0