1

わかりました、コードを書き直します。問題は、コンパイルすると、「ノードaの複数定義」というエラーが発生することです。私のコードで何が間違っているのか誰にもわかりません。モデルが多くの定数を持たないように、変数 a、b、および c を作成します。 

  model{



  for(i in 1:n){


            a <- (k[1] + step(s1[i]-.9)*k[2] + step(s1[i]*.5-.9)*k[3]) 
            b <- (r[1] + step(s2[i]-.9)*r[2] + step(s2[i]*.5-.9)*r[3])
            c <- (s[1] + step(s3[i]-.9)*s[2] + step(s3[i]*.5-.9)*s[3])


    dummy[i] <- 0

      dummy[i] ~ dloglik(logLike[i])

      # This is the log transformation of the 3-variate poisson 

      logLike[i] <- -theta12[i] + a*log(theta12[i]) - logfact(a) -theta13[i] - b*log(theta13[i]) - logfact(b)-theta23[i] - c*log(theta23[i]) - logfact(c)-theta1[i] + (y1[i]-a-b)*log(theta1[i]) + logfact(y1[i]-a-b)-theta2[i] + (y2[i]-a-c)*log(theta2[i]) + logfact(y2[i]-a-c)-theta3[i] + (y3[i]-b-c)*log(theta3[i]) + logfact(y3[i]-b-c)





            log(theta1[i]) <- alpha1 + beta1*log(L[i])
            log(theta2[i]) <- alpha2 + beta2*log(L[i])
            log(theta3[i]) <- alpha3 + beta3*log(L[i])
            log(theta12[i]) <- alpha4 + beta4*log(L[i])
            log(theta23[i]) <- alpha5 + beta5*log(L[i])
            log(theta13[i]) <- alpha6 + beta6*log(L[i])
    }

            alpha1 ~ dgamma(.001, .001)
            alpha2 ~ dgamma(.001, .001)
            alpha3 ~ dgamma(.001, .001)
            alpha4 ~ dgamma(.001, .001)
            alpha5 ~ dgamma(.001, .001)
            alpha6 ~ dgamma(.001, .001)
            beta1 ~ dgamma(.001, .001)
            beta2 ~ dgamma(.001, .001)
            beta3 ~ dgamma(.001, .001)
            beta4 ~ dgamma(.001, .001)
            beta5 ~ dgamma(.001, .001)
            beta6 ~ dgamma(.001, .001)




    }


    data
    list(y1 =c(0,1,2,3,2),
    y2 = c(1,2,3,0,2),
    y3 = c(1,3,2,3,2),
    k =   c(0,1,2,3,4),
    s =   c(0,1,2,3,4), 
    r =    c(0,1,2,3,4), 
    L =   c(1,2,3,4,5), 
    s1= c(0,1,2,0,2), 
    s2= c(0,0,0,0,0),
    s3=  c(1,1,0,0,0),n=5)


    inits
    list(beta1=0,beta2=0,beta3=0,beta4=0,beta5=0,beta6=0,alpha1=0,alpha2=0,alpha3=0,
    alpha4=0,alpha5=0,alpha6=0)
4

1 に答える 1

8

for() ループ内で a、b、c を定義したため、このエラーが発生します。つまり、定数ノードを複数回 (n 回) 定義しました。a、b、c を a[i]、b[i]、c[i] に置き換える必要があります。BUGS は、プログラミング言語ではなく、モデル記述言語と考えてください。モデル内のすべての量は、他の量の固定関数またはランダム関数として 1 つの定義を持つ必要があります。

于 2013-01-28T09:28:46.230 に答える