model {
for (i in 1:23) {
y[i] ~ dnorm(mu[i], tau)
mu[i] <- beta0 + beta[1]*MAN[i]
}
beta0 ~ dnorm(0, 0.0001)
beta[1] ~ dnorm(0, 0.0001)
tau <- 1/pow(sigma, 2)
sigma ~ dunif(0, 100)
dummy <- AUTO[1] + DIST[1]
}
Inits:
list(beta0 = 0, beta = c(0), sigma = 1)
Data:
MAN[] AUTO[] y[] DIST[]
-15.76 1.09 3.19 1
0.98 0.62 -3.45 1
3.71 0.61 0.04 1
-5.37 -1.01 6.62 1
-10.23 -0.76 3.61 1
-8.32 1.91 2.67 1
-7.80 0.40 -2.45 1
6.77 -1.71 9.31 1
-8.81 -0.76 15.29 1
-9.56 -1.34 3.68 1
-2.06 -1.71 8.63 2
-0.76 -1.82 10.82 2
-6.30 -4.91 -0.50 2
39.40 -2.65 -11.00 2
-10.79 0.11 2.05 2
-8.16 0.52 11.80 2
-2.82 -2.54 -2.02 2
-16.19 -0.07 0.94 3
-11.00 -0.83 4.42 3
-14.60 0.98 -0.86 3
-17.96 -3.41 -0.92 3
0.76 2.97 2.61 3
-10.77 2.92 1.58 3
END
node mean sd MC error 2.5% median 97.5% start sample
beta[1] -0.1761 0.1096 0.001085 -0.3921 -0.1765 0.04073 5001 10000
beta0 1.97 1.362 0.01394 -0.7489 1.963 4.632 5001 10000
sigma 5.873 0.9813 0.01054 4.297 5.74 8.136 5001 10000