My old rig has fallen, so I just finished basically formatting everything and reinstalling CosmoMC (and all deps).
I have now run 3 separate runs with the different datasets.
Just as a reference:
Code: Select all
#general settings
#Bicep-Keck-Planck, varying cosmological parameters
DEFAULT(batch2/BKPlanck.ini)
#Planck 2015, default just include native likelihoods (others require clik)
#DEFAULT(batch2/plik_dx11dr2_HM_v18_TT.ini)
DEFAULT(batch2/lowTEB.ini)
DEFAULT(batch2/lowl.ini)
DEFAULT(batch2/lensing.ini)
#Other Likelihoods
#DEFAULT(batch2/BAO.ini)
#DEFAULT(batch2/WiggleZ_MPK.ini)
#DEFAULT(batch2/MPK.ini)
#DEFAULT(batch2/WL.ini)
#general settings
DEFAULT(batch2/common.ini)
#e.g. to vary r in addition to standard 6:
#(for r>0 also need compute_tensors=T)
compute_tensors = T
param[r] = 0.06 0.0 0.5 0.04 0.04
param[logA] = 3.1 1.5 4.5 0.04 0.04
param[ns] = 0.96 0.9 1.05 0.04 0.04
param[nrun]=-0.09 -0.5 0.5 0.04 0.04
param[nrunrun]=0
param[tau]=0.07
#limits[r02]=0 N
limits[r]=0 N
The results for [tex]n_{run}[/tex] are always slightly positive, where in my previous experience, when [tex]n_{run}[/tex] is used (and [tex]n_{run,run}[/tex] is set to 0 ) , one usually gets a slightly negative value.
This is also the case with the most recent Planck collaboration analysis.
Has the convention for CosmoMC changed between the last stable version and the current master branch?
If not, any advise on what is happening?
Thanks
-- Ira