model {            
            for (i in 1:N2) {
            
                    y[i] ~ dnorm(mean.y[i],tau.y)
                    mean.y[i]<-b0+b1*x1.mcar[i]+b2*x2.mcar[i]+b3*x3.mcar[i]
                    logit(mean.x1[i])<-mu.x1
                    x1.mcar[i]~dbern(mean.x1[i])
                    mean.x2[i]<-mu.x2
                    mean.x3[i]<-mu.x3
                    x2.mcar[i]~dnorm(mean.x2[i],tau.x2)
                    x3.mcar[i]~dnorm(mean.x3[i],tau.x3)                            
            } 
            tau.y<-pow(sigma.y,-2)
            sigma.y~dunif(.00001,100)
            tau.x2<-pow(sigma.x2,-2)
            sigma.x2~dunif(.00001,100)
            tau.x3<-pow(sigma.x3,-2)
            sigma.x3~dunif(.00001,100)
            b0 ~ dnorm(0,.01)
            b1 ~ dnorm(0,.01)
            b2 ~ dnorm(0,.01)
            b3 ~ dnorm(0,.01)
            mu.x1~dnorm(0,.01)
            mu.x2~dnorm(0,.01)
            mu.x3~dnorm(0,.01)
        }