用于在 JAGS 中定义分布的 if/else 语句

Question

在 JAGS 中，我想为参数 w[i] 定义泊松分布，如果另一个参数 e[i] 大于 0，它也会被截断（大于或等于 2）。

基本上我希望它代表：

w[i] ~ ifelse( e[i] > 0, dpois(mu) T(2,) , dpois(mu) )

我尝试通过调整响应其他人的帖子而给出的代码来使用 step 函数，该帖子要求类似的内容： 根据 WinBugs/JAGS 中的 if - else 条件选择不同的发行版

但这似乎不起作用？

谢谢

Answer 1

也许是这样的？

pois1 ~ dpois(mu) T(2,)
pois2 ~ dpois(mu)
for(i in 1:N){
indicator1[i] <- ifelse(e[i] > 0, 1, 0)
indicator2[i] <- ifelse(e[i] <= 0, 1, 0)
w[i] <- (pois1 * indicator1[i]) + (pois2 * indicator2[i])
}

当e[i]大于 1 时w[i]从 中取值pois1。如果它不是w[i]从pois2.

编辑：或者，您可以只定义一个指标变量并像这样进行操作。

pois1 ~ dpois(mu) T(2,)
pois2 ~ dpois(mu)
for(i in 1:N){
indicator[i] <- ifelse(e[i] > 0, 1, 0)
w[i] <- (pois1 * indicator[i]) + (pois2 * (1 - indicator[i]))
}