PTMCMCSampler icon indicating copy to clipboard operation
PTMCMCSampler copied to clipboard

MSTI feature added, temp replaced by beta

Open amandakaris opened this issue 6 months ago • 1 comments

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 avatar Jul 31 '24 15:07 amandakaris