Dear BUGS-users,
who can help me generating truncated valuations?
My problem is similar to the Birats-case but the
estimates for beta[2] have to be restricted by
(0.02, 0.66)
I tried as follows:
model
{
for( i in 1 : N ) {
beta[i , 1] ~ dmnorm(mu.beta[], R[ , ])
beta[i , 2] ~ dmnorm(mu.beta[], R[ , ])I(0.02,
0.66)
for( j in 1 : T ) {
Y[i , j] ~ dnorm(mu[i , j], tauC)
mu[i , j] <- beta[i , 1] + beta[i , 2] * x[j]
}
}
mu.beta[1:2] ~ dmnorm(mean[],prec[ , ])
R[1:2 , 1:2] ~ dwish(Omega[ , ], 2)
tauC ~ dgamma(0.001, 0.001)
sigma <- 1 / sqrt(tauC)
}
It didn't work. Where is the mistake, or how can
I construct the "one-by-one" sampling method from
Devroye (1986)? Would this help?
Regards,
Marc Werkmeister
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|