# linear regression
model {
for (i in 1:20) {
a[i] ~ dnorm(mu.a[i], tau)
mu.a[i] <- alpha * deltaV[i]
}
#Priors for the regression parameters
alpha ~ dnorm(0.1, 1.0E-6)
#prior for the precision parameter
tau ~ dgamma (1.0E-3, 1.0E-3)
#monitor the standard deviation
sigma <- 1/sqrt(tau)
}
#Data
list(
a=c(-0.49,1.008,-2.302,0.912,-2.354,0.976,-2.238,-0.336,-0.457,-1.994,1.57,-1.91,1.52,-1.906,-0.048,1.507,-1.906,1.488,4.072,1.233),
deltaV=c(0.59,3.47,3.34,1.38,0.13,0.17,0.49,0.36,0.02,0.03,0.04,0.02,0.01,0.02,0.02,0.02,0.01,-1.53,-1.52,-1.38),
)
# initials values
list(alpha=0, tau=1)