首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >节点W[1]的多个定义

节点W[1]的多个定义
EN

Stack Overflow用户
提问于 2014-06-15 21:24:20
回答 1查看 154关注 0票数 1

我试图在条件和( deltaj,1:5)=0的条件下,对每个"j“估计deltaj,k。但是当我编译代码时,软件输出是“节点W1的多个定义”。有人能帮我吗?

代码语言:javascript
运行
复制
model {
for (j in 1:p){
for (k in 1:5){
Z[j, k]<- sum(delta[j,1:k])
}
for (i in 1:n){
Y[i , j] ~ dcat ( prob [i , j , 1: 5])
}}
for (i in 1:n){
theta [i] ~ dnorm (0.0 , 1.0)
}
for (i in 1:n){
for (j in 1:p){
for (k in 1:5){
eta[i , j , k] <- alpha [j] * (k*theta [i] - k*beta [j]+Z[j, k])
psum [i , j , k] <- sum(eta[i , j , 1: k])
exp.psum[i , j , k]<- exp( psum [i , j , k])
prob [i , j , k] <- exp.psum[i , j , k] / sum(exp.psum [i , j , 1:5])
}  } }
for (j in 1:p){
W[j] <- sum(delta [j, 1:5])
W[j]<- 0
alpha [j] ~ dlnorm (0.83 , pr.alpha)
beta [j] ~ dnorm (-1.73 , pr.beta ) 
delta[j,1] <- 0.0
for (k in 2:5){
delta [j , k] ~ dnorm (0.02 , pr.delta )
} }
pr.alpha <- pow(1.2 , -2)
pr.beta <- pow(0.7, -2)
pr.delta <- pow(1.3, -2)
}

谢谢

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2014-06-16 09:14:24

BUGS不允许您覆盖确定性节点.在上一个W[j] <-循环中有两次for

我想有很多方法来编写符合你条件的代码。例如,您可以对delta使用不同的发行版,或者将delta[1]设置为在模拟delta[2]delta[5]之后将所有delta之和为0所需的余数。

票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/24234247

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档