model
{
   for(i in 1:n) 
   {
      logit(mu[i]) <- sbeta0 + sbeta1*sbirthweight[i]
      BPD[i] ~ dbern(mu[i])	
   }	
   sbeta0 ~ dnorm(0,1.0E-8)
   sbeta1 ~ dnorm(0,1.0E-8)
   beta1 <- sbeta1/sdBirthweight
   beta0 <- sbeta0 - beta1*meanBirthweight
}