flat-mcmc: Painless general-purpose sampling.
flat-mcmc* is a Haskell library for painless, efficient, general-purpose sampling from continuous distributions.
flat-mcmc* uses an ensemble sampler that is invariant to affine transformations of space. It wanders a target probability distribution's parameter space as if it had been "flattened" or "unstretched" in some sense, allowing many particles to explore it locally and in parallel.
In general this sampler is useful when you want decent performance without dealing with any tuning parameters or local proposal distributions.
flat-mcmc* exports an
mcmcfunction that prints a trace to stdout, as well as a
flattransition operator that can be used more generally.
import Numeric.MCMC.Flat import Data.Vector (Vector, toList, fromList) rosenbrock :: Vector Double -> Double rosenbrock xs = negate (5 *(x1 - x0 ^ 2) ^ 2 + 0.05 * (1 - x0) ^ 2) where [x0, x1] = toList xs ensemble :: Ensemble ensemble = fromList [ fromList [negate 1.0, negate 1.0] , fromList [negate 1.0, 1.0] , fromList [1.0, negate 1.0] , fromList [1.0, 1.0] ] main :: IO () main = withSystemRandom . asGenIO $ mcmc 12500 ensemble rosenbrock
|Versions [faq]||0.1.0.0, 1.0.0, 1.0.1, 1.1.1, 1.2.1, 1.2.2, 1.3.0, 1.4.0, 1.4.1, 1.4.2, 1.5.0, 1.5.1, 1.5.2|
|Dependencies||base (<5), mcmc-types (>=1.0.1 && <2), monad-par, monad-par-extras, mwc-probability (>=1.0.1 && <2), pipes (==4.*), primitive, transformers, vector [details]|
|Source repo||head: git clone http://github.com/jtobin/flat-mcmc.git|
|Uploaded||by JaredTobin at 2016-04-06T13:34:10Z|
|Distributions||LTSHaskell:1.5.0, NixOS:1.5.2, Stackage:1.5.0|
|Downloads||6950 total (19 in the last 30 days)|
|Rating||(no votes yet) [estimated by Bayesian average]|
Docs available [build log]
Last success reported on 2016-04-06 [all 1 reports]
For package maintainers and hackage trustees