Use in practice of importance sampling for repeated MCMC for poisson models

Dorota Gajda, Chantal Guihenneuc-Jouyaux, Judith Rousseau, Kerry Mengersen, Darfiana Nur

Research output: Contribution to journalArticlepeer-review

4 Citations (Scopus)

Abstract

The Importance Sampling method is used as an alternative approach to MCMC in repeated Bayesian estimations. In the particular context of numerous data sets, MCMC algorithms have to be called on several times which may become computationally expensive. Since Importance Sampling requires a sample from a posterior distribution, our idea is to use MCMC to generate only a certain number of Markov chains and use them later in the subsequent IS estimations. For each Importance Sampling procedure, the suitable chain is selected by one of three criteria we present here. The first and second criteria are based on the L1 norm of the difference between two posterior distributions and their Kullback-Leibler divergence respectively. The third criterion results from minimizing the variance of IS estimate. A supplementary automatic selection procedure is also proposed to choose those posterior for which Markov chains will be generated and to avoid arbitrary choice of importance functions. The featured methods are illustrated in simulation studies on three types of Poisson model: simple Poisson model, Poisson regression model and Poisson regression model with extra Poisson variability. Different parameter settings are considered.

Original languageEnglish
Pages (from-to)361-383
Number of pages23
JournalElectronic Journal of Statistics
Volume4
DOIs
Publication statusPublished - 2010
Externally publishedYes

Fingerprint

Dive into the research topics of 'Use in practice of importance sampling for repeated MCMC for poisson models'. Together they form a unique fingerprint.

Cite this