Ошибки в модели ошибок

Question

Я пытаюсь запустить эту модель BUGS, но сталкиваюсь с некоторыми ошибками. Любая помощь будет принята с благодарностью.

model {
  for (i in 1:n){
    y[i] ~ dnorm (a[state[i]] + theta*treat[i] + b*hispanic, tau.y)
  }
    theta ~ dnorm (0, .0001)
    b ~ dnorm (0, 1000)
  for (j in 1:J){
    a[j] ~ rnorm (mu.a, tau.a^2)
  }
  mu.a ~ dnorm (0, .0001)
  tau.a <- pow (sigma.a, -2)
  sigma.a ~ dunif (0, 100)
  tau.y <- pow (sigma.y, -2)
  sigma.y <- dunif (0, 100)
}