In [ ]:

    

%matplotlib inline
import numpy as np
import sys
import corner


    

In [ ]:

    

data = np.loadtxt("/home/desvignes/projets/1906/new_data/results/omega-fixed/chains-post_equal_weights.dat", usecols=(0,1,2,3,17))


    

In [ ]:

    

figure = corner.corner(data, labels=[r"$\alpha$", r"$\zeta$", r"$\phi_0$", r"$\psi_0$", r"$i$"], 
                         bins=64,
                         quantiles=[0.16, 0.5, 0.84],
                         #quantiles=[0.023, 0.5, 0.977],
                         show_titles=True, title_args={"fontsize": 12},
                      plot_datapoints=False, levels=1.0 - np.exp(-0.5 * np.arange(1, 3.1, 1.) ** 2)
                      )


    

In [ ]:

    

#figure.savefig('1257-RVM.pdf', format='pdf')


    

In [ ]: