小编Cam*_*des的帖子

我需要帮助才能在我的WINBUGS代码中找到错误(v.1.4.3).

在"模型规范"步骤中,模型看起来在语法上是正确的.但是,在我尝试加载数据时,我收到此错误:

数组索引大于phi3的数组上限

有人可以帮帮我吗？我的模型如下:

model {

        for(w in 1: W){
        m[w] <- n[w]-y1[w]
        h[w] <- n[w]-y1[w]-y2[w]
        y1[w] ~ dbin(delta[w],n[w])
        y2[w] ~ dbin(theta[w],m[w])
        y3[w] ~ dbin(eta[w],h[w])
        y4[w] <- n[w]-y1[w]-y2[w]-y3[w]
        logit(delta[w]) <- mu1+theta1[a[w]]+phi1[p[w]]+psi1[c[w]]
        logit(theta[w]) <- mu2+theta2[a[w]]+phi2[p[w]]+psi2[c[w]]
        logit(eta[w])   <- mu3+theta3[a[w]]+phi3[p[w]]+psi3[c[w]]
    }

    ## Autoregressive prior model for p effects

    phi1mean[1] <- 0.0
    phi1prec[1] <- tauphi1*1.0E-6
    phi1mean[2] <- 0.0
    phi1prec[2] <- tauphi1*1.0E-6

    phi2mean[1] <- 0.0
    phi2prec[1] <- tauphi2*1.0E-6
    phi2mean[2] <- 0.0
    phi2prec[2] <- tauphi2*1.0E-6

    phi3mean[1] <- 0.0
    phi3prec[1] <- tauphi3*1.0E-6
    phi3mean[2] <- 0.0
    phi3prec[2] <- …