timothydmorton / isochrones

Pythonic stellar model grid access; easy MCMC fitting of stellar properties
http://isochrones.readthedocs.org
MIT License
120 stars 62 forks source link

No samples being generated when adjusting MCMC parameters #72

Open atw5179 opened 6 years ago

atw5179 commented 6 years ago

When attempting to use the MCMC fit with parameters other than those set by default, no samples are generated and this error occurs when the corner plot function is called:

Traceback (most recent call last): File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/isochrones/starmodel.py", line 980, in corner fig = corner.corner(df[params], labels=params, priors=priors, **kwargs) File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/corner/corner.py", line 147, in corner assert xs.shape[0] <= xs.shape[1], "I don't believe that you want more " \ AssertionError: I don't believe that you want more dimensions than samples!

During handling of the above exception, another exception occurred:

Traceback (most recent call last): File "isochrone_distance_models_no_gaia_mcmc.py", line 45, in mod_b.corner_plots(basename='/gpfs/group/ebf11/default/atw5179/distance_plots_no_gaia_mist_mcmc/'+i+'_bound_no_gaia_corner'); File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/isochrones/starmodel.py", line 991, in corner_plots fig1, fig2 = self.corner_physical(kwargs), self.corner_observed(kwargs) File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/isochrones/starmodel.py", line 1010, in corner_physical return self.corner(props, range=rng, kwargs) File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/isochrones/starmodel.py", line 983, in corner fig = corner.corner(df[params], labels=params, kwargs) File "/gpfs/group/ebf11/default/anaconda3/lib/python3.6/site-packages/corner/corner.py", line 147, in corner assert xs.shape[0] <= xs.shape[1], "I don't believe that you want more " \ AssertionError: I don't believe that you want more dimensions than samples!

For this particular run, I believe nwalkers = 1000, niter = 10000, and nburn = 5000. However, I have also used lower and higher values of the number of walkers, iterations, and burn-ins. The fits were attempted using star_ini files like the following:

maxAV = 0.162 RA = 287.698000 dec = 42.338718 Teff = 5507.7, 60 feh = 0.108, 0.04 logg = 4.469, 0.1

[twomass] J = 11.252, 0.021 H = 10.910, 0.019 K = 10.871, 0.013

[Lick] resolution = 0.5 separation_1 = 3.839 PA_1 = 53.229 H_1 = 4.343, 0.1 K_1 = 4.105, 0.1

Which work fine when using Multinest.