Model
{
for (i in 1:n)
{
z[i]<-1-alpha/(A[i]*pow(D[i], beta));
theta[i]<-(0.9*step(z[i])+0.1)*exp(-0.5*abs(z[i]));
x[i]~dbern(theta[i]);
}
1/alpha~dunif (0.1, 100);
beta~dunif (0.1, 2);
}