In [1]:
%pylab
Using matplotlib backend: <object object at 0x7f57101c0a10>
Populating the interactive namespace from numpy and matplotlib
In [20]:
%matplotlib inline
In [2]:
import h5py

Non cir pol¶

In [3]:
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')
In [4]:
bf220 = f220.attrs['log_evidence']
bf330 = f330.attrs['log_evidence']
bf221 = f221.attrs['log_evidence']
In [5]:
bf330 - max(bf221,bf220)
Out[5]:
3.7490169680677354
In [6]:
np.exp(bf330 - max(bf221,bf220))
Out[6]:
42.47930295709339

Circular polarization runs¶

In [7]:
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')
In [8]:
cir_bf220 = cir_f220.attrs['log_evidence']
cir_bf330 = cir_f330.attrs['log_evidence']
cir_bf221 = cir_f221.attrs['log_evidence']
In [9]:
cir_bf330 - max(cir_bf221,cir_bf220)
Out[9]:
4.00806698994711
In [10]:
np.exp(cir_bf330 - max(cir_bf221,cir_bf220))
Out[10]:
55.040374070835256

Circular pol vs Noncircular pol¶

In [11]:
cir_bf220 - bf220
Out[11]:
-0.24397759791463614
In [12]:
np.exp(cir_bf220 - bf220)
Out[12]:
0.7835051862163263
In [13]:
cir_bf330 - bf330
Out[13]:
0.015072423964738846
In [14]:
np.exp(cir_bf330-bf330)
Out[14]:
1.015186585790832

221¶

In [15]:
cir_bf221 - bf221
Out[15]:
0.10032191732898355
In [16]:
np.exp(cir_bf221 - bf221)
Out[16]:
1.105526749016665

¶

In [17]:
cir_bf330 - max(bf220,cir_bf221)
Out[17]:
3.7640893920324743
In [18]:
np.exp(cir_bf330 - max(bf220,cir_bf221))
Out[18]:
43.12441853578603

Plot¶

In [35]:
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()
Out[35]:
<matplotlib.legend.Legend at 0x7f568cecba58>
In [50]:
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()
Out[50]:
<matplotlib.legend.Legend at 0x7f568d5a9908>
In [ ]: