getdist.py cannot export *.margestats,etc
-
Tian Qiu
- Posts: 12
- Joined: June 01 2017
- Affiliation: University of Science and Technology of China
Post
by Tian Qiu » March 18 2018
I have added some parameters to CosmoMC and it finished running. However, when I want to analyze the data via GetDist.py, it cannot export the likestats and margestats files. The problem shows below,
Code: Select all
Number of chains used = 4
var(mean)/mean(var), remaining chains, worst e-value: R-1 = 0.00634
RL: Thin for Markov: 55
RL: Thin for indep samples: 57
RL: Estimated burn in steps: 265 ( 87 rows)
using 69328 rows, 62 parameters; mean weight 3.04083487191, tot weight 210815.0
Approx indep samples (N/corr length): 2758.0
Equiv number of single samples (sum w)/max(w): 3833.0
Effective number of weighted samples (sum w)^2/sum(w^2): 37368
Best fit sample -log(Like) = 1229.291000
mean(-Ln(like)) = 1241.424637
-Ln(mean like) = 1236.414926
Traceback (most recent call last):
File "python/GetDist.py", line 357, in <module>
main(args)
File "python/GetDist.py", line 194, in main
mc._setDensitiesandMarge1D(writeDataToFile=not no_plots and plot_data_dir, meanlikes=plot_meanlikes)
File "/home/Tian/program/0-CosmoMC/python/getdist/mcsamples.py", line 2077, in _setDensitiesandMarge1D
meanlikes=meanlikes)
File "/home/Tian/program/0-CosmoMC/python/getdist/mcsamples.py", line 1265, in get1DDensityGridData
if paramrange <= 0: raise MCSamplesError('Parameter range is <= 0: ' + par.name)
getdist.mcsamples.MCSamplesError: Parameter range is <= 0: thetad
I am not sure where the problem occured.
Thank you very much.
-
Antony Lewis
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
-
Contact:
Post
by Antony Lewis » March 19 2018
Suggests maybe the parameters are mixed up, I'd check the .paramnames, .ranges and chain files (and see why thetad is wrong).
-
Tian Qiu
- Posts: 12
- Joined: June 01 2017
- Affiliation: University of Science and Technology of China
Post
by Tian Qiu » March 20 2018
Antony Lewis wrote:Suggests maybe the parameters are mixed up, I'd check the .paramnames, .ranges and chain files (and see why thetad is wrong).
Thanks for your reply.
I have canceled the running of Omegabh2, Omegach2 and theta to accelerate the running and added 5 more parameters, calculated the tansor-scalar ratio. Does it matter?
I have uploaded the two files.
N55k0.05-9.paramnames
N55k0.05-9.ranges
The chain files is too large to upload. This is the first five rows.
Code: Select all
0.1000000E+01 0.7038408E+04 0.9330470E-01 0.6153850E+00 0.5911068E-01 0.1501714E+02 0.2100146E+01 0.2010086E+01 0.5920694E-01 0.3203897E+01 0.1566531E+01 0.9998750E+00 0.7617339E+02 0.2787716E-01 0.2807137E+01 0.2293336E+03 0.1790998E+02 0.2631529E+02 0.8849017E+02 0.1138627E+01 0.6493253E+01 0.1206891E+02 0.1961263E+02 0.9621628E+02 0.6597487E-01 0.3369599E-01 0.7139685E-01 0.9093314E-01 0.1943887E+00 0.6607681E+00 0.2185664E+00 0.1851540E+00 0.3010104E+00 0.2202591E+00 0.4482271E+00 0.1795984E+01 0.9996879E+00 0.9946687E+00 0.6718994E+02 0.6838066E+00 0.3161934E+00 0.1427451E+00 0.6451439E-03 0.9591038E-01 0.8194297E+00 0.4607742E+00 0.6144689E+00 0.9996767E+00 0.2494488E+01 0.1139069E+02 0.2219795E+01 0.1841914E+01 0.1287712E+04 0.5580961E+04 0.2319546E+04 0.8017290E+03 0.2268662E+03 0.9600000E+00 0.2452681E+00 0.2465941E+00 0.2642598E+01 0.1381948E+02 0.1090262E+04 0.1446375E+03 0.1041309E+01 0.1388997E+02 0.1059322E+04 0.1473915E+03 0.1403399E+00 0.1611637E+00 0.3395751E+04 0.1036418E-01 0.8138427E+00 0.4498719E+00 0.7136085E-01 0.9280797E+02 0.1393361E+04 0.6772173E+00 0.4776219E+00 0.6084566E+00 0.1365313E+00 0.1406799E+00 0.1179698E+01 0.6859150E-01 0.3253391E+00 0.2699559E+00 0.4820857E+02 0.1398452E+05 0.2201244E+02 0.8777867E-01 0.9265765E+00 0.2861379E+01 0.1240028E+01 0.1696297E+02 0.1405474E+05 0.5115763E+01
0.1000000E+01 0.6032280E+04 0.7957435E-01 0.6150554E+00 0.5914784E-01 0.1501313E+02 0.2100818E+01 0.2010467E+01 0.6707525E-01 0.3226953E+01 0.1564720E+01 0.9950958E+00 0.4450664E+02 0.7129291E+00 0.8105211E+01 0.1777296E+03 0.3251322E+02 0.6327366E+02 0.1338806E+03 0.1092477E+01 0.1097519E+02 0.9739167E+01 0.2315989E+02 0.1164276E+03 0.8109535E-01 0.3925033E-01 0.9580603E-01 0.8891009E-01 0.2081114E+00 0.3919916E+00 0.2568888E+00 0.2357087E+00 0.4799147E+00 0.1857674E+00 0.5382995E+00 0.2274546E+01 0.9994271E+00 0.9946802E+00 0.6718994E+02 0.6838066E+00 0.3161934E+00 0.1427451E+00 0.6451439E-03 0.9591038E-01 0.8184628E+00 0.4602305E+00 0.6137439E+00 0.9984971E+00 0.2497795E+01 0.1017609E+02 0.2219795E+01 0.1893196E+01 0.1319027E+04 0.5746690E+04 0.2347093E+04 0.8239561E+03 0.2333956E+03 0.9600000E+00 0.2452681E+00 0.2465941E+00 0.2642598E+01 0.1381948E+02 0.1090262E+04 0.1446375E+03 0.1041309E+01 0.1388997E+02 0.1059322E+04 0.1473915E+03 0.1403399E+00 0.1611637E+00 0.3395751E+04 0.1036418E-01 0.8138427E+00 0.4498719E+00 0.7136085E-01 0.9280797E+02 0.1393361E+04 0.6772173E+00 0.4770574E+00 0.6077387E+00 0.1365401E+00 0.1407021E+00 0.1181276E+01 0.6933288E-01 0.3258528E+00 0.2779099E+00 0.4882879E+02 0.1193688E+05 0.2628909E+02 0.8777867E-01 0.9265765E+00 0.2861379E+01 0.1240028E+01 0.4744716E+02 0.1201200E+05 0.5115763E+01
0.1000000E+01 0.5491155E+04 0.6897591E-01 0.6148492E+00 0.5917426E-01 0.1500967E+02 0.2100230E+01 0.2009844E+01 0.6794995E-01 0.3235347E+01 0.1567871E+01 0.9864056E+00 0.3173089E+02 0.9804691E+00 0.6214583E+01 0.1815517E+03 0.1869505E+02 0.6008731E+02 0.1484988E+03 0.1157677E+01 0.1295497E+02 0.9842983E+01 0.8779803E+01 0.8737904E+02 0.8634765E-01 0.3636871E-01 0.7906379E-01 0.9811919E-01 0.2655550E+00 0.2977050E+00 0.1950853E+00 0.2546530E+00 0.3952159E+00 0.1873733E+00 0.4744679E+00 0.1972243E+01 0.9992038E+00 0.9910233E+00 0.6718994E+02 0.6838066E+00 0.3161934E+00 0.1427451E+00 0.6451439E-03 0.9591038E-01 0.8173233E+00 0.4595897E+00 0.6128894E+00 0.9971070E+00 0.2495660E+01 0.9190369E+01 0.2219795E+01 0.1933754E+01 0.1340031E+04 0.5858512E+04 0.2382985E+04 0.8397117E+03 0.2378998E+03 0.9600000E+00 0.2452681E+00 0.2465941E+00 0.2642598E+01 0.1381948E+02 0.1090262E+04 0.1446375E+03 0.1041309E+01 0.1388997E+02 0.1059322E+04 0.1473915E+03 0.1403399E+00 0.1611637E+00 0.3395751E+04 0.1036418E-01 0.8138427E+00 0.4498719E+00 0.7136085E-01 0.9280797E+02 0.1393361E+04 0.6772173E+00 0.4763926E+00 0.6068925E+00 0.1364964E+00 0.1406585E+00 0.1179008E+01 0.6986135E-01 0.3251148E+00 0.2832207E+00 0.4917851E+02 0.1081762E+05 0.3137045E+02 0.8777867E-01 0.9265765E+00 0.2861379E+01 0.1240028E+01 0.7902403E+02 0.1089817E+05 0.5115763E+01
0.8000000E+01 0.4870358E+04 0.7037061E-01 0.6157047E+00 0.5920063E-01 0.1501671E+02 0.2100068E+01 0.2009402E+01 0.7710049E-01 0.3240423E+01 0.1557970E+01 0.9863944E+00 0.3166919E+02 0.2031408E+00 0.5942209E+01 0.1852994E+03 0.3031205E+00 0.3219802E+02 0.1386859E+03 0.1063818E+00 0.1184527E+02 0.1288806E+02 0.7692999E+01 0.8865802E+02 0.1164125E+00 0.3869284E-01 0.6307779E-01 0.9887212E-01 0.2662149E+00 0.3195717E+00 0.2004348E+00 0.2116982E+00 0.4202694E+00 0.1324725E+00 0.4662539E+00 0.1893711E+01 0.9978488E+00 0.9909347E+00 0.6718994E+02 0.6838066E+00 0.3161934E+00 0.1427451E+00 0.6451439E-03 0.9591038E-01 0.8191851E+00 0.4606366E+00 0.6142855E+00 0.9993783E+00 0.2491350E+01 0.9324646E+01 0.2219795E+01 0.1928367E+01 0.1334334E+04 0.5828056E+04 0.2437377E+04 0.8381463E+03 0.2370805E+03 0.9600000E+00 0.2452681E+00 0.2465941E+00 0.2642598E+01 0.1381948E+02 0.1090262E+04 0.1446375E+03 0.1041309E+01 0.1388997E+02 0.1059322E+04 0.1473915E+03 0.1403399E+00 0.1611637E+00 0.3395751E+04 0.1036418E-01 0.8138427E+00 0.4498719E+00 0.7136085E-01 0.9280797E+02 0.1393361E+04 0.6772173E+00 0.4774779E+00 0.6082749E+00 0.1364489E+00 0.1406231E+00 0.1177453E+01 0.6976528E-01 0.3246095E+00 0.2819929E+00 0.4922631E+02 0.9561826E+04 0.3075131E+02 0.8777867E-01 0.9265765E+00 0.2861379E+01 0.1240028E+01 0.9379786E+02 0.9641803E+04 0.5115763E+01
0.2000000E+01 0.4812519E+04 0.6749549E-01 0.6153008E+00 0.5921701E-01 0.1501806E+02 0.2099383E+01 0.2009478E+01 0.8241080E-01 0.3207296E+01 0.1547025E+01 0.9872390E+00 0.2716322E+02 0.3227418E+00 0.5192261E+01 0.1964299E+03 0.8573368E+00 0.3395110E+02 0.1435689E+03 0.9884349E+00 0.1193140E+02 0.1221385E+02 0.7298816E+01 0.9046884E+02 0.1103312E+00 0.4027035E-01 0.5813927E-01 0.1012253E+00 0.2572858E+00 0.3295250E+00 0.2181826E+00 0.2154254E+00 0.4282159E+00 0.1213092E+00 0.4436232E+00 0.1867950E+01 0.9979328E+00 0.9904908E+00 0.6718994E+02 0.6838066E+00 0.3161934E+00 0.1427451E+00 0.6451439E-03 0.9591038E-01 0.8194110E+00 0.4607636E+00 0.6144549E+00 0.9996539E+00 0.2490900E+01 0.9049988E+01 0.2219795E+01 0.1939488E+01 0.1340605E+04 0.5860536E+04 0.2457607E+04 0.8430676E+03 0.2384385E+03 0.9600000E+00 0.2452681E+00 0.2465941E+00 0.2642598E+01 0.1381948E+02 0.1090262E+04 0.1446375E+03 0.1041309E+01 0.1388997E+02 0.1059322E+04 0.1473915E+03 0.1403399E+00 0.1611637E+00 0.3395751E+04 0.1036418E-01 0.8138427E+00 0.4498719E+00 0.7136085E-01 0.9280797E+02 0.1393361E+04 0.6772173E+00 0.4776093E+00 0.6084426E+00 0.1364446E+00 0.1406259E+00 0.1177365E+01 0.6990758E-01 0.3245811E+00 0.2835942E+00 0.4946242E+02 0.9450907E+04 0.3115665E+02 0.8777867E-01 0.9265765E+00 0.2861379E+01 0.1240028E+01 0.8839599E+02 0.9531526E+04 0.5115763E+01