Skip to content
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

add stepping stone method for calculating evidences #20

Draft
wants to merge 1 commit into
base: dev
Choose a base branch
from

Conversation

meyers-academic
Copy link

@meyers-academic meyers-academic commented Apr 20, 2022

Add stepping stone method for calculating evidences. Uses la_forge.slices.SliceCore object (which is still only on dev branch). Include "moving block bootstrap" for calculating errors on evidences.

Example usage is fairly simple:

from la_forge import slices, bfacs
m3a_slice = slices.SlicesCore(label='pat', slicedirs=['12p5yr_ti/12p5yr_ti_3a/'], pt_chains=True, pars2pull=['lnlike'])
evidence, evidence_error = bfacs.calc_stepping_stone(m3a_slice, bs_iterations=100)

Still to finish:

  • Unit test comparing method to analytic evidence calculation
  • Proper documentation and warnings about MBB method for getting error (may need to maximize over block)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant