%pylab
Using matplotlib backend: <object object at 0x7f57101c0a10> Populating the interactive namespace from numpy and matplotlib
%matplotlib inline
import h5py
f220 = h5py.File('./noncirc-marpol//220/result.hdf','r')
f330 = h5py.File('./noncirc-marpol//330/result.hdf','r')
f221 = h5py.File('./noncirc-marpol//221/result.hdf','r')
bf220 = f220.attrs['log_evidence']
bf330 = f330.attrs['log_evidence']
bf221 = f221.attrs['log_evidence']
bf330 - max(bf221,bf220)
3.7490169680677354
np.exp(bf330 - max(bf221,bf220))
42.47930295709339
cir_f220 = h5py.File('../github/posteriors/kerr/220/KERR-220-06MS.hdf','r')
cir_f330 = h5py.File('../github/posteriors/kerr/220_330/KERR-220_330-06MS.hdf','r')
cir_f221 = h5py.File('../github/posteriors/kerr/220_221/KERR-220_221-06MS.hdf','r')
cir_bf220 = cir_f220.attrs['log_evidence']
cir_bf330 = cir_f330.attrs['log_evidence']
cir_bf221 = cir_f221.attrs['log_evidence']
cir_bf330 - max(cir_bf221,cir_bf220)
4.00806698994711
np.exp(cir_bf330 - max(cir_bf221,cir_bf220))
55.040374070835256
cir_bf220 - bf220
-0.24397759791463614
np.exp(cir_bf220 - bf220)
0.7835051862163263
cir_bf330 - bf330
0.015072423964738846
np.exp(cir_bf330-bf330)
1.015186585790832
cir_bf221 - bf221
0.10032191732898355
np.exp(cir_bf221 - bf221)
1.105526749016665
cir_bf330 - max(bf220,cir_bf221)
3.7640893920324743
np.exp(cir_bf330 - max(bf220,cir_bf221))
43.12441853578603
plt.bar(1,np.exp(cir_bf220-bf220),label='220: cir vs noncir')
plt.bar(2,np.exp(cir_bf221-bf221),label='221: cir vs noncir')
plt.bar(3,np.exp(cir_bf330-bf330),label='330: cir vs noncir')
plt.title('Comparison of Bays factor')
plt.xticks([])
plt.legend()
<matplotlib.legend.Legend at 0x7f568cecba58>
plt.bar(1,np.exp(bf221-bf220),label='Noncir 221 vs 220')
plt.text(0.9,5,'%.2f' % np.exp(bf221-bf220))
plt.bar(2,np.exp(bf330-bf220),label='Noncir 330 vs 220')
plt.text(1.9,43,'%.2f' % np.exp(bf330-bf220))
plt.title('Comparison of Bays factor')
plt.xticks([])
plt.legend()
<matplotlib.legend.Legend at 0x7f568d5a9908>