model{
    for(i in 1:n){
      y[i] ~ dnorm(mu[i],tau)
      mu[i] <- b0 +b[1]*X[i,1]+b[2]*X[i,2]+b[3]*X[i,3]+b[4]*X[i,4]+
                  b[5]*X[i,5]+b[6]*X[i,6]+b[7]*X[i,7]+b[8]*X[i,8] 
    }
    b0 ~ dnorm(0,0.001)
    tau ~ dgamma(0.001,0.001)
    for(j in 1:8){
       b[j] ~ dnorm(0,0.001)  
    }
}