hnuts

No U-Turn Sampling in Haskell.
git clone git://git.jtobin.io/hnuts.git
Log | Files | Refs | README | LICENSE

commit 0dab29b17a260d9d18e1a3a23ea99056e839e0af
parent 25c15b9f2eaaa13dc844d6ad85a25a17944cdfa3
Author: Jared Tobin <jared@jtobin.ca>
Date:   Tue, 15 Oct 2013 16:06:42 +1300

Remove Debug.Trace stuff.

Diffstat:
Msrc/Numeric/MCMC/NUTS.hs | 4+---
1 file changed, 1 insertion(+), 3 deletions(-)

diff --git a/src/Numeric/MCMC/NUTS.hs b/src/Numeric/MCMC/NUTS.hs @@ -12,8 +12,6 @@ import System.Random.MWC import System.Random.MWC.Distributions hiding (gamma) import Statistics.Distribution.Normal -import Debug.Trace - type Parameters = [Double] type Density = Parameters -> Double type Gradient = Parameters -> Parameters @@ -152,7 +150,7 @@ nutsKernelDualAvg lTarget glTarget e eAvg h m daParams t g = do logEm = mu daParams - sqrt (fromIntegral m) / gamma daParams * hm logEbarM = (1 - zeta) * log eAvg + zeta * logEm - trace (show eAvgNext) $ return (eNext, eAvgNext, hNext, nextPosition) + return (eNext, eAvgNext, hNext, nextPosition) -- | A single iteration of NUTS. nutsKernel