Comments (5)
A paper from David van der Spoel's group truncates the window when the empirical ACF goes below 0.05 http://cmb.bio.uni-goettingen.de/pub/Hub_JCTC2010.pdf
from pyhmc.
For batch means, batch sizes ~n^(1/3) are apparently optimal (page 10) http://projecteuclid.org/download/pdfview_1/euclid.aos/1266586622. They also talk about the tukey-hanning window. Code is in CRAN http://cran.r-project.org/web/packages/mcmcse/mcmcse.pdf
from pyhmc.
In case you haven't seen Geyer's chapter from the Handbook of Markov chain Monte Carlo, it's has some nice discussion of this as well in Section 1.10 "Variance Estimation".
I'd love to see the initial convex sequence (ICS) estimator be implemented as well (discussed in Section 1.10.2). The canonical reference for this and related methods seems to be:
Geyer, C. J. 1992. Practical Markov chain Monte Carlo (with discussion). Statistical Science, 7:473–511.
PDF available here.
from pyhmc.
Yes, I've read the paper and cited it in the IPS docs. Adding the ICS is a good idea too.
from pyhmc.
Also there's another ACF from http://www.sciencedirect.com/science/article/pii/0370269389915633 that I just saw, where you truncate the summation after a fixed number of terms (e.g. 20), and then add in a remainder based on the ratio of the final two values included in the sum, which looks like an analytical correction based on the assumption that the terms are decaying exponentially.
from pyhmc.
Related Issues (6)
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from pyhmc.