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