Displaying 1 result from an estimated 1 matches for "sigma_dep".
Did you mean:
sigma_d
2008 Jul 23
1
R2WinBUGS problem
...with hope and some anxiety.
Best wishes
Troels
cat("
model {
for (i in 1:126) {
dep[i] ~ dnorm(mu.dep[i],tau.dep)
mu.dep[i] <- alpha[ID[i]]*NAK.cor[i]+beta[ID[i]]*TBW.cor[i]
fit[i] <- mu.dep[i]+CONST[i]
}
tau.dep <- pow(sigma_dep,-2)
sigma_dep ~ dunif(0,100)
for (j in 1:16) {
alpha[lev[j]] <- AB[lev[j],1]
beta[lev[j]] <- AB[lev[j],2]
AB[lev[j],1:2] ~ dmnorm(AB.hat[lev[j],],tau.AB[,])I(lower[],
upper[])
AB.hat[lev[j],1] <- mu.alpha
AB.hat[lev[j],2] <- m...