私はこのコードをwinBUGSで書いています:
n <- 100
x1 <- rbinom(n,1,.7)
x2 <- rbinom(n,1,.5)
sum(x1)
sum(x2)
model{
x1 ~ dbin(p1, n) x2 ~ dbin(p2, n) p1 ~ dbeta(a1, b1) p2 ~ dbeta(a2,b2)
diff <- p1 - p2 p.value <- step(diff)
} list(n = 100, x1 = 70, x2 = 54, a1 = 1, b1 = 1, a2 = 1, b2 = 1)
R/JAGS でこれを実行する方法に問題があります。実際のところ、このコードが何をしようとしているのか完全にはわかりません (事後数を計算すると思いますか?)。これまでに winBUGS を使用したことがなく、R も初めてです。これは私の最初のベイズ クラスでもあり、コードが導入されるとかなり迷ってしまいます。
さらに、比率の差の事後平均と標準偏差をどのように計算すればよいでしょうか? p1またはよりも大きい事後確率を見つける方法と、p2それが有意であるかどうか?