model
		{
			for( j in 1 : q ) {
				u[j] ~ dnorm(0.0,tau)
			}
			for( j in 1 : p ) {
#				b[j] ~ dnorm(0.0,1.0E-6)
				b[j] ~ dunif(-10,10)
			}			
			for( i in 1 : N ) {			 
			  logit(P[i]) <- inprod(x[i,],b[])+inprod(z[i,],u[])
			  y[i] ~ dbin(P[i],1)
			}
			tau ~ dgamma(0.001,0.001)I(.003,20000)
			sigma <- 1 / sqrt(tau)
			logsigma <- log(sigma)
		}
