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)
}