Skip to main content

The R package MitISEM: Efficient and Robust Simulation Procedures for Bayesian Inference

Grassi, Stefano, Bastürk, Nalan, Hoogerheide, Lennart, Opschoor, Anne, Van Dijk, Herman K. (2017) The R package MitISEM: Efficient and Robust Simulation Procedures for Bayesian Inference. Journal of Statistical Software, 79 (1). ISSN 1548-7660. (doi:10.18637/jss.v079.i01) (KAR id:51624)

PDF Publisher pdf
Language: English


Creative Commons Licence
This work is licensed under a Creative Commons Attribution 4.0 International License.
Download (3MB) Preview
[img]
Preview
PDF Author's Accepted Manuscript
Language: English


Creative Commons Licence
This work is licensed under a Creative Commons Attribution 4.0 International License.
Download (3MB) Preview
[img]
Preview
Official URL
http://dx.doi.org/10.18637/jss.v079.i01

Abstract

This paper presents the R-package MitISEM (mixture of t by importance sampling weighted expectation maximization) which provides an automatic and flexible two-stage method to approximate a non-elliptical target density kernel -- typically a posterior density kernel -- using an adaptive mixture of Student-t densities as approximating density. In the first stage a mixture of Student-t densities is fitted to the target using an expectation maximization (EM) algorithm where each step of the optimization procedure is weighted using importance sampling. In the second stage this mixture density is a candidate density for efficient and robust application of importance sampling or the Metropolis-Hastings (MH) method to estimate properties of the target distribution. The package enables Bayesian inference and prediction on model parameters and probabilities, in particular, for models where densities have multi-modal or other non-elliptical shapes like curved ridges. These shapes occur in research topics in several scientific fields. For instance, analysis of DNA data in bio-informatics, obtaining loans in the banking sector by heterogeneous groups in financial economics and analysis of education's effect on earned income in labor economics. The package MitISEM provides also an extended algorithm, 'sequential MitISEM', which substantially decreases computation time when the target density has to be approximated for increasing data samples.

Item Type: Article
DOI/Identification number: 10.18637/jss.v079.i01
Uncontrolled keywords: finite mixtures, Student-t densities, importance sampling, MCMC, Metropolis-Hastings algorithm, expectation maximization, Bayesian inference, R-software
Subjects: H Social Sciences > HA Statistics
Divisions: Faculties > Social Sciences > School of Economics
Depositing User: Stefano Grassi
Date Deposited: 10 Nov 2015 10:31 UTC
Last Modified: 29 May 2019 16:19 UTC
Resource URI: https://kar.kent.ac.uk/id/eprint/51624 (The current URI for this page, for reference purposes)
  • Depositors only (login required):

Downloads

Downloads per month over past year