对于先前的约束度量,我试图在[-1,1]之间延伸β分布,“[a] s描述为Barnard,McCulloch& Meng(2000)”(根据{{3 }})。
具体来说,我正在尝试实施这个建议:
rho_half_with ~ dbeta(1, 1)
# shifting and streching rho_half_with from [0,1] to [-1,1]
rho ~ 2 * rho_half_with - 1
然而,我总是得到
syntax error on line (...) near "2"
this tutorial或BUGS手册中没有条目我发现处理分布操作(作为随机关系分配的来源)。是否确实可以将基本算术运算应用于BUGS / JAGS随机关系(遵循~
运算符),如果是,如何?
答案 0 :(得分:2)
您发布的代码的问题是您在非随机关系中使用~
,而JAGS希望您使用<-
。以下应该有效:
rho_half_with ~ dbeta(1, 1)
# shifting and streching rho_half_with from [0,1] to [-1,1]
rho <- 2 * rho_half_with - 1
关于您在评论中提到的错误消息,因为您尝试启动非随机变量(rho
)。删除该初始化或切换到初始化rho_half_with
以解决该问题。