Print

Print


Hi,
Does anyone know how to code bivariate probit model in WinBugs? Basically, 
the model is as follows:

Y1 and Y2 are observed binary variables.  Z is a vector of latent 
variables Z1 and Z2. Zj>0 if Yj=1, Zj<=0 if Yj=0. (j=1,2).  Z~MVN
(mu,sigma).  

I tried code:
for( i in 1 : N ) {
   Z[i , 1:2] ~ dmnorm(mu[ ], Sigma[ , ]) 
   for(j in 1:2){
    y[i,j]<-step(Z[i,j])
   } 
  }
Winbugs won't compile the above code.  The error message says "multiple 
definitions of node y[1,1]"
 
Anyone has the code or examples? Thanks. 

-------------------------------------------------------------------
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