model{
for(i in 1:26){
e[i]<-(mexp[i]+fexp[i])/10
count[i]~dpois(mu[i])
log(mu[i])<-log(e[i])+bet0+bet1*dep[i]+u[i]
u[i]~dnorm(0,tauU)
smr1[i]<-SMR[i]
}
bet0~dunif(-1000,1000)
bet1~dnorm(0,tauB1)
tauU<-3
tauB1<-3

}

list( mexp=c( 20.1743,14.0029,8.453,20.3617,20.1327,17.5309,21.1234,23.7393,18.5391,19.754,19.7998,16.0856,21.374,13.1552,22.362,17.0192,14.71,23.3018,13.6725,19.8907,23.6726,6.2829,9.2587,18.7117,23.7979,5.1447 ), 
fexp=c( 8.04,6.9128,2.9689,7.0274,9.4097,7.4324,11.6288,9.9962,7.5676,7.6736,9.0426,7.6734,7.2171,4.8009,10.9576,7.95,9.7088,12.6952,6.9551,6.1989,8.9045,3.3625,5.4941,9.3766,13.1685,1.7782 ),
 dep=c( 3.04,-1.36,-2.51,3,0.96,0.38,1,-0.17,2.74,1.89,0.69,1.19,-3.23,-3.92,-0.57,-2.39,-0.2,-0.56,-0.77,-2.74,0.69,-3.82,-3.03,0.95,1.04,-3.55 ), 
count=c( 3,4,1,2,6,2,2,2,2,2,3,3,2,1,3,1,3,2,3,3,1,1,3,1,3,1 ), 
SMR=c( 1.06329,1.91244,0.87551,0.73022,2.03098,0.80118,0.61065,0.59285,0.76609,0.72919,1.04014,1.26268,0.69952,0.55691,0.90037,0.40049,1.22856,0.5556,1.45436,1.14988,0.30696,1.03676,2.03351,0.35602,0.81155,1.44448 ) )
##### inits ##################
list(bet0=-0.1,bet1=0.1,u=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0))