Примечание. Изначально это сообщение было отправлено на доску обсуждений JAGS на SourceForge .
У меня возникли некоторые проблемы, и я не знаю, почему. Я также должен отметить, что это мой первый проект, использующий сэмплер dnormmix. Код, который, кажется, имеет ошибку:
for( i in 1:N ){
y[i] ~ dnormmix( mu[i,], mtau, pi[i,] )
mu[i,1] <- 0
mu[i,2] <- h[i]*alpha + theta[w[i],x1[i]] - theta[w[i],x2[i]]
pi[i,1] <- 1-pi[i,2]
pi[i,2] <- 0.5 + 0.95*abs(0.5-rmix[i])
rmix[i] ~ dbeta( 1, 1 )
}
Пользовательские переменные y
, h
, w
, x1
и x2
.
Отдельная переменная alpha
и матричная переменная theta
имеют нормальные априорные значения.
Переменная mtau
- это первые два элемента из вектора tau
переменных, заданных
mtau <- tau[1:2]
for( i in 1:5 ){
tau[i] ~ dunif(0,1000)
}
Наконец, сообщение об ошибке состоит в том, что узел y[1]
имеет несовместимых родителей.