SSブログ

負の二項分布のベイズ推定 [統計]

負の二項分布のパラメータをJAGSで推定するテスト。

データは
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は整数でなくてもいいみたい。



negbin.png


nice!(0)  コメント(0)  トラックバック(1) 
共通テーマ:学問

nice! 0

コメント 0

コメントを書く

お名前:
URL:
コメント:
画像認証:
下の画像に表示されている文字を入力してください。

Facebook コメント

トラックバック 1