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

MSTI feature added, temp replaced by beta #56

Open
wants to merge 9 commits into
base: master
Choose a base branch
from

Conversation

amandakaris
Copy link

Model-Switch Thermodynamic Integration (MSTI) feature has been added, allowing users to sample across the posteriors of two models, instead of sampling between each model's prior and posterior, individually. Beta=0 is the posterior of model 0 and beta=1 is the posterior of model 1. A beta in between 0 and 1 is a mixture of the two models’ posterior distributions. Changes include the ability to give PTMCMC a tuple of prior and likelihood functions, a new parameter for a tuple of lists of indices for each model’s parameters used to sift through the combined set of parameters in p0, and a parameter to indicate the preferred ladder shape (only geometric or linear as of right now). When using MSTI, four additional columns are written in the chain file -- the log posterior and log likelihood of the two models being compared.

Instances of 1/temperature have been replaced by beta. Changes include PTMCMC accepting beta max and beta min (while still accepting Tmin and Tmax), chain files named with beta instead of temperature, and a parameter to indicate whether the temperature naming convention should be use instead.

Temperature scaling has been commented out since for hot chains it produces massive and frequently invalid jumps.

amandakaris and others added 8 commits July 30, 2024 13:31
Model-Switch Thermodynamic Integration (MSTI) feature has been added, allowing users to sample across the posteriors of two models, instead of sampling between each model's prior and posterior, individually. Beta=0 is the posterior of model 0 and beta=1 is the posterior of model 1. A beta in between 0 and 1 is a mixture of the two models’ posterior distributions. Changes include the ability to give PTMCMC a tuple of prior and likelihood functions, a parameter of a tuple of lists of indices for each model’s parameters in MSTI used to sift through the combined set of parameters in p0, and a parameter to indicate the preferred ladder shape (only geometric or linear as of right now). When using MSTI, four additional columns are written in the chain file -- the log posterior and log likelihood of the two models being compared.  

Instances of 1/temperature have been replaced by beta. Changes include PTMCMC accepting beta max and beta min (while still accepting Tmin and Tmax), chain files named with beta instead of temperature, and a parameter to indicate whether the temperature naming convention should be use.

Temperature scaling has been commented out since for hot chains it produces massive and frequently invalid jumps.
E501, E711, E262, E203, F821, E266
@kdolum kdolum self-requested a review August 6, 2024 20:32
Copy link
Collaborator

@kdolum kdolum left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good to me.

@kdolum
Copy link
Collaborator

kdolum commented Aug 7, 2024

Looks good to me. @vhaasteren , @AaronDJohnson do you have any comments?

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.

2 participants