Displaying 1 result from an estimated 1 matches for "eco3".
Did you mean:
eco
2009 Dec 03
2
Help R2WinBUGS
...model.file="oecd1.bug",
n.chains=3,
n.iter=1000,
debug=TRUE)
# This code works in WinBUGS.
model {
# dummy variable
for (i in 1:N){
D.eco2[i] <- equals(eco[i],2)
D.eco3[i] <- equals(eco[i],3)
}
# model's likelihood
for (i in 1:N){
y[i] ~ dnorm (mu[i], tau.y)
mu[i] <- beta.0 + beta.x*x[i] + beta.eco[2]*D.eco2[i] +
beta.eco[3]*D.eco3[i]
}
#
beta.eco[1] <- 0
# priors
beta.0 ~ dnorm (0,...