model {

	psi ~ dunif(0,1)
	for(i in 1:4){
	  b[i]~dnorm(0,0.1)
	  g[i]~dnorm(0,0.1)
	}
        tau~dgamma(0.001,0.001)
        pie[1]~dunif(0,1)
        pie[2]~dunif(0,1)

	for(i in 1:M){
	   z[i]~dbin(psi,1)
           logit(p[i])<-b[1]+b[2]*sex[i]+b[3]*age[i]+b[4]*weight[i]
           mu[i]<-g[1]+g[2]*sex[i]+g[3]*age[i]
           weight[i]~dnorm(mu[i],tau)
           
           # Model for sex and age
           sex[i]~dbern(pie[1])
           age[i]~dbern(pie[2])
            
           zp[i]<-p[i]*z[i]
	   x[i] ~ dbin(zp[i],T)
	}
	N<-sum(z[]) 
	
}