Hi, I am interested in generating replicates from the posterior predictive distribution of a r.v. X, which has been programmed in WinBUGS using the zeros trick. Can someone suggest whether I am doing it correctly? So, zeros.rep[] in the code below are my replicates. K <- 1000 ll[i] <- ...... ## Log-likelihood of X zeros[i] < - 0 zeros[i] ~ dpois(phi[i]) ## Zeros trick phi[i] <- - ll[i] + K zeros.rep[i] ~ dpois(phi[i]) ## Replicates of y[i] ? Thanks a bunch, Dipankar *______________________________ * * Dipankar Bandyopadhyay, PhD * *Associate Professor,* *Division of Biostatistics, School of Public Health,** **University of Minnesota * *Minneapolis, MN 55455 USA* *http://www.biostat.umn.edu/~dipankar*<http://www.biostat.umn.edu/~dipankar> ------------------------------------------------------------------- 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