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

Question

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

基本上我希望它代表:

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

我尝试通过调整响应其他人的帖子而给出的代码来使用 step 函数，该帖子要求类似的内容:Choosing Different Distributions based on if - else condition in WinBugs/JAGS

但这似乎不起作用？

谢谢

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]))
}