Electronic Journal of Statistics
- Electron. J. Statist.
- Volume 12, Number 2 (2018), 3295-3311.
Geometric ergodicity of Pólya-Gamma Gibbs sampler for Bayesian logistic regression with a flat prior
The logistic regression model is the most popular model for analyzing binary data. In the absence of any prior information, an improper flat prior is often used for the regression coefficients in Bayesian logistic regression models. The resulting intractable posterior density can be explored by running Polson, Scott and Windle’s (2013) data augmentation (DA) algorithm. In this paper, we establish that the Markov chain underlying Polson, Scott and Windle’s (2013) DA algorithm is geometrically ergodic. Proving this theoretical result is practically important as it ensures the existence of central limit theorems (CLTs) for sample averages under a finite second moment condition. The CLT in turn allows users of the DA algorithm to calculate standard errors for posterior estimates.
Electron. J. Statist., Volume 12, Number 2 (2018), 3295-3311.
Received: February 2018
First available in Project Euclid: 5 October 2018
Permanent link to this document
Digital Object Identifier
Primary: 60J05: Discrete-time Markov processes on general state spaces
Secondary: 62F15: Bayesian inference
Wang, Xin; Roy, Vivekananda. Geometric ergodicity of Pólya-Gamma Gibbs sampler for Bayesian logistic regression with a flat prior. Electron. J. Statist. 12 (2018), no. 2, 3295--3311. doi:10.1214/18-EJS1481. https://projecteuclid.org/euclid.ejs/1538705039