model
    {
        for (i in 1 : 100) {
            cases[i] ~ dpois(mu[i])
            mu[i]<-pop[i]*lambda[i]
            log(lambda[i]) <- theta[i]
            theta[i] ~ dnorm(-6.5, 10)
        }
        
    }