10

JAGS でロジスティック回帰モデルを当てはめようとしていますが、バイナリ変数ではなく (# 成功 y、# 試行 n) の形式のデータがあります。R では、"weights" 引数を指定して glm(y/n ~ ) を使用することで、モデルをこれらのようなデータに適合させることができますが、これを JAGS に適合させる方法がわかりません。

これは、私が尋ねようとしていることに対処するための簡単な例です。rjags パッケージを使用していることに注意してください。助けてくれてありがとう!

y <- rbinom(10, 500, 0.2)
n <- sample(500:600, 10)
p <- y/n
x <- sample(0:100, 10) # some covariate

data <- data.frame(y, n, p, x)

model <- "model{
# Specify likelihood
for(i in 1:10){
    y[i] ~ dbin(p[i], n[i])
    logit(p[i]) <- b0 + b1*x
}

# Specify priors
b0 ~ dnorm(0, 0.0001)
b1 ~ dnorm(0, 0.0001)
}"
4

1 に答える 1