[精彩WinBUGS答问]Error Message in WinBUGS:Array index is less than one

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

for (i in 1:160) {
    situps ~ dnorm(mu,prec[gender+1]);
    mu <- b0+b1[gender]+b2[health]+b1b2[health,gender];
}
b0 ~ dnorm(0,.001);
for (i in 1:4) {
    b2 ~ dnorm(0,.001);
}
for (j in 1:2) {
    b1[j] ~ dnorm(0,.001);
    prec[j] <- pow(sigma[j],-2);
    sigma[j] ~ dunif(0,10);
}
for (i in 1:4) {
    for (j in 1:2) {
        b1b2[i,j] ~ dnorm(0,.001);
    }
}
I'm not sure why this gives me the error "array index is less than one" because I clearly have all of my arrays starting at 1. Any ideas?

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：winbugs Message WINBUG array Error Error

How is gender coded? Your use of prec[gender[i]+1] looks problematic. Try something like

situps[i]~dnorm(mu[i],tau[i])
tau[i] <- prec[gender[i]] # or some similar function
mu[i]<-b0+b1[gender[i]]+b2[health[i]]+b1b2[health[i],gender[i]]
NB you don't need the semi-colons at the end of each line

Indexing with gender[i] + 1 doesn't seem as problematic as indexing with gender[i] and health[i] in this case

你好！
请问你解答出来了吗？我在编译的过程中也遇到这个问题。array index is less than one
麻烦再加QQ327871324讨论一下呢。谢谢！