Event Title

A Comparison of Programming Languages for Markov chain Monte Carlo Applications

Presenter Information

Bethany Faz

Presentation Type

Poster Presentation/Art Exihibt

College

College of Natural Sciences

Major

Mathematics

Location

Event Center BC

Faculty Mentor

Dr. Jorge Carlos Roman

Start Date

5-18-2017 11:00 AM

End Date

5-18-2017 12:00 PM

Abstract

Markov chain Monte Carlo (MCMC) methods are utilized to generate samples from intractable distributions. For example, one MCMC method commonly used is the Gibbs sampler, which draws approximate samples from conditional distributions. The prevalence of the Gibbs sampler and other MCMC methods in various disciplines has created a need for efficient programs to run MCMC algorithms. For our research, we used the following five programming languages to compare computation times for Gibbs samplers in several statistical models: R, C++ (using the Rcpp package in R), Just Another Gibbs Sampler (JAGS), Julia, and MATLAB. From our results, we see that each language had its benefits and limitations. However, with its simple syntax and fast computation time, we contend that Julia is the most optimal programming language for MCMC algorithms.

Share

COinS
 
May 18th, 11:00 AM May 18th, 12:00 PM

A Comparison of Programming Languages for Markov chain Monte Carlo Applications

Event Center BC

Markov chain Monte Carlo (MCMC) methods are utilized to generate samples from intractable distributions. For example, one MCMC method commonly used is the Gibbs sampler, which draws approximate samples from conditional distributions. The prevalence of the Gibbs sampler and other MCMC methods in various disciplines has created a need for efficient programs to run MCMC algorithms. For our research, we used the following five programming languages to compare computation times for Gibbs samplers in several statistical models: R, C++ (using the Rcpp package in R), Just Another Gibbs Sampler (JAGS), Julia, and MATLAB. From our results, we see that each language had its benefits and limitations. However, with its simple syntax and fast computation time, we contend that Julia is the most optimal programming language for MCMC algorithms.