Is there a recommended default prior for a covariance matrix? I’m asking because the Prior Choice Recommendations page has contradictory information which is virtually certain to confuse a reader.

I would say there is a consensus now to decompose a covariance matrix into a correlation matrix and something else. But there is less consensus on whether the something else should be standard deviations or variances and less consensus on what the prior should be.

I think people should do what **rstanarm** does, but that is somewhat complicated.

I’ll add that to the wiki, unless someone beats me to it.

Hi,

Would you have some explanation/reference on why decomposing the covariance into a correlation matrix and something else is better? I noticed in my simulations that I get way more precise results, and faster computations using LKJ priors (following the manual’s reparametrization advices) instead of any sort of invert Wishart or Wishart prior (though I have not tried any reparametrization using them). I would like to understand why!

Thank you!

The inverse wishart shows correlation between correlation and variance of the cov, which make it hard to sample from. There is a wide discussion about in literature.

I like Ben’s paper on the topic (with a bunch of co-authors):

http://www.stat.columbia.edu/~gelman/research/unpublished/Visualization.pdf

This is far better than the hack found in the Gelman and Hill regression book, which rescales a Wishart prior!