ＢＵＧＳのモデルを載せるの忘れてました。。

model
{
for(i in 1:N)
{
lwage[i] ~ dnorm (mu[i], tau)
mu[i] <- b[1] + b[2]*exper[i] + b[3]*expersq[i] + b[4]*educ[i] + b[5]*age[i] + b[6]*kidslt6[i] + b[7]*kidsge6[i] + b[8]*city[i]
}

for(k in 1:8)
{
b[k] ~ dnorm(0, 0.0001)
}

sigma ~ dunif(0, 1000)
tau <-1/(sigma*sigma)
}