model{
#
#  Stevens
#  Breyer
#  Ginsburg
#  Souter
#  Kennedy
#  OConnor
#  Rehnquist
#  Scalia
#  Thomas
#
for (i in 1:9){
# Change #
	  for (j in 1:510){
		y[i,j] ~ dbern(p[i,j])
		logit(p[i,j]) <- (x[i]-b[j])*(x[i]-b[j])-(x[i]-q[j])*(x[i]-q[j])
	     }
	  }
   
  ## priors
  for(k in 1:9) {x[k] ~ dnorm(0.0, 1.0)}

# Change #	
	for(j in 1:510){
	   b[j] ~ dnorm(0.0,1.0)
	   q[j] ~ dnorm(0.0,1.0)
	}

}