2015-03-02 2 views
0

У меня проблема в коде ниже, проблема возникает при компиляции процесса «множественное определение узла xi1 [1,1]», кто-нибудь поможет мне решить эту проблему, пожалуйста. много спасибо заранеемножественное определение узла xi1 [1,1]

model { 
for(i in 1:N){ 
#measurement equation model 
for(j in 1:P){y[i,j]~dnorm(mu[i,j],psi[j])I(thd[j,z[i,j]],thd[j,z[i,j]+1])} 
xi[i,1]<-mu[i,1]+lam[1]*mu[i,2]+lam[2]*mu[i,3]+lam[3]*mu[i,4]+lam[4]*mu[i,5]+lam[5]*mu[i,6] 
xi[i,2]<-mu[i,1]+lam[1]*mu[i,2]+lam[2]*mu[i,3]+lam[3]*mu[i,4]+lam[4]*mu[i,5]+lam[5]*mu[i,6]+xi[i,1] 
#structural equation model 
xi[i,1:2]~dmnorm(zero2[1:2],phi[1:2,1:2]) 
eta[i]~dnorm(nu[i],psd1) 
nu[i]<-gam[1]*xi[i,1] 
mu[i,1:6]~dnorm(0.8,4.0) 
}# end of i 
for(i in 1:2){zero2[i]<-0.0} 
#priors on loadings and coefficients 
for(i in 1:5){lam[i]~dnorm(0.8,4.0)} 
for(i in 1:1){gam[i]~dnorm(0.6,4.0)} 
#priors on precisions 
for(j in 1:P){ 
psi[j]~dgamma(10,8) 
sgm[j]<-1/psi[j]} 
psd~dgamma(10,8) 
sgd<-1/psd 
phi[1:2,1:2]~dwish(R[1:2,1:2], 30) 
phx[1:2,1:2]<-inverse(phi[1:2,1:2]) 
} #end of model 

ответ

0

у вас есть эти две строки:

xi[i,1]<-mu[i,1]+lam[1]*mu[i,2]+lam[2]*mu[i,3]+lam[3]*mu[i,4]+lam[4]*mu[i,5]+lam[5]*mu[i,6] 
... 
xi[i,1:2]~dmnorm(zero2[1:2],phi[1:2,1:2]) 

Так как оба в цикле я, вы определили хх [1,1] дважды: один раз в детерминированной количестве и один раз в качестве случайной переменной (данные?).

+0

Большое вам спасибо за ваш ответ, что мне теперь делать, если я удалю второе определение xi [i, 1: 2], ошибка будет «ожидаемым многомерным узлом», поэтому я не могу удалить ни одного из них. любая помощь будет оценена по достоинству. – thanoon

+0

Ошибка теперь представляет собой «ожидаемый многомерный узел», не удаляя ни одного из них. – thanoon

Смежные вопросы