CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

[CosmoMC] (custom-run) chain not found
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Mason Ng



Joined: 17 May 2017
Posts: 25
Affiliation: The University of Auckland

PostPosted: June 19 2017  Reply with quote

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:
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:
[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


One difference is that in params.ini, I've used, for the file root:
Code:

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.
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1293
Affiliation: University of Sussex

PostPosted: June 19 2017  Reply with quote

You'll need to set chain_dir to a list of directories to search for the chains.
Back to top
View user's profile [ Hidden ] Visit poster's website
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.