I'm new at using cosmomc and I have some problems with it. I'm currently using the 2015 version. My goal is to produce the plot of w - H0 for CMB and BAO likelihoods using a certain w(α) model and parametrization.

More specifically the test.ini file I'm using is,

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.03 0 2 0.04 0.04
#high for new runs
MPI_Max_R_ProposeUpdate = 30
propose_matrix= planck_covmats/base_BAO_TT_lowTEB_plik.covmat
```

Code: Select all

```
#param[omegabh2] = 0.0221
param[omegabh2] = 0.023 0.013 0.033 0.0001 0.0001
#param[omegach2] = 0.12
param[omegach2] = 0.12 0.001 0.99 0.001 0.0005
#param[theta] = 1.0411
param[theta] = 1.0411 0.5 10 0.0004 0.0002
#param[tau] = 0.09
param[tau] = 0.09 0.01 0.8 0.01 0.005
num_massive_neutrinos=1
param[mnu] = 0.06
param[meffsterile] = 0
param[omegak] = 0
param[w] = -1 -2 0 0.001 0.001
param[nrun] = 0
param[nrunrun] = 0
param[r] = 0
param[wa] = 0
param[nnu] = 3.046
param[yhe] = 0.24
param[alpha1] = 0
param[deltazrei] = 0.5
param[Alens] = 1
param[Alensf]=-1
param[fdm] = 0
#param[ns] = 0.96
param[ns] = 0.96 0.7 1.3 0.004 0.002
#log[10^10 A_s]
#param[logA] = 3.1
param[logA] = 3.1 1.7 5 0.001 0.001
```

I'm open to any and all suggestions as to where my mistake is.

EDIT: Could it be because the default sample number is not enough for 2 likelihoods?

Thank you in advance