model{

for(i in 1:N){
	Nfrat[i] ~ dpois(lambda[i])
	log(lambda[i]) <- log(Npop[i]) + beta0 + beta.ses[Ses3b[i]] + beta.dur*dur[i] + psi[IndexConc[i]]
}

for(j in 1:278){
	psi[j] ~ dnorm(0,tau.psi)
}

tau.psi ~ dgamma(0.5,0.005)
beta0 ~ dnorm(0,0.001)
beta.ses[1] <- 0
beta.ses[2] ~ dnorm(0,0.001)
beta.ses[3] ~ dnorm(0,0.001)
beta.dur ~ dnorm(0,0.001)

}
