model{
	psi ~ dunif(0,1)
        logit(p[1])<-b[1]
        for(j in 2:K){
          logit(p[j])<-b[1]+sum(absb[2:j])
          absb[j]<-abs(b[j])
        }
        for(j in 1:K){
          b[j]~dnorm(0,0.1)
        }
 	pie[1:K]~ddirch(alpha[])
	for(i in 1:M){
		z[i]~dbin(psi,1)
		class[i]~dcat(pie[])
		zp[i]<-p[class[i]]*z[i]
		x[i] ~ dbin(zp[i],T)
	}
	N<-sum(z[]) 
}