Hi all. I just wanted to let youall know that problems can occur with seemingly-noninformative inverse-gamma prior distributions for variance components. In many cases, it's more reliable to simply use uniform prior distributions on the standard deviation parameters. We have an example in Appendix C of the second edition of our book, "Bayesian Data Analysis." To get this appendix, just go to the bugs.R page (http://www.stat.columbia.edu/~gelman/bugsR/) and near the top of the page, there's a place to click to download Appendix C. For example, in a simple 1-way data structure (the "8 schools" example from "Bayesian Data Analysis"), we use the following model: model { for (j in 1:J){ y[j] ~ dnorm (theta[j], tau.y[j]) theta[j] ~ dnorm (mu.theta, tau.theta) tau.y[j] <- pow(sigma.y[j], -2) } mu.theta ~ dnorm (0, 1.0E-6) tau.theta <- pow(sigma.theta, -2) sigma.theta ~ dunif (0, 1000) } This is a uniform prior distribution on the the sd parameter (sigma.theta) and works fine. An alternative uses the inverse-gamma prior, replacing the last 2 lines in the above model by, tau.theta ~ dgamma (1, 1) sigma.theta <- 1/sqrt(tau.theta) This does not work well at all (see Figure C.3 on page 597). - Andrew Gelman P.S. I'm not trying to knock the inverse-gamma model in general. It can be very effective, especially when applied hierarchically to several variance parameters. I'm just making a recommendation for what to do when a noninformative distribution is required. ------------------------------------------------------------------- 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