list(d=c(0,1,0,1,0,1,0,1,0,1), ws=c(1,2,3,4,5,6,7,8,9,10), N=10)  #data

list(beta0=0.5, beta=0.5, sigma=0)#  3õÖμ


model
        {
                for (i in 1 : N) {
                        d[i]   ~ dbern(p[i])        
                        logit(p[i]) <-  beta0 + beta * ws[i] + sigma       
                }                                       
                beta0 ~ dnorm(0.0, 0.00001);
                beta   ~ dnorm(0.0, 0.00001);
                sigma ~ dnorm(0.0,0.00001);
        }

数据在这