Poisson log-linear models are ubiquitous in many applications, and one of the most popular approaches for parametric count regression. In the Bayesian context, however, there are no sufficient specific computational tools for efficient sampling from the posterior distribution of parameters, and standard algorithms, such as random walk Metropolis-Hastings or Hamiltonian Monte Carlo algorithms, are typically used. Herein, we developed an efficient Metropolis-Hastings algorithm and importance sampler to simulate from the posterior distribution of the parameters of Poisson log-linear models under conditional Gaussian priors with superior performance with respect to the state-of-the-art alternatives. The key for both algorithms is the introduction of a proposal density based on a Gaussian approximation of the posterior distribution of parameters. Specifically, our result leverages the negative binomial approximation of the Poisson likelihood and the successful Polya-gamma data augmentation scheme. Via simulation, we obtained that the time per independent sample of the proposed samplers is competitive with that obtained using the successful Hamiltonian Monte Carlo sampling, with the Metropolis-Hastings showing superior performance in all scenarios considered. Supplementary materials for this article are available online.

D'Angelo, L., Canale, A. (2023). Efficient Posterior Sampling for Bayesian Poisson Regression. JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS, 32(3), 917-926 [10.1080/10618600.2022.2123337].

Efficient Posterior Sampling for Bayesian Poisson Regression

D'Angelo, L
;
2023

Abstract

Poisson log-linear models are ubiquitous in many applications, and one of the most popular approaches for parametric count regression. In the Bayesian context, however, there are no sufficient specific computational tools for efficient sampling from the posterior distribution of parameters, and standard algorithms, such as random walk Metropolis-Hastings or Hamiltonian Monte Carlo algorithms, are typically used. Herein, we developed an efficient Metropolis-Hastings algorithm and importance sampler to simulate from the posterior distribution of the parameters of Poisson log-linear models under conditional Gaussian priors with superior performance with respect to the state-of-the-art alternatives. The key for both algorithms is the introduction of a proposal density based on a Gaussian approximation of the posterior distribution of parameters. Specifically, our result leverages the negative binomial approximation of the Poisson likelihood and the successful Polya-gamma data augmentation scheme. Via simulation, we obtained that the time per independent sample of the proposed samplers is competitive with that obtained using the successful Hamiltonian Monte Carlo sampling, with the Metropolis-Hastings showing superior performance in all scenarios considered. Supplementary materials for this article are available online.
Articolo in rivista - Articolo scientifico
Count data; Horseshoe prior; Importance sampling; Log-linear models; Metropolis-Hastings; Pólya-gamma;
English
24-ott-2022
2023
32
3
917
926
reserved
D'Angelo, L., Canale, A. (2023). Efficient Posterior Sampling for Bayesian Poisson Regression. JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS, 32(3), 917-926 [10.1080/10618600.2022.2123337].
File in questo prodotto:
File Dimensione Formato  
DAngelo-2023-J Computat Graph Stats-VoR.pdf

Solo gestori archivio

Tipologia di allegato: Publisher’s Version (Version of Record, VoR)
Licenza: Tutti i diritti riservati
Dimensione 1.42 MB
Formato Adobe PDF
1.42 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/10281/396651
Citazioni
  • Scopus 7
  • ???jsp.display-item.citation.isi??? 3
Social impact