model
   {
      for( i in 1 : doses ) {
         for( j in 1 : plates ) {
            y[i , j] ~ dpois(mu[i , j])
            log(mu[i , j]) <- alpha + beta * log(x[i] + 10) + 
               gamma * x[i] + lambda[i , j]
            lambda[i , j] ~ dnorm(0.0, tau)   
         }
      }
      alpha ~ dnorm(0.0,1.0E-6)
      beta ~ dnorm(0.0,1.0E-6)
      gamma ~ dnorm(0.0,1.0E-6)
      tau ~ dgamma(0.001, 0.001)
      sigma <- 1 / sqrt(tau)
   }