model{
    for( i in 1:N){
        y1[i]~dpois(lambda1[i])
        log(lambda1[i])<-delta1[i]       
        delta1[i]~dnorm(mu1,tau1)
    }
            
        mu1~dnorm(0,1.0E-6)
        
        tau1~dgamma(0.01,0.01)
            
}