Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Mason Ng
Posts: 25
Joined: May 17 2017
Affiliation: The University of Auckland

Sorry for the many questions - I'm just about mostly set-up for CosmoMC.

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&#40;chain_dir='/Users/masonng/Documents/University/Honours Project/cosmomc/chains/'&#41;
roots = &#91;'base_plikHM_TT_lowl', 'second_planck_best_fit', 'seventh_planck_best_fit'&#93;
#TT_lowl uses both the high-l likelihood from plik
#and the low-l likelihood from Commander
params = &#91;'theta', 'omegabh2', 'omegach2', 'logA', 'ns', 'tau'&#93;
param_3d = None
g.export&#40;'/Users/masonng/Documents/University/Honours Project/tests.pdf'&#41;

Code: Select all

&#91;local host information&#93;&#58;chains masonng&#36; cd second
&#91;local host information&#93;&#58;second masonng&#36; 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
&#91;local host information&#93;&#58;second masonng&#36; pwd
/Users/masonng/Documents/University/Honours Project/cosmomc/chains/second
One difference is that in params.ini, I've used, for the file root:

Code: Select all

file_root = planck_best_fit
I did manually rename the chains (added prefix second_ for example). Is this the cause of the problem? The end goal is to scp the files from the cluster to my desktop (particularly into the /Users/masonng/Documents/University/Honours Project/cosmomc/chains/ folder) and then use Python to generate the joint posterior distributions and parameter posterior distributions as before.

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