Comments (5)
Was this resolved by #470 when the MBAR object is created an carried through, or did you want to actually store the f_k
array to file to be used as an initial guess next time analyze
has been called?
from yank.
At some point, I experimented with storing f_k
each iteration, initializing the MBAR estimate from the previous iteration (if one was available). There were some practical iterations:
- It could slow things down unless we figure out how to do this asynchronously (e.g. using dynamic process management in
mpi4py
with MPI-2) - Re-analyzing the trajectory for equilibration start and statistical inefficiency is slow, so we could either use all the data for the online estimate of
f_k
or update the estimates of equilibration start and statistical inefficiency eachnupdate
iterations. - Estimates after 1 or a few iterations don't make much sense, though I think they can still be computed
from yank.
In my own experience, MBAR initializes faster if you feed it an f_k
even from just a few samples, and even if you add more states on after the fact. That said, the more samples you have, the slower it gets, so I really think we would have to do some type of parallelization to have a practical implementation.
re-analyzing the trajectory for equilibration start and statistical inefficiency is slow,
Yes, but for the on-the-fly analysis we can probably skip that part. It will change the answer, but more by bias than anything. We just have to tell the user that the numbers are inaccurate and to run proper statistical analysis on it later.
As a side thought, would the timesereies
module benefit from some attempted cython speedup (I say attempted so as not to make it mandatory for pymbar
).
Estimates after 1 or a few iterations don't make much sense, though I think they can still be computed
Both are true, low sampling is garbage for stats. It could just be as simple as adding a if iteration < min_iteration
or something.
from yank.
Depends on analyze.py refactor, touched for 1.0
from yank.
Implemented as part of the online analysis in #712, and cached for the API calls on #660
from yank.
Related Issues (20)
- Load() missing one argument
- pymbar 4 compatibility HOT 1
- I got the following error. HOT 6
- multichain receptor for 'restrained_receptor_atoms'
- Issue with Ligand Coordinates During Absolute-Binding Free Energy
- NetCDF upstream changes breaking Yank (and can't update OpenMMtools)
- Failures due to API changes in `openmmtools` HOT 3
- problem during absolute binding energy calculation HOT 3
- support for python 3.10? HOT 2
- Issue in Yank's calls to OpenMMTools MultiStateReporter HOT 1
- Question: How to measure simulation time and how much is necessary?
- Silent GAFF2 failures can happen easily
- Stalling on multiple SLURM nodes HOT 8
- support Membrane Protein Systems ?
- Varying free energy estimates depending on sampler used
- unallowed value hdac11-ft895
- leap returned non-zero exit status 31
- Clarification on why vacuum simulation needed for hydration free energy HOT 2
- Calculate Free binding energy without schrodingers suite
- YANK Support for Multiple Conformational States
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 yank.