Published by the Foundation for Open Access Statistics Editors-in-chief: Bettina Grün, Torsten Hothorn, Rebecca Killick, Edzer Pebesma, Achim Zeileis    ISSN 1548-7660; CODEN JSSOBK
[test by reto]
Authors: Alireza S. Mahani, Mansour T. A. Sharabiani
Title: Multivariate-From-Univariate MCMC Sampler: The R Package MfUSampler
Abstract: The R package MfUSampler provides Markov chain Monte Carlo machinery for generating samples from multivariate probability distributions using univariate sampling algorithms such as the slice sampler and the adaptive rejection sampler. The multivariate wrapper performs a full cycle of univariate sampling steps, one coordinate at a time. In each step, the latest sample values obtained for other coordinates are used to form the conditional distributions. The concept is an extension of Gibbs sampling where each step involves, not an independent sample from the conditional distribution, but a Markov transition for which the conditional distribution is invariant. The software relies on proportionality of conditional distributions to the joint distribution to implement a thin wrapper for producing conditionals. Examples illustrate basic usage as well as methods for improving performance. By encapsulating the multivariate-from-univariate logic, package MfUSampler provides a reliable package for rapid prototyping of custom Bayesian models while allowing for incremental performance optimizations such as taking advantage of conditional independence, and high-performance implementation of function evaluations. Utility functions for MCMC diagnostics as well as sample-based construction of predictive posterior distributions are provided in MfUSampler.

Page views:: 1746. Submitted: 2014-12-28. Published: 2017-06-13.
Paper: Multivariate-From-Univariate MCMC Sampler: The R Package MfUSampler     Download PDF (Downloads: 2263)
MfUSampler_1.0.4.tar.gz: R source package Download (Downloads: 120; 579KB)
v78c01.R: R replication code Download (Downloads: 166; 8KB)

DOI: 10.18637/jss.v078.c01

This work is licensed under the licenses
Paper: Creative Commons Attribution 3.0 Unported License
Code: GNU General Public License (at least one of version 2 or version 3) or a GPL-compatible license.