An efficient method of posterior sampling for Poisson INGARCH models

Abstract

We develop an efficient posterior sampling scheme for the Poisson INGARCH models. The proposed method is based on the approximation of the posterior density that exploits the Poisson limit of the negative binomial distribution. It allows us to rewrite the model in a form amenable to Pólya-Gamma data augmentation scheme, which yields simple conditionally Gaussian updates for the autoregressive coefficients. Sampling from the approximate posterior is straightforward via Gibbs-type iterations and remains numerically stable even under strong temporal dependence. Using this sampler as a proposal distribution will enhance the efficiency in Metropolis-Hastings algorithm and adaptive importance sampling. Numerical simulations indicate accurate posterior estimates, high effective sample sizes, and rapidly mixing chains.

Figures (16)

• Figure 1: Sample path and acf of Scenario A1.
• Figure 3: Posterior distribution of the samples fitted to the simulated dataset illustrated in Fig. \ref{['A1']}.
• Figure 4: Posterior distribution of the samples fitted to the simulated dataset illustrated in Fig. \ref{['A2']}.
• Figure 5: Sample path and acf of Scenario B1.
• Figure 7: Posterior distribution of the samples fitted to the simulated dataset illustrated in Fig. \ref{['B1']}.