I am estimating a stochastic frontier with a mixed model. So far the half normal distribution worked good but I need a truncated normal distribution. It does not work, and I receive the error „Expected collection operator c“. I am using R2WinBUGS and as you can see in the model I have tried OpenBUGS and WinBUGS. Any suggestion? My model looks like:
for (i in 1:N)
{
mu <- alpha + x[i,1]*beta[1] + x[i,2]*beta[2] + x[i,3]*beta[3] +
x[i,4]*beta[4] + x[i,5]*beta[5] + u0[county] + u1[county]*x[i,1] +
u2[county]*x[i,2] + u3[county]*x[i,3] + u4[county]*x[i,4] +
u5[county]*x[i,5] - z[ID]
y ~ dnorm(mu, tau)
}
for (i in 1:220) {
z ~ djl.dnorm.trunc(rho,lambda,0,1000)
z ~ dnorm(rho,lambda)T(0,100) #For openbugs
eff <- exp(-z)
}
prior for rho~dnorm(0,0.027)
I would appreciate your help! Regards