grd349 / PBjam

A repo for our peak baggin code and tips on jam
MIT License
17 stars 6 forks source link

CPnest default settings need to be changed #254

Closed nielsenmb closed 1 month ago

nielsenmb commented 3 years ago

At the moment the maxmcmc and Nlive points values are probably too small for a high dimensional parameter space that we are using it for in asy_peakbag.

@alexelthomas did some testing and found that for high dimensions (~10), and low maxmcmc (~100) produced a normally distributed marginalized posterior for parameters that had uniform priors and likelihoods. The reason for this is most likely because the walkers that CPnest uses in the sampling process are started in a tight ball, and a low maxmcmc doesn't allow them to explore the parameter space sufficiently.

Chris suggested that parallelization should work for CPnest despite it not being beneficial for EMCEE. Also, apparently Dynesty is the new hot nested sampler in Grav Waves and better maintained compared to CPnest.

alexelthomas commented 3 years ago

Chris also suggested a minimum of 1000 live points and probably more like 2000