Bayesian indirect inference for models with intractable normalizing functions

Research output: Contribution to journalArticlepeer-review

Abstract

Inference for doubly intractable distributions is challenging because the intractable normalizing functions of these models include parameters of interest. Previous auxiliary variable MCMC algorithms are infeasible for multi-dimensional models with large data sets because they depend on expensive auxiliary variable simulation. We develop a fast Bayesian indirect algorithm by replacing an expensive auxiliary variable simulation from a probability model with a computationally cheap simulation from a surrogate model. We learn the relationship between the surrogate model parameters and the probability model parameters using Gaussian process approximations. We apply our methods to challenging examples, and illustrate that the algorithm addresses both computational and inferential challenges for doubly intractable distributions. Especially for a large social network model with 10 parameters, we show that our method can reduce computing time from about 2 weeks to 5 hours, compared to the previous method.

Original languageEnglish
Pages (from-to)300-315
Number of pages16
JournalJournal of Statistical Computation and Simulation
Volume91
Issue number2
DOIs
Publication statusPublished - 2021

Bibliographical note

Funding Information:
Jaewoo Park was partially supported by the Yonsei University Research Fund of 2019- 22-0194 and the National Research Foundation of Korea (NRF-2020R1C1C1A01003868). The author is grateful to the anonymous reviewers for their careful reading and valu- able comments.

Publisher Copyright:
© 2020 Informa UK Limited, trading as Taylor & Francis Group.

Copyright:
Copyright 2021 Elsevier B.V., All rights reserved.

All Science Journal Classification (ASJC) codes

  • Statistics and Probability
  • Modelling and Simulation
  • Statistics, Probability and Uncertainty
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Bayesian indirect inference for models with intractable normalizing functions'. Together they form a unique fingerprint.

Cite this