19,822 research outputs found
BayesBinMix: an R Package for Model Based Clustering of Multivariate Binary Data
The BayesBinMix package offers a Bayesian framework for clustering binary
data with or without missing values by fitting mixtures of multivariate
Bernoulli distributions with an unknown number of components. It allows the
joint estimation of the number of clusters and model parameters using Markov
chain Monte Carlo sampling. Heated chains are run in parallel and accelerate
the convergence to the target posterior distribution. Identifiability issues
are addressed by implementing label switching algorithms. The package is
demonstrated and benchmarked against the Expectation-Maximization algorithm
using a simulation study as well as a real dataset.Comment: Accepted to the R Journal. The package is available on CRAN:
https://CRAN.R-project.org/package=BayesBinMi
A Bayesian Approach to the Detection Problem in Gravitational Wave Astronomy
The analysis of data from gravitational wave detectors can be divided into
three phases: search, characterization, and evaluation. The evaluation of the
detection - determining whether a candidate event is astrophysical in origin or
some artifact created by instrument noise - is a crucial step in the analysis.
The on-going analyses of data from ground based detectors employ a frequentist
approach to the detection problem. A detection statistic is chosen, for which
background levels and detection efficiencies are estimated from Monte Carlo
studies. This approach frames the detection problem in terms of an infinite
collection of trials, with the actual measurement corresponding to some
realization of this hypothetical set. Here we explore an alternative, Bayesian
approach to the detection problem, that considers prior information and the
actual data in hand. Our particular focus is on the computational techniques
used to implement the Bayesian analysis. We find that the Parallel Tempered
Markov Chain Monte Carlo (PTMCMC) algorithm is able to address all three phases
of the anaylsis in a coherent framework. The signals are found by locating the
posterior modes, the model parameters are characterized by mapping out the
joint posterior distribution, and finally, the model evidence is computed by
thermodynamic integration. As a demonstration, we consider the detection
problem of selecting between models describing the data as instrument noise, or
instrument noise plus the signal from a single compact galactic binary. The
evidence ratios, or Bayes factors, computed by the PTMCMC algorithm are found
to be in close agreement with those computed using a Reversible Jump Markov
Chain Monte Carlo algorithm.Comment: 19 pages, 12 figures, revised to address referee's comment
Accelerating MCMC via Parallel Predictive Prefetching
We present a general framework for accelerating a large class of widely used
Markov chain Monte Carlo (MCMC) algorithms. Our approach exploits fast,
iterative approximations to the target density to speculatively evaluate many
potential future steps of the chain in parallel. The approach can accelerate
computation of the target distribution of a Bayesian inference problem, without
compromising exactness, by exploiting subsets of data. It takes advantage of
whatever parallel resources are available, but produces results exactly
equivalent to standard serial execution. In the initial burn-in phase of chain
evaluation, it achieves speedup over serial evaluation that is close to linear
in the number of available cores
The effects of LIGO detector noise on a 15-dimensional Markov-chain Monte-Carlo analysis of gravitational-wave signals
Gravitational-wave signals from inspirals of binary compact objects (black
holes and neutron stars) are primary targets of the ongoing searches by
ground-based gravitational-wave (GW) interferometers (LIGO, Virgo, and
GEO-600). We present parameter-estimation results from our Markov-chain
Monte-Carlo code SPINspiral on signals from binaries with precessing spins. Two
data sets are created by injecting simulated GW signals into either synthetic
Gaussian noise or into LIGO detector data. We compute the 15-dimensional
probability-density functions (PDFs) for both data sets, as well as for a data
set containing LIGO data with a known, loud artefact ("glitch"). We show that
the analysis of the signal in detector noise yields accuracies similar to those
obtained using simulated Gaussian noise. We also find that while the Markov
chains from the glitch do not converge, the PDFs would look consistent with a
GW signal present in the data. While our parameter-estimation results are
encouraging, further investigations into how to differentiate an actual GW
signal from noise are necessary.Comment: 11 pages, 2 figures, NRDA09 proceeding
- …