New getdist no longer produces data for plots?

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Sunny Vagnozzi
Posts: 55
Joined: August 15 2016
Affiliation: Kavli Institute for Cosmology (KICC), University of Cambridge

New getdist no longer produces data for plots?

Post by Sunny Vagnozzi » November 07 2019

Hi Antony and all,

I've run a couple of chains with the latest version of CosmoMC and analyzed them with the latest getdist as usual:

Code: Select all

python python/GetDist.py distparams.ini chains/root
Getdist works fine in producing the usual useful files (e.g. root.margestats) and convergence diagnostics. However, I also noticed the message
WARNING:root:plot_data_dir is no longer supported by getdist, value ignored
As a result, Getdist seems to no longer produce the data used to produce the contour plots. For instance, files such as root_p_param.dat containing the (normalized) 1D posterior density for the parameter "param".

Is it indeed the case that these files are no longer produced (or have I done something wrong in using Getdist), and if so is there a way to revert to the previous behaviour where these files were produced, other than using an older version of Getdist? While the old Getdist used to produce a lot of ancillary files, some of them were actually useful beyond plotting purposes (e.g. the root_p_param.dat files I mentioned).

Cheers,
Sunny

Antony Lewis
Posts: 1936
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: New getdist no longer produces data for plots?

Post by Antony Lewis » November 07 2019

It was removed, you should be able to make all relevant plots and get relevant data via the python interface, e.g. using
https://getdist.readthedocs.io/en/latest/mcsamples.html#getdist.mcsamples.MCSamples.get1DDensity

Jun Wang
Posts: 2
Joined: March 03 2022
Affiliation: School of Physics and Astronomy, Sun Yat-Sen University

Re: New getdist no longer produces data for plots?

Post by Jun Wang » April 06 2022

You could write a python script

Code: Select all

from getdist import plots
g=plots.get_single_plotter(chain_dir=r''./PLA'')
samples=g.sample_analyser.samples_for_root('base_plikHM_TTTEEE_lowl_lowE_lensing')
print(samples.getMargeStats())
print(samples.getLikeStats())
print(samples.getConvergeTests())
print(samples.getInlineLatex("omegabh2",limit=1))
You could also use getdist-gui >> Data to check those informations

Post Reply