悬赏 5 个论坛币 未解决
初接触bugs,遇到麻烦希望大家帮忙解决。
weight和g[],wd[],tq[]都是设定的数据,需要估出theta[j], bg[i,1], bwd[i,1], btq[i,1]
我是想设bg[i,1]服从正态,bg[i,2]不需要估计,设定其为= 0
但compiling时,提示multiple definitions of node bg[1,1],水平太差,请达人指教!
model{
for (j in 1:nperson){
theta[j] ~ dnorm(0, 1);
for (i in 1:nitem){
Y[j, i] ~ dbern(prob[j, i]);
ppg[j,i] <- exp(weightg[i]*theta[j]- bg[i,g[i]]);
ppwd[j,i] <- exp(weightwd[i]*theta[j]- bwd[i,wd[i]]);
pptq[j,i] <- exp(weighttq[i]*theta[j]- btq[i,tq[i]]);
num[j,i] <- (pw[i] + pg[i]*ppg[j,i] + pwd[i]*ppwd[j,i] + ptq[i]*pptq[j,i])
den[j,i] <- (1+weightg[i]*ppg[j,i] + weightwd[i]*ppwd[j,i] + weighttq[i]*pptq[j,i]);
prob[j, i] <- num[j,i]/den[j,i];
for(k in 1:nt){
bg[i,k] ~ dnorm(0, .5);
bwd[i,k] ~ dnorm(0, .5);
btq[i,k] ~ dnorm(0, .5);
}
bg[i,2] <- 0;
bwd[i,2] <- 0;
btq[i,2] <-0;
}
}
}