model{
  theta ~ dgamma(alpha,beta)
  for(i in 1:n){
    y[i] ~ dpois(theta)
  }
  ynew ~ dpois(theta)
}