Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory
We provide a Bayesian analysis of ocean circulation based on data collected in the South Atlantic Ocean. The analysis incorporates a reaction-diffusion partial differential equation that is not solvable in closed form. This leads to an intractable likelihood function. We describe a novel Markov chai...
Main Authors: | , |
---|---|
Format: | Article in Journal/Newspaper |
Language: | unknown |
Subjects: | |
Online Access: | http://hdl.handle.net/10.1080/01621459.2014.914439 |
id |
ftrepec:oai:RePEc:taf:jnlasa:v:109:y:2014:i:507:p:944-954 |
---|---|
record_format |
openpolar |
spelling |
ftrepec:oai:RePEc:taf:jnlasa:v:109:y:2014:i:507:p:944-954 2023-05-15T18:21:00+02:00 Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory Radu Herbei L. Mark Berliner http://hdl.handle.net/10.1080/01621459.2014.914439 unknown http://hdl.handle.net/10.1080/01621459.2014.914439 article ftrepec 2020-12-04T13:33:02Z We provide a Bayesian analysis of ocean circulation based on data collected in the South Atlantic Ocean. The analysis incorporates a reaction-diffusion partial differential equation that is not solvable in closed form. This leads to an intractable likelihood function. We describe a novel Markov chain Monte Carlo approach that does not require a likelihood evaluation. Rather, we use unbiased estimates of the likelihood and a Bernoulli factory to decide whether or not proposed states are accepted. The variates required to estimate the likelihood function are obtained via a Feynman-Kac representation. This lifts the common restriction of selecting a regular grid for the physical model and eliminates the need for data preprocessing. We implement our approach using the parallel graphic processing unit (GPU) computing environment. Article in Journal/Newspaper South Atlantic Ocean RePEc (Research Papers in Economics) |
institution |
Open Polar |
collection |
RePEc (Research Papers in Economics) |
op_collection_id |
ftrepec |
language |
unknown |
description |
We provide a Bayesian analysis of ocean circulation based on data collected in the South Atlantic Ocean. The analysis incorporates a reaction-diffusion partial differential equation that is not solvable in closed form. This leads to an intractable likelihood function. We describe a novel Markov chain Monte Carlo approach that does not require a likelihood evaluation. Rather, we use unbiased estimates of the likelihood and a Bernoulli factory to decide whether or not proposed states are accepted. The variates required to estimate the likelihood function are obtained via a Feynman-Kac representation. This lifts the common restriction of selecting a regular grid for the physical model and eliminates the need for data preprocessing. We implement our approach using the parallel graphic processing unit (GPU) computing environment. |
format |
Article in Journal/Newspaper |
author |
Radu Herbei L. Mark Berliner |
spellingShingle |
Radu Herbei L. Mark Berliner Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
author_facet |
Radu Herbei L. Mark Berliner |
author_sort |
Radu Herbei |
title |
Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
title_short |
Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
title_full |
Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
title_fullStr |
Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
title_full_unstemmed |
Estimating Ocean Circulation: An MCMC Approach With Approximated Likelihoods via the Bernoulli Factory |
title_sort |
estimating ocean circulation: an mcmc approach with approximated likelihoods via the bernoulli factory |
url |
http://hdl.handle.net/10.1080/01621459.2014.914439 |
genre |
South Atlantic Ocean |
genre_facet |
South Atlantic Ocean |
op_relation |
http://hdl.handle.net/10.1080/01621459.2014.914439 |
_version_ |
1766200005728665600 |