model{
 for(k in 1:38){
   kills[k] ~ dpois(lambda[k]*trials[k])
   log(lambda[k]) <- beta[1]+beta[2]*(age[k]-mean(age[]))+
   beta[3]*(age[k]-mean(age[]))**2+eta[k] }
 for(i in 1:38){ eta[i] ~ dnorm(0,tau) }
 for(i in 1:11){
    log(ExKill[i]) <-  beta[1] + beta[2]*(pages[i]-mean(age[])) + 
    beta[3]*(pages[i]-mean(age[]))**2 + 1/(2*tau)
   d[i] ~ dnorm(0,tau)
   log(r[i]) <- beta[1] + beta[2]*(pages[i]-mean(age[])) + 
   beta[3]*(pages[i]-mean(age[]))**2  + d[i]
 }
 beta[1] ~ dnorm(0,0.0001)
 beta[2] ~ dnorm(0,0.0001)
 beta[3] ~ dnorm(0,0.0001)
 tau ~ dgamma(0.0001,0.0001)
}