JAGS中的概率ifelse陈述

Question

我希望JAGS模型中的变量被赋值为概率为p [i]的值和另一个概率为1-p [i]的值。 p [i]是模型中的变量。

我尝试了以下但我不确定它是否合适。我担心单个值（0或1）将归属于主题i，并且该值将始终保持不变：对于给定主题i，向量a [i]中的值将始终等于20或30我想要的是一个包含20和30的向量，其中20的概率等于p [i]。

  for (i in 1:n) {
  m[i] ~ dbin(p[i],1) 
  a[i] <- ifelse(m[i]==1, 20, 30)     
  } 

Answer 1

你可以做类似的事情。

for(i in 1:n){
m[i] ~ dbern(p[i])
a[i] <- (m[i] * 20) + ((1 - m[i]) * 30)
}

当Bernoulli路径以概率p取1时，a变为20.当它取0时变为30。