Hello All,

 

For some reason, OpenBUGS has trouble sampling from a multivariate normal distribution with boundary constraints.

 

The snippet of code below runs in WinBUGS but it does not run in OpenBUGS.

 

for (i in 1:n) {

    Z[i,1:2]        ~ dmnorm(mu[i,1:2],

                             Omega[1:2,1:2]) I(lo[i,1:2],up[i,1:2])

 

The upper and lower bounds (lo[i,1:2],up[i,1:2]) are read in as data.

 

Any suggestions.

 

Sincerely,

 

Ziad Elmously