model 3PLM;
{

   for (k in 1:I) {
      a[k] ~ dlnorm(0, 4);
      b[k] ~ dnorm(0, .25);
      c[k] ~ dbeta(5,17);
   }

   for (j in 1:N) {
      for (k in 1:I) {
         p[j,k] <- c[k]+((1-c[k])*(exp(a[k]*(theta[j]-b[k]))/(1+exp(a[k]*(theta[j]-b[k])))));
         r[j,k] ~ dbern(p[j,k]);
       }
       theta[j] ~ dnorm(0,1);
   } 
}
list(N=50,  I=5, r=structure(.Data= c(
 1, 1, 1, 1, 0,
 1, 1, 1, 1, 1,
 1, 1, 0, 0, 0,
 1, 0, 0, 1, 0,
 1, 0, 1, 0, 1,
 1, 1, 0, 0, 1,
 1, 0, 0, 0, 0,
 1, 1, 0, 0, 0,
 1, 0, 0, 0, 0,
 1, 1, 1, 0, 0,
 0, 1, 0, 0, 0,
 1, 0, 0, 0, 0,
 1, 1, 1, 1, 1,
 1, 1, 1, 1, 1,
 1, 1, 0, 1, 1,
 1, 1, 1, 0, 1,
 1, 1, 1, 0, 1,
 0, 1, 0, 0, 0,
 1, 1, 0, 0, 1,
 1, 0, 0, 0, 1,
 1, 0, 0, 1, 0,
 1, 0, 0, 0, 0,
 1, 1, 1, 0, 1,
 1, 0, 0, 1, 1,
 1, 1, 1, 1, 0,
 1, 0, 0, 0, 0,
 1, 0, 0, 0, 0,
 1, 1, 1, 1, 1,
 0, 0, 1, 1, 0,
 1, 1, 1, 0, 1,
 1, 0, 1, 0, 0,
 1, 0, 1, 1, 0,
 1, 1, 1, 1, 0,
 1, 0, 1, 1, 0,
 1, 1, 1, 0, 0,
 1, 0, 1, 0, 0,
 1, 1, 1, 1, 0,
 1, 0, 1, 0, 0,
 1, 1, 0, 0, 0,
 0, 0, 1, 0, 0,
 0, 0, 0, 0, 0,
 1, 1, 1, 0, 0,
 1, 0, 0, 1, 1,
 1, 0, 1, 1, 1,
 1, 1, 0, 0, 0,
 1, 1, 1, 1, 0,
 1, 1, 1, 1, 0,
 1, 0, 1, 0, 0,
 0, 1, 0, 1, 0,
 1, 1, 1, 1, 1), .Dim=c(50,5)))