I have run (as a test) CosmoMC on the cluster and got back some chains. I am now trying to use GetDist on Python to analyze the chains. However, I am getting the error that the chain is not found even though I've basically gone back to an earlier thread (http://cosmocoffee.info/viewtopic.php?t=2836&highlight=) I asked in and followed the instructions:
Below is the Python code used
Code: Select all
g = plots.getSubplotPlotter(chain_dir='/Users/masonng/Documents/University/Honours Project/cosmomc/chains/')
roots = ['base_plikHM_TT_lowl', 'second_planck_best_fit', 'seventh_planck_best_fit']
#TT_lowl uses both the high-l likelihood from plik
#and the low-l likelihood from Commander
params = ['theta', 'omegabh2', 'omegach2', 'logA', 'ns', 'tau']
param_3d = None
g.triangle_plot(roots, params, plot_3d_with_param=param_3d, filled_compare=True, shaded=False)
g.export('/Users/masonng/Documents/University/Honours Project/tests.pdf')
Code: Select all
[local host information]:chains masonng$ cd second
[local host information]:second masonng$ ls
second_planck_best_fit.converge_stat second_planck_best_fit_2.data second_planck_best_fit_4.txt second_planck_best_fit_7.data
second_planck_best_fit.inputparams second_planck_best_fit_2.log second_planck_best_fit_5.chk second_planck_best_fit_7.log
second_planck_best_fit.likelihoods second_planck_best_fit_2.txt second_planck_best_fit_5.data second_planck_best_fit_7.txt
second_planck_best_fit.paramnames second_planck_best_fit_3.chk second_planck_best_fit_5.log second_planck_best_fit_8.chk
second_planck_best_fit.ranges second_planck_best_fit_3.data second_planck_best_fit_5.txt second_planck_best_fit_8.data
second_planck_best_fit_1.chk second_planck_best_fit_3.log second_planck_best_fit_6.chk second_planck_best_fit_8.log
second_planck_best_fit_1.data second_planck_best_fit_3.txt second_planck_best_fit_6.data second_planck_best_fit_8.txt
second_planck_best_fit_1.log second_planck_best_fit_4.chk second_planck_best_fit_6.log
second_planck_best_fit_1.txt second_planck_best_fit_4.data second_planck_best_fit_6.txt
second_planck_best_fit_2.chk second_planck_best_fit_4.log second_planck_best_fit_7.chk
[local host information]:second masonng$ pwd
/Users/masonng/Documents/University/Honours Project/cosmomc/chains/second
Code: Select all
file_root = planck_best_fit