model {


  for (i in 1:N) {
    Y[i] ~ dbern(p[i])
    logit(p[i]) <- a0[pid[i]] + a1[pid[i]]*Age[i]
  }


  for (j in 1:J) {
    a0[j] ~ dnorm(b0, tau02inv)
    a1[j] ~ dnorm(b1, tau12inv)
  }


  b0 ~ dnorm(0, 0.0001)

  b1 ~ dnorm(0, 0.0001)

  tau02inv <- pow(tau0,-2)
  tau12inv <- pow(tau1,-2)

  tau0 ~ dunif(0, 5)
  tau1 ~ dunif(0, 5)

}