8000 change the subsampling method to use pymbar.timeseries.subsampleCorrelatedData() · Issue #39 · alchemistry/alchemlyb · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

change the subsampling method to use pymbar.timeseries.subsampleCorrelatedData() #39

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
shuail opened this issue Nov 14, 2017 · 1 comment
Assignees
Milestone

Comments

@shuail
Copy link
Collaborator
shuail commented Nov 14, 2017

Proposed to change the subsampling method to avoid using the pandas slicing for subsampling, this issue is corresponding to the pull request #38

The pandas slicing method requires the steps to be integer which is rounded from the statistical inefficiency of a time series. This will sometimes make different time series indistinguishable from each other in terms of their subsampling positions. Proposed to switch to the pymbar subsampleCorrelatedData function.

@orbeckst
Copy link
Member
orbeckst commented Nov 16, 2017

The discussion on PR #38 lead to the following consensus (see comments by @davidlmobley @shuail and myself leading up to #38 (comment)):

  • use pymbar.timeseries.subsampleCorrelatedData() for fractional g or ceil(g)
  • add a keyword to allow the user to switch between conservative=True and False
  • make the default conservative=True

@orbeckst orbeckst changed the title change the subsamping method? change the subsamping method to use pymbar.timeseries.subsampleCorrelatedData() Nov 16, 2017
@orbeckst orbeckst changed the title change the subsamping method to use pymbar.timeseries.subsampleCorrelatedData() change the subsampling method to use pymbar.timeseries.subsampleCorrelatedData() Nov 16, 2017
orbeckst pushed a commit that referenced this issue Nov 22, 2017
- change the subsampling code to use the pymbar module to get the
  subsampled indexes, the original method use the slicing function in
  pandas which requires the slicing step to be integer
orbeckst added a commit that referenced this issue Nov 22, 2017
- kwarg conservative=True (default)
- test that series for statistical ineffiency calculation and the data
  are sampled on the same times (otherwise the slicing would be wrong)
- added tests
- added notes and warning to the docs
- fixes #39
orbeckst added a commit that referenced this issue Nov 23, 2017
change subsampling with statistical_inefficiency() (#39)
@orbeckst orbeckst added this to the release 0.2.0 milestone Jan 17, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants
0