model{
  theta ~ dexp(0.05)
  for(i in 1:n){
    y[i] ~ dpois(theta)
  }
}