model
   {
      for( i in 1 : N )
      {
      x[i] ~ dgpz(alpha, theta)      
      }      
   # Prior distributions of the model parameters         
   # Gamma Prior      
         alpha ~ dgamma(0.001, 0.001)
         theta~ dgamma(0.001, 0.001)               
   }