itadd mcmc methods - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations Err mx1.adamsgaard.dk 70 hgit clone git://src.adamsgaard.dk/cosmo URL:git://src.adamsgaard.dk/cosmo mx1.adamsgaard.dk 70 1Log /src/log.gph mx1.adamsgaard.dk 70 1Files /src/files.gph mx1.adamsgaard.dk 70 1Refs /src/refs.gph mx1.adamsgaard.dk 70 1README /src/file/README.md.gph mx1.adamsgaard.dk 70 1LICENSE /src/file/LICENSE.gph mx1.adamsgaard.dk 70 i--- Err mx1.adamsgaard.dk 70 1commit 4769723624cf47e9a1c827ebe1388971f8616e2c /src/commit/4769723624cf47e9a1c827ebe1388971f8616e2c.gph mx1.adamsgaard.dk 70 1parent c61d24deab6ce81dcb4e8b538cab571d035beb2e /src/commit/c61d24deab6ce81dcb4e8b538cab571d035beb2e.gph mx1.adamsgaard.dk 70 hAuthor: Anders Damsgaard URL:mailto:anders.damsgaard@geo.au.dk mx1.adamsgaard.dk 70 iDate: Fri, 27 Nov 2015 15:34:25 +0100 Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 iadd mcmc methods Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 iDiffstat: Err mx1.adamsgaard.dk 70 i M pages/methods.html | 26 ++++++++++++++++++++++++-- Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 i1 file changed, 24 insertions(+), 2 deletions(-) Err mx1.adamsgaard.dk 70 i--- Err mx1.adamsgaard.dk 70 1diff --git a/pages/methods.html b/pages/methods.html /src/file/pages/methods.html.gph mx1.adamsgaard.dk 70 it@@ -52,13 +52,35 @@ Err mx1.adamsgaard.dk 70 i
Err mx1.adamsgaard.dk 70 i

Err mx1.adamsgaard.dk 70 i Markov-Chain Monte Carlo (MCMC) basics

Err mx1.adamsgaard.dk 70 i-

bla bla

Err mx1.adamsgaard.dk 70 i+

The inversion problem of turning observed TCN Err mx1.adamsgaard.dk 70 i+ concentrations into erosion histories is handled using a Err mx1.adamsgaard.dk 70 i+ conventional Metropolis-Hastings MCMC approach where the Err mx1.adamsgaard.dk 70 i+ model parameters are constrained between fixed model Err mx1.adamsgaard.dk 70 i+ parameter bounds specified by the user. Erosion rates Err mx1.adamsgaard.dk 70 i+ (εint, εgla), which Err mx1.adamsgaard.dk 70 i+ may vary over several orders of magnitude, are specified Err mx1.adamsgaard.dk 70 i+ with uniform probability across the logarithmic parameter Err mx1.adamsgaard.dk 70 i+ interval. The temporal parameter (tdegla) Err mx1.adamsgaard.dk 70 i+ and climate record threshold value Err mx1.adamsgaard.dk 70 i+ (δ18Othreshold) is specified Err mx1.adamsgaard.dk 70 i+ with uniform probability across the linear Err mx1.adamsgaard.dk 70 i+ parameter interval. The user specifies the bounds Err mx1.adamsgaard.dk 70 i+ of the model parameters, which define the model space. Err mx1.adamsgaard.dk 70 i+

Err mx1.adamsgaard.dk 70 i
Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 i
Err mx1.adamsgaard.dk 70 i

Err mx1.adamsgaard.dk 70 i What is a MCMC walker?

Err mx1.adamsgaard.dk 70 i-

bla bla

Err mx1.adamsgaard.dk 70 i+

Err mx1.adamsgaard.dk 70 i+ forward responses are computed based on an initial set of Err mx1.adamsgaard.dk 70 i+ model parameters that is proposed using the Err mx1.adamsgaard.dk 70 i+ Metropolis-Hastings technique. A burn-in phase of 1000 Err mx1.adamsgaard.dk 70 i+ iterations is first used to make a crude initial search of Err mx1.adamsgaard.dk 70 i+ the model space. This step is followed by a more detailed Err mx1.adamsgaard.dk 70 i+ and local search of the model space based on the best-fit Err mx1.adamsgaard.dk 70 i+ model parameters from the burn-in phase. Err mx1.adamsgaard.dk 70 i+

Err mx1.adamsgaard.dk 70 i
Err mx1.adamsgaard.dk 70 i Err mx1.adamsgaard.dk 70 i
Err mx1.adamsgaard.dk 70 .