model{
  for(i in 1:30){
    Conc[i] ~ dnorm(mu[i], tau)
    mu[i] <- beta[Treatment[i]] + b[Dog[i]]
  }
  for(i in 1:10){ b[i] ~ dnorm(0,taub) }
  for(i in 1:3){ beta[i] ~ dnorm(0,0.01) } 
  taub <- 1/(sigmab*sigmab)
  tau ~ dgamma(0.01,0.01)
  sigmab ~ dunif(0,100)
  sigma <- 1/sqrt(tau)
  diff12 <- beta[2] - beta[1]
  diff13 <- beta[3] - beta[1]
  diff23 <- beta[3] - beta[2]
}