初接触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*theta[j]- bg[i,g]);
ppwd[j,i] <- exp(weightwd*theta[j]- bwd[i,wd]);
pptq[j,i] <- exp(weighttq*theta[j]- btq[i,tq]);
num[j,i] <- (pw + pg*ppg[j,i] + pwd*ppwd[j,i] + ptq*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;
}
}
}