model {
  for (i in 1:N) {
  ataques[i] ~ dnorm(FirstGroup[X[i]], tau[X[i]])
  }
  for (i in 1:n) {
  FirstGroup[i] ~ dnorm(mq, tau1)
  tau[i] <- 1/(sdw[i]*sdw[i])
  sdw[i] ~ dunif(0, 100)
  }

  mq ~ dnorm(0, 0.010E-6)
  tau1 <- 1/(sdb*sdb)
  sdb ~ dunif(0,100)
  diff2 <- FirstGroup[C]-FirstGroup[D]
  z1 <- step(diff2)

  }