Непоследовательные родители, использующие dnormmix в JAGS - PullRequest
0 голосов
/ 24 февраля 2020

Примечание. Изначально это сообщение было отправлено на доску обсуждений 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] имеет несовместимых родителей.

...