choderalab / pymbar

Python implementation of the multistate Bennett acceptance ratio (MBAR)
http://pymbar.readthedocs.io
MIT License
238 stars 92 forks source link

Do correlated timeseries bias the PMF or only the uncertainty estimates? #208

Closed John-Jumper closed 9 years ago

John-Jumper commented 9 years ago

The documentation for MBAR strenuously points out the need for uncorrelated input configurations. I am only interested in unbiased estimates of observables for many simulations and do not need to measure uncertainties (these are handled in another part of my procedure). Will correlated data bias expectation values and observables? If it does bias observables, am I okay as long as the correlation time is the roughly the same for all thermodynamic states?

Best, John Jumper

mrshirts commented 9 years ago

As long as you have good uncertainties somewhere (like boostrapping), then correlated data will not bias the expectations PROVIDED you are decorrelating over time. If you discover the correlation time is 50% of your simulation, you probably haven't visited all of the relevant states, and you will have bias simply because of lack of sampling. So I certainly would suggest looking at the autocorrelation time just to be sure.

jchodera commented 9 years ago

Our experience has been that using correlated data to compute expectations does not lead to significant bias, though it can be much slower than using uncorrelated data.

John-Jumper commented 9 years ago

Great, thanks. I have good control over statistical uncertainty later in the pipeline, although I will definitely have to monitor equilibration.