Bayesian Inference for Polya Inverse Gamma Models

05/29/2019
by   Christopher Glynn, et al.
0

Probability density functions that include the gamma function are widely used in statistics and machine learning. The normalizing constants of gamma, inverse gamma, beta, and Dirichlet distributions all include model parameters as arguments in the gamma function; however, the gamma function does not naturally admit a conjugate prior distribution in a Bayesian analysis, and statistical inference of these parameters is a significant challenge. In this paper, we construct the Polya-inverse Gamma (P-IG) distribution as an infinite convolution of Generalized inverse Gaussian (GIG) distributions, and we represent the reciprocal gamma function as a scale mixture of normal distributions. As a result, the P-IG distribution yields an efficient data augmentation strategy for fully Bayesian inference on model parameters in gamma, inverse gamma, beta, and Dirichlet distributions. To illustrate the applied utility of our data augmentation strategy, we infer the proportion of overdose deaths in the United States attributed to different opioid and prescription drugs with a Dirichlet allocation model.

READ FULL TEXT
research
06/03/2021

Bayesian Inference for Gamma Models

We use the theory of normal variance-mean mixtures to derive a data augm...
research
05/01/2021

Bayesian Inference of a Dependent Competing Risk Data

Analysis of competing risks data plays an important role in the lifetime...
research
03/12/2019

The Inverse first passage time method for a two compartment model as a tool to relate Inverse Gaussian and Gamma spike distributions

In a previous paper (Lansky, Sacerdote, Zucca (2016)) we related the sto...
research
08/23/2016

Softplus Regressions and Convex Polytopes

To construct flexible nonlinear predictive distributions, the paper intr...
research
06/18/2015

Dependent Multinomial Models Made Easy: Stick Breaking with the Pólya-Gamma Augmentation

Many practical modeling problems involve discrete data that are best rep...
research
11/23/2022

A Generator for Generalized Inverse Gaussian Distributions

We propose a new generator for the generalized inverse Gaussian (GIG) di...
research
07/21/2020

Evaluating the Informativeness of the Besag-York-Mollié CAR Model

The use of the conditional autoregressive framework proposed by Besag, Y...

Please sign up or login with your details

Forgot password? Click here to reset