Winbugs中的“多重定义节点a”错误

好吧我只是重写我的代码。现在的问题是，当我编译它，我得到一个错误“节点a的多重定义”。有谁知道我的代码有什么问题。我为模型创建了变量a，b和c，但没有多少常量。Winbugs中的“多重定义节点a”错误

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)

来源

2013-01-24 Ansel Stars

我回滚了你的编辑，因为如果你删除你的代码并用咆哮取代它，任何人都不可能帮助你。 – josh3736

@ josh3736，谢谢你。 –

你得到这个错误是因为你已经在for（）循环中定义了a，b，c。换句话说，您已经多次定义了一个常量节点（n次）。你应该用a [i]，b [i]，c [i]代替a，b，c。将BUGS看作模型描述语言，而不是编程语言。模型中的每个数量都必须具有一个单一的定义，作为其他数量的固定或随机函数。

来源

2013-01-28 09:28:46

Winbugs中的“多重定义节点a”错误

回答

相关问题