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[i] <- 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[i]] + u1[county[i]]*x[i,1] + u2[county[i]]*x[i,2] + u3[county[i]]*x[i,3] + u4[county[i]]*x[i,4] + u5[county[i]]*x[i,5] - z[ID[i]] y[i] ~ dnorm(mu[i], tau) } for (i in 1:220) { z[i] ~ djl.dnorm.trunc(rho,lambda,0,1000) z[i] ~ dnorm(rho,lambda)T(0,100) #For openbugs eff[i] <- exp(-z[i]) }
prior for rho~dnorm(0,0.027)
I would appreciate your help!
Regards,
Daniel
Contents
hide
Best Answer
For truncating, I use the operator I
i.e.:
x ~ dnorm(mean, tau)I(low, high)
Similar Posts:
- Solved – How to estimate the precision of a normal using a Gibbs sampler
- Solved – How to estimate the precision of a normal using a Gibbs sampler
- Solved – get predictions from Winbugs/OpenBUGS
- Solved – get predictions from Winbugs/OpenBUGS
- Solved – the difference between a truncated normal distribution and a half normal distribution in a Stochastic Frontier Analysis