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