I want to calculate posterior distributions in a case where I have a four-dimensional normal prior and a likelihood that is a bivariate binomial distribution: Y1~binomial(n,p1) Y2~binomial(n,p2) log(-log(p1))=a1-exp(-b1)*x log(-log(p2))=a2-exp(-b2)*x (a1,b1,a2,b2) ~ four-dimensional normal distribution Y1, Y2, x and n are observed. The posterior distribution of (a1,b1,a2,b2) is of interest. As far as I can see it is not possible to fit this model in classic BUGS (running on a Linux platform)? Or ...? Can anyone suggest another program/code for solving this problem?? Preferably for the Linux platform, but other platforms are of interest, too. Regards, Marianne ------------------------------------------------------------------- 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