Copyright	(c) 2015 Jared Tobin
License	MIT
Maintainer	Jared Tobin <jared@jtobin.ca>
Stability	unstable
Portability	ghc
Safe Haskell	None
Language	Haskell2010

Numeric.MCMC.Metropolis

Contents

Re-exported

Description

This implementation uses spherical Gaussian proposals to implement a reliable and computationally inexpensive sampling routine. It can be used as a baseline from which to benchmark other algorithms for a given problem.

The mcmc function streams a trace to stdout to be processed elsewhere, while the metropolis transition can be used for more flexible purposes, such as working with samples in memory.

Synopsis

Documentation

mcmc :: (MonadIO m, PrimMonad m, Traversable f, Show (f Double)) => Int -> Double -> f Double -> (f Double -> Double) -> Gen (PrimState m) -> m () #

Trace n iterations of a Markov chain and stream them to stdout.

>>> let rosenbrock [x0, x1] = negate (5  *(x1 - x0 ^ 2) ^ 2 + 0.05 * (1 - x0) ^ 2)
>>> withSystemRandom . asGenIO $ mcmc 3 1 [0, 0] rosenbrock
0.5000462419822702,0.5693944056267897
0.5000462419822702,0.5693944056267897
-0.7525995304580824,1.2240725505283248

chain :: (PrimMonad m, Traversable f) => Int -> Double -> f Double -> (f Double -> Double) -> Gen (PrimState m) -> m [Chain (f Double) b] #

Trace n iterations of a Markov chain and collect the results in a list.

>>> let rosenbrock [x0, x1] = negate (5  *(x1 - x0 ^ 2) ^ 2 + 0.05 * (1 - x0) ^ 2)
>>> results <- withSystemRandom . asGenIO $ chain 3 1 [0, 0] rosenbrock
>>> mapM_ print results
0.0,0.0
1.4754117657794871e-2,0.5033208261760778
3.8379699517007895e-3,0.24627131099479127

metropolis :: (Traversable f, PrimMonad m) => Double -> Transition m (Chain (f Double) b) #

A generic Metropolis transition operator.

Re-exported

module Data.Sampling.Types

create :: PrimMonad m => m (Gen (PrimState m)) #

Create a generator for variates using a fixed seed.

createSystemRandom :: IO GenIO #

Seed a PRNG with data from the system's fast source of pseudo-random numbers. All the caveats of withSystemRandom apply here as well.

withSystemRandom :: PrimBase m => (Gen (PrimState m) -> m a) -> IO a #

Seed a PRNG with data from the system's fast source of pseudo-random numbers ("/dev/urandom" on Unix-like systems or RtlGenRandom on Windows), then run the given action.

This is a somewhat expensive function, and is intended to be called only occasionally (e.g. once per thread). You should use the Gen it creates to generate many random numbers.

asGenIO :: (GenIO -> IO a) -> GenIO -> IO a #

Constrain the type of an action to run in the IO monad.