Liang et al. (2008) Mixtures of g-Priors for Bayesian Variable Selection

g-priorvariable-selectionbayes-factormodel-selectionbayesian-model-averagingshrinkagemarginal-likelihood

Summary

Zellner's g-prior is the popular conventional prior for Bayesian variable selection because its marginal likelihoods are closed-form functions of the model R2R^2 — but a single fixed gg creates consistency pathologies. Liang, Paulo, Molina, Clyde and Berger study mixtures of g-priors (a prior on gg, integrated out) as a default that resolves these pathologies while keeping tractability. They show the Zellner-Siow Cauchy prior is a special case, introduce the hyper-g and hyper-g/n prior families with marginal likelihoods in closed form via the Gaussian hypergeometric function 2F1{}_2F_1, and establish the theoretical properties (paradox resolution, model-selection consistency) with real and simulated comparisons against fixed-gg, empirical-Bayes, and other default procedures.

Key Claims

Concepts Introduced or Extended

Entities Mentioned

Quotes

"Zellner's g-prior remains a popular conventional prior for use in Bayesian variable selection, despite several undesirable consistency issues … we study mixtures of g-priors as an alternative … that resolve many of the problems with the original formulation, while maintaining the computational tractability that has made the g-prior so popular."

My Take

This is the paper that made the g-prior safe to use as a default for model selection rather than just estimation. Its contribution is precise diagnosis plus a tractable cure: it names the two failure modes of fixed gg (Bartlett's and the information paradox), shows they are exactly what a mixture repairs, and — crucially for adoption — delivers the hyper-g family with closed-form marginals so practitioners lose nothing computationally. It underwrites the marginal-likelihood machinery in BAS (same authorship lineage — Clyde, Berger) and slots directly beside Maruyama-George (2011), which pushes the same fully-Bayes-gg idea into the p>np>n regime. The honest caveat is that "default" still involves a choice among hyper-g, hyper-g/n and Zellner-Siow, whose differences matter most exactly when nn is small.