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)
}