In [2]:
import numpy as np
import pandas as pd
import emcee
In [5]:
DF = pd.read_table('../grid/SG_US_step.dat', sep='\t')
In [6]:
DF.head()
Out[6]:
id
M
Y
Z
alpha
diffusion
settling
eta
overshoot
undershoot
...
r13_58
r13_59
r10_57
r10_58
r02_61
r01_59
r01_60
r13_60
r10_59
r02_60
0
20000
0.75408
0.258478
0.000279
2.220276
0.009972
0.009972
0
0.35059
0.083513
...
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
1
20000
0.75408
0.258478
0.000279
2.220276
0.009972
0.009972
0
0.35059
0.083513
...
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
2
20000
0.75408
0.258478
0.000279
2.220276
0.009972
0.009972
0
0.35059
0.083513
...
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
3
20000
0.75408
0.258478
0.000279
2.220276
0.009972
0.009972
0
0.35059
0.083513
...
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
4
20000
0.75408
0.258478
0.000279
2.220276
0.009972
0.009972
0
0.35059
0.083513
...
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
NaN
5 rows × 348 columns
In [22]:
DF = DF[np.logical_or.reduce((
DF['diffusion'] > 0.8,
DF['diffusion'] < 1.2,
DF['settling'] > 0.8,
DF['settling'] < 1.2,
DF['overshoot'] < 0.2,
DF['undershoot'] < 0.2))]
In [70]:
DF = DF.drop(DF.index[np.isnan(DF['dnu02'])])
In [7]:
KIC = pd.read_table('../regression/perturb/final/7970740_perturb.dat')
In [27]:
KIC.head()
Out[27]:
Teff
Fe/H
nu_max
Dnu0
dnu02
r02
r01
r10
0
5309.000000
-0.540000
4197.000000
173.379034
8.536836
0.049290
0.014690
0.014803
1
5406.247480
-0.572623
4259.500565
173.461297
8.348929
0.048154
0.014845
0.014886
2
5291.495691
-0.513386
4179.294972
173.189313
8.578192
0.049516
0.014402
0.014621
3
5326.634508
-0.440346
4233.747382
173.299814
8.426018
0.048641
0.014673
0.014663
4
5274.490584
-0.492826
4141.516938
173.399443
8.282791
0.047827
0.014722
0.015179
In [41]:
Teff = KIC['Teff'][0]
dTeff = np.std(KIC['Teff'])
FeH = KIC['Fe/H'][0]
dFeH = np.std(KIC['Fe/H'])
nu_max = KIC['nu_max'][0]
dnu_max = np.std(KIC['nu_max'])
Dnu0 = KIC['Dnu0'][0]
dDnu0 = KIC['Dnu0'][0]
dnu02 = KIC['dnu02'][0]
ddnu02 = np.std(KIC['dnu02'])
In [91]:
min_Y = np.min(DF['Y'])
max_Y = np.max(DF['Y'])
min_Z = np.min(DF['Z'])
max_Z = np.max(DF['Z'])
min_a = np.min(DF['alpha'])
max_a = np.max(DF['alpha'])
def lnprior(theta):
M, Y, Z, alpha, age = theta
if 0.7 < M < 0.9 and min_Y < Y < max_Y and \
min_Z < Z < max_Z and min_a < alpha < max_a and \
8 < age < 13:
return 0.0
return -np.inf
def lnlike(theta):
M, Y, Z, alpha, age = theta
idx = (np.abs(DF['M']-M) + \
np.abs(DF['Y']-Y) + \
np.abs(DF['Z']-Z) + \
np.abs(DF['alpha']-alpha) + \
np.abs(DF['age']-age)).values.argmin()
row = DF.iloc[idx]
chi2 = (row['Teff'] - Teff)**2 / dTeff + \
(row['Fe_H'] - FeH)**2 / dFeH + \
(row['nu_max'] - nu_max)**2 / dnu_max + \
(row['Dnu0'] - Dnu0)**2 / dDnu0 + \
(row['dnu02'] - dnu02)**2 / ddnu02
return(-chi2/2)
def lnprob(theta):
lp = lnprior(theta)
if not np.isfinite(lp):
return -np.inf
return lp + lnlike(theta)
In [93]:
lnprob([0.75408, 0.258478, 0.000279, 2.220276, 9])
#lnprob([7.21748641e-01, 2.47988905e-01, 3.90556787e-03, 1.54462387e+00,
# 1.07337299e+01])
#lnprob([0.82510679, 0.23599072, 0.01876741, 1.52110795,10.73980817])
#M, Y, Z, alpha, age = [7.21748641e-01, 2.47988905e-01, 3.90556787e-03, 1.54462387e+00,
# 1.07337299e+01]
#idx = (np.abs(DF['M']-M) + \
# np.abs(DF['Y']-Y) + \
# np.abs(DF['Z']-Z) + \
# np.abs(DF['alpha']-alpha) + \
# np.abs(DF['age']-age)).values.argmin()
#print(idx)
#print(DF.shape)
#row = DF.iloc[idx]
#row
Out[93]:
-14768.421605072388
In [94]:
ndim, nwalkers = 5, 100
pos = [[7.601373789062498831e-01,
2.480115351562501291e-01,
7.114554687500006941e-03,
1.557601917968750138e+00,
1.070942461134874435e+01] + 1e-4*np.random.randn(ndim)
for i in range(nwalkers)]
In [95]:
sampler = emcee.EnsembleSampler(nwalkers, ndim, lnprob)
sampler.run_mcmc(pos, 500)
Out[95]:
(array([[7.13090434e-01, 2.66569250e-01, 1.43809643e-02, 1.37375424e+00,
1.07025513e+01],
[7.11343815e-01, 2.33124751e-01, 3.15593137e-02, 1.38310917e+00,
1.06742072e+01],
[8.36411909e-01, 2.78383652e-01, 2.64391857e-02, 1.56469230e+00,
1.07279507e+01],
[7.19065176e-01, 2.55899131e-01, 8.56046024e-03, 1.39378670e+00,
1.07092428e+01],
[7.29060450e-01, 2.26846355e-01, 1.43732221e-02, 1.58331367e+00,
1.06264320e+01],
[7.48392728e-01, 3.20868179e-01, 3.86751454e-02, 1.51110822e+00,
1.07291327e+01],
[8.92390236e-01, 2.43017880e-01, 3.03166153e-02, 1.72783563e+00,
1.07718002e+01],
[7.01807776e-01, 2.36327028e-01, 1.31985288e-02, 1.41302581e+00,
1.07107101e+01],
[7.02166975e-01, 2.33383770e-01, 8.05570317e-03, 1.78034166e+00,
1.09293734e+01],
[7.08311921e-01, 2.21620966e-01, 2.67323905e-02, 1.34497127e+00,
1.06457416e+01],
[7.95879481e-01, 2.79233145e-01, 2.80060182e-02, 1.55389113e+00,
1.06725538e+01],
[7.90769087e-01, 2.38639116e-01, 9.08707447e-03, 1.65306003e+00,
1.06079756e+01],
[7.23985527e-01, 3.18509085e-01, 1.06280969e-02, 1.76386887e+00,
1.09322156e+01],
[7.14701170e-01, 2.46054622e-01, 1.20738565e-02, 1.37779891e+00,
1.06536988e+01],
[7.12779882e-01, 2.67333195e-01, 1.85693687e-02, 1.36539166e+00,
1.06645439e+01],
[8.06053683e-01, 3.13153791e-01, 1.43328317e-02, 1.59299249e+00,
1.07165962e+01],
[7.02309676e-01, 2.73232279e-01, 9.73357185e-03, 1.40661509e+00,
1.06708732e+01],
[8.27881979e-01, 3.26340232e-01, 1.88002804e-02, 1.89066752e+00,
1.07676050e+01],
[7.04950151e-01, 2.30832219e-01, 2.39102413e-02, 1.37884921e+00,
1.06629831e+01],
[7.06926609e-01, 2.20559707e-01, 2.43378629e-02, 1.39818678e+00,
1.06831578e+01],
[7.12172464e-01, 2.61494921e-01, 7.89167040e-03, 1.40235225e+00,
1.07030726e+01],
[7.02108591e-01, 2.53644101e-01, 9.18229498e-03, 1.33722217e+00,
1.06930989e+01],
[7.09370808e-01, 2.74496524e-01, 9.88813759e-03, 1.37827687e+00,
1.06599820e+01],
[7.02533117e-01, 2.49433918e-01, 1.46638146e-02, 1.38226588e+00,
1.06587177e+01],
[7.35904897e-01, 2.81524363e-01, 1.48406283e-02, 1.75149778e+00,
1.09271010e+01],
[7.06738748e-01, 2.33183443e-01, 3.95745174e-03, 1.78236221e+00,
1.09491055e+01],
[7.63749760e-01, 2.55394032e-01, 1.40939040e-02, 1.43034654e+00,
1.06629877e+01],
[7.91698687e-01, 2.36729232e-01, 3.14647256e-03, 1.65760993e+00,
1.04520742e+01],
[8.20921923e-01, 3.10234978e-01, 1.60700944e-02, 1.54476241e+00,
1.07014839e+01],
[7.67909093e-01, 2.83421240e-01, 1.68441088e-02, 1.23389005e+00,
1.06001058e+01],
[7.11231690e-01, 2.20896596e-01, 1.38610183e-02, 1.36712705e+00,
1.06902291e+01],
[7.37170434e-01, 3.10877529e-01, 3.91192274e-02, 1.51365459e+00,
1.07302040e+01],
[7.13062683e-01, 2.55889602e-01, 1.35958607e-02, 1.38426870e+00,
1.07233856e+01],
[7.33888895e-01, 2.59705879e-01, 1.33314810e-02, 1.59390116e+00,
1.06218766e+01],
[8.24166458e-01, 2.48432748e-01, 2.49684444e-02, 1.78011174e+00,
1.07759832e+01],
[7.49630085e-01, 2.24564467e-01, 3.19551641e-03, 1.41237509e+00,
1.07005965e+01],
[8.13889200e-01, 2.69045646e-01, 1.35713715e-02, 1.81226349e+00,
1.08197789e+01],
[8.49416250e-01, 2.27722499e-01, 3.94201701e-02, 2.15221772e+00,
1.09371707e+01],
[8.32058617e-01, 2.76048285e-01, 2.31297627e-02, 1.75755824e+00,
1.07931813e+01],
[7.68007790e-01, 3.07584015e-01, 2.28990105e-02, 1.50229068e+00,
1.07208770e+01],
[7.17337594e-01, 2.40101879e-01, 2.17426617e-02, 1.41626604e+00,
1.06916917e+01],
[7.60042390e-01, 3.04773221e-01, 3.67834747e-03, 1.50701162e+00,
1.07152381e+01],
[7.65472167e-01, 2.67804086e-01, 1.62136942e-03, 1.26122542e+00,
1.05971190e+01],
[7.18290199e-01, 2.40394993e-01, 4.86086930e-03, 1.04916640e+00,
1.06162739e+01],
[8.16579750e-01, 3.08425036e-01, 2.16943726e-02, 1.83123802e+00,
1.07817317e+01],
[7.75117276e-01, 2.68309862e-01, 1.02979527e-02, 1.27990443e+00,
1.05877157e+01],
[7.11276629e-01, 2.32219722e-01, 5.97321162e-03, 1.02002372e+00,
1.05924594e+01],
[7.08030811e-01, 2.25107916e-01, 3.13307134e-03, 1.41047598e+00,
1.06841101e+01],
[7.63647444e-01, 2.30210150e-01, 1.99708603e-02, 1.65650625e+00,
1.06317582e+01],
[7.20723329e-01, 2.24056466e-01, 1.59342066e-02, 1.40068123e+00,
1.06925368e+01],
[7.13214470e-01, 2.23365406e-01, 6.72719273e-03, 1.38025469e+00,
1.06840653e+01],
[7.43952724e-01, 2.86186015e-01, 3.89794480e-02, 1.89522655e+00,
1.08462306e+01],
[7.08443874e-01, 2.62610123e-01, 3.90371933e-02, 1.40403108e+00,
1.06936775e+01],
[7.14189089e-01, 2.20275075e-01, 5.33745121e-03, 1.36275784e+00,
1.06822957e+01],
[7.10322568e-01, 2.38103135e-01, 3.44835470e-02, 1.35917914e+00,
1.06776586e+01],
[7.06478519e-01, 2.44032665e-01, 1.12019228e-02, 1.36860080e+00,
1.06966411e+01],
[7.08730380e-01, 2.27804059e-01, 2.71656056e-02, 1.38838866e+00,
1.06858700e+01],
[7.13667761e-01, 2.44974706e-01, 9.58714115e-03, 1.37208144e+00,
1.07121457e+01],
[7.37356194e-01, 2.49710383e-01, 2.67617189e-02, 1.64272012e+00,
1.06290659e+01],
[7.10829572e-01, 3.13717158e-01, 3.11643056e-02, 1.52553706e+00,
1.07073593e+01],
[7.18812641e-01, 2.41798463e-01, 3.54726964e-02, 1.62570218e+00,
1.05630845e+01],
[7.17156510e-01, 2.38885098e-01, 1.85226157e-02, 1.39409567e+00,
1.06735085e+01],
[7.00919980e-01, 2.68593193e-01, 1.55747766e-02, 1.43415074e+00,
1.06904106e+01],
[7.03588821e-01, 2.45531218e-01, 1.94910034e-02, 1.40405894e+00,
1.06929410e+01],
[7.01446131e-01, 2.30036382e-01, 2.12300772e-02, 1.37785756e+00,
1.07243261e+01],
[7.07539873e-01, 2.59411805e-01, 1.12012077e-02, 1.37843834e+00,
1.06232658e+01],
[8.63310685e-01, 2.88719990e-01, 7.47233188e-03, 1.56699192e+00,
1.07126785e+01],
[7.98402627e-01, 2.94747284e-01, 3.72332340e-02, 1.91878278e+00,
1.08133855e+01],
[8.83847929e-01, 3.13553635e-01, 2.87288279e-02, 1.72528836e+00,
1.07461980e+01],
[7.07549645e-01, 2.22126580e-01, 6.69698278e-03, 1.37824868e+00,
1.07234607e+01],
[8.01622264e-01, 3.16799989e-01, 3.06169155e-02, 1.52889605e+00,
1.07248024e+01],
[7.45606241e-01, 2.28971152e-01, 1.92114238e-03, 1.42553214e+00,
1.07071073e+01],
[7.53031787e-01, 2.85652539e-01, 2.19840109e-02, 1.86357738e+00,
1.08181501e+01],
[7.04822341e-01, 2.54762251e-01, 3.57686342e-02, 1.43269995e+00,
1.06831304e+01],
[7.20576406e-01, 3.05676835e-01, 2.25118251e-03, 1.53802546e+00,
1.07369217e+01],
[8.63088595e-01, 2.67575045e-01, 2.72312665e-02, 1.68326832e+00,
1.07480962e+01],
[7.62385252e-01, 3.39693886e-01, 3.43136085e-02, 1.69009036e+00,
1.06474340e+01],
[8.04038217e-01, 2.55726170e-01, 1.21731415e-02, 1.39391194e+00,
1.06645210e+01],
[7.60260211e-01, 3.09564691e-01, 2.94151844e-02, 1.55871474e+00,
1.07124168e+01],
[7.06351083e-01, 2.62487439e-01, 3.30843197e-02, 1.38219377e+00,
1.07018887e+01],
[7.01004026e-01, 2.69797016e-01, 3.12793205e-02, 1.38809840e+00,
1.06891613e+01],
[7.09365790e-01, 2.49260593e-01, 1.83172322e-02, 1.36054350e+00,
1.06271136e+01],
[7.18187260e-01, 2.39323768e-01, 5.98102526e-03, 1.41855044e+00,
1.06971454e+01],
[7.25892674e-01, 2.86625147e-01, 2.18201152e-02, 1.73977074e+00,
1.09618811e+01],
[7.14702738e-01, 2.38000470e-01, 2.71586632e-02, 1.35239516e+00,
1.06913391e+01],
[7.00794406e-01, 2.33276519e-01, 1.18765613e-02, 1.39291663e+00,
1.06753310e+01],
[7.32414775e-01, 3.31996463e-01, 2.39304287e-02, 1.70017654e+00,
1.06789232e+01],
[7.27981949e-01, 2.90870053e-01, 1.76468264e-03, 1.62241015e+00,
1.09973052e+01],
[8.83866349e-01, 2.94741506e-01, 1.05887021e-02, 1.73960075e+00,
1.07619093e+01],
[7.06857740e-01, 2.55773647e-01, 2.73363113e-02, 1.35891973e+00,
1.06368971e+01],
[7.21308079e-01, 2.98591921e-01, 5.54310663e-03, 1.67161178e+00,
1.10807051e+01],
[7.08303287e-01, 2.66719312e-01, 3.24843309e-02, 1.43130027e+00,
1.06787966e+01],
[7.12138359e-01, 2.26732349e-01, 6.58651021e-03, 1.39872907e+00,
1.07186675e+01],
[7.19091409e-01, 2.46483820e-01, 9.50243177e-03, 1.40925745e+00,
1.06590226e+01],
[7.50552921e-01, 3.07653250e-01, 3.34766223e-02, 1.52732369e+00,
1.07226943e+01],
[7.87686178e-01, 2.66408468e-01, 9.11301989e-03, 1.23232306e+00,
1.06023671e+01],
[7.29714148e-01, 2.65968055e-01, 2.69841951e-02, 1.40901893e+00,
1.06794295e+01],
[7.76585929e-01, 2.26803978e-01, 1.02068315e-02, 1.62348919e+00,
1.04363545e+01],
[7.41271283e-01, 3.21765290e-01, 1.34083972e-02, 1.54726637e+00,
1.07055109e+01],
[7.34599236e-01, 2.43467147e-01, 3.16197930e-03, 1.89178432e+00,
1.08369085e+01]]),
array([ -170.84838319, -170.84838319, -1756.55801553, -170.84838319,
-452.2599463 , -1311.52728189, -701.99040092, -170.84838319,
-1419.68527819, -170.84838319, -1756.55801553, -452.2599463 ,
-274.2649876 , -170.84838319, -170.84838319, -1756.55801553,
-170.84838319, -371.47217196, -170.84838319, -170.84838319,
-170.84838319, -170.84838319, -170.84838319, -170.84838319,
-274.2649876 , -1419.68527819, -887.32983878, -181.62430464,
-1756.55801553, -910.1303167 , -170.84838319, -1311.52728189,
-170.84838319, -452.2599463 , -621.4564915 , -887.32983878,
-763.39529045, -63.50671656, -621.4564915 , -1311.52728189,
-170.84838319, -1311.52728189, -910.1303167 , -163.69949141,
-763.39529045, -910.1303167 , -163.69949141, -170.84838319,
-452.2599463 , -170.84838319, -170.84838319, -327.78531054,
-170.84838319, -170.84838319, -170.84838319, -170.84838319,
-170.84838319, -170.84838319, -452.2599463 , -1311.52728189,
-452.2599463 , -170.84838319, -170.84838319, -170.84838319,
-170.84838319, -170.84838319, -1756.55801553, -500.33526589,
-701.99040092, -170.84838319, -1311.52728189, -887.32983878,
-327.78531054, -170.84838319, -1311.52728189, -701.99040092,
-1142.19337916, -1354.40675988, -1311.52728189, -170.84838319,
-170.84838319, -170.84838319, -170.84838319, -274.2649876 ,
-170.84838319, -170.84838319, -1142.19337916, -956.37840856,
-701.99040092, -170.84838319, -462.11599855, -170.84838319,
-170.84838319, -170.84838319, -1311.52728189, -910.1303167 ,
-170.84838319, -181.62430464, -1311.52728189, -327.78531054]),
('MT19937', array([2354650672, 2141310750, 2937627131, 3887778756, 1942614139,
138983963, 2623666375, 3587133513, 3844697145, 2186184993,
1799864891, 2166308755, 3092160559, 804237428, 1866794912,
2933086344, 3038334293, 1055332369, 1944940624, 2843475482,
493889579, 2809797318, 2281222168, 705894738, 3569740330,
24810295, 2886994166, 2600302973, 1726625092, 4250855207,
4211679265, 453271181, 2727514084, 3832349945, 1728081393,
4103794096, 2063802445, 1523956943, 3724851066, 765035171,
3680180966, 4187632401, 1695799162, 728105803, 2846417261,
1117714600, 418333148, 3016408067, 746328115, 473765952,
2395499343, 4137237708, 3357502939, 4023482257, 2082636578,
3137283699, 1775771684, 2425769787, 2939358796, 2007659448,
631200233, 1568523610, 3569008197, 3007581311, 1079224628,
2816461629, 842970254, 1166360367, 1210309705, 2277805089,
1177142320, 1339124609, 2273841844, 3934938688, 2747515865,
3904385534, 1756239732, 2441794154, 3905968397, 1092796712,
301933310, 1200820315, 1205188095, 3654764603, 1252176823,
1185515485, 4071017349, 334347969, 4196765783, 1651630800,
1065442593, 2642151070, 2298361638, 2062637290, 838873422,
2578317281, 125552211, 443353747, 192670710, 1909695198,
638789967, 4260090017, 2185335617, 2776081428, 2557611553,
3868391071, 3849583452, 496616036, 2447661973, 3866829626,
1943374077, 1901185614, 2991943218, 1829637081, 726345131,
2672005896, 59985607, 3698137977, 3019131028, 954451785,
367314508, 4063030238, 1448632778, 307065777, 1653621494,
2764146993, 2860364473, 906597518, 285850190, 3103472398,
1695820432, 2992678314, 4221383409, 666970654, 2738586734,
1531821831, 3227112886, 4246166225, 3836331500, 2973055772,
3724128235, 3004162876, 2325420803, 2100252627, 1164001416,
515860724, 3686412069, 1632379218, 465072812, 781418325,
2515425561, 1115868069, 3645161681, 2510082276, 3069104873,
2934873800, 1699749134, 4196430282, 78886344, 3341134791,
873611978, 3794740852, 3069118679, 1580861597, 2666430804,
2082087204, 2228369375, 1959834331, 3454656379, 108403112,
169646919, 1915781817, 4281722224, 4169803169, 2921933641,
3419042249, 903861970, 2229097989, 3176094085, 2802988435,
3703459221, 1091161771, 3397307663, 1050300875, 1590564541,
2764663771, 1505564318, 4161689881, 2910836781, 2215015792,
1477617342, 2318036581, 662258524, 3323413123, 1396536234,
4142439468, 1476435870, 2816372204, 2004611802, 1074379391,
3321864069, 420212612, 2491278406, 1588963187, 4202619573,
4166456893, 4151067534, 2601402032, 694053005, 2579149731,
2974426352, 2333119116, 1548145118, 2932436117, 3615938033,
26527279, 1452077869, 4190236145, 2103024852, 3882393467,
2317316805, 111005747, 102676277, 1176454950, 3058394691,
1904422329, 1931665570, 2817245938, 2179060674, 4179109146,
3732194043, 3279814475, 1454637641, 771044230, 308331685,
3719296928, 2732082214, 2792237151, 3831530409, 669777799,
1227841635, 1365352397, 2287779353, 1787422641, 3582074963,
3319872668, 3912834953, 210704227, 1431052044, 1383027151,
3178711739, 1084612969, 4116513395, 1608119997, 1140438384,
2140182991, 3790697439, 19257600, 2067333061, 3692837240,
3363857386, 804872720, 1436221993, 1588951596, 4111606362,
870835700, 3913016709, 209714466, 1275637007, 1552760609,
2583847206, 3891797480, 2341057403, 1266098533, 3091162588,
453747572, 1266228459, 1256044217, 684383282, 851035686,
1362411930, 900883215, 3519625417, 3129374313, 1567096026,
1953579763, 2059632467, 4275098921, 483202983, 2035850546,
567968637, 3545869051, 2726012626, 3118179970, 552562968,
3397017266, 3727050953, 1850293199, 3707024021, 2013942532,
4285636351, 2274937782, 840014701, 1827787305, 2793364101,
4027735186, 624946864, 4261988010, 4174662242, 2442116510,
1884845989, 15374606, 2953849132, 1415754425, 768494729,
2303280392, 2708507401, 3433284621, 2451858676, 7467970,
3136816905, 2794048297, 3006827605, 2057556422, 12107465,
2584670923, 522873654, 1179576010, 2668588222, 936342882,
2940988895, 1725541674, 3430938540, 1111188448, 3055440002,
1523236540, 795418489, 3039984950, 3043038409, 2215405157,
606896927, 652070257, 2916576590, 2100950171, 4256763311,
849786360, 2981029759, 3447171754, 361744720, 358268529,
3854937103, 127697435, 3177601396, 3014984597, 1682272909,
2899157669, 4099393380, 4273664135, 3757869518, 2399532902,
4253009668, 455975597, 2787395230, 1704975193, 2302377729,
832314404, 617456512, 1235856695, 384308388, 4000336356,
3574659824, 3633444627, 41101286, 308900618, 3334515470,
499023274, 1499663345, 1458315604, 122152275, 1336697802,
2663115876, 3266976783, 2908013717, 1928026950, 1227118173,
2532169619, 1818975441, 210883970, 4100684182, 1979203969,
3596822213, 3726299578, 3820498293, 89009277, 293478117,
1442704749, 2168806066, 2991114712, 1981737910, 2795817438,
550759880, 115912720, 2858900255, 3221451841, 1651611851,
1684155626, 1604117545, 2779137317, 1078305801, 2769602676,
391065195, 172651684, 3111811157, 3869287084, 3261046661,
3256017528, 830500001, 3346139553, 2671538985, 3911018301,
4062923565, 3212989812, 1789457690, 3346068426, 2310297444,
3181597358, 1066039956, 981141567, 901225105, 2130734089,
3864510354, 2364572731, 3744148551, 3087731123, 4010395850,
127739479, 3805565510, 1379148986, 4248868154, 2186835460,
3409129123, 3839389460, 235098230, 4211274133, 347105611,
3856415194, 1102209303, 1747047158, 3964636763, 4287242217,
179932101, 2009779217, 961304193, 2853160436, 2506590693,
1656480572, 3506177828, 2863537385, 87159197, 2798497906,
2340356225, 4083024395, 3209778851, 1050962909, 4287846012,
4015077225, 739865597, 954873907, 328079449, 3738919065,
647942246, 620098765, 2831432318, 210760219, 4294760946,
3717358870, 854437415, 3359743252, 3073856523, 1490746675,
1209576080, 941862132, 1099571209, 2564650248, 723239818,
2628276675, 118166931, 1974835034, 3892992033, 670929420,
1235290116, 2577447319, 2351974438, 4177622667, 363445112,
2243618868, 3647069871, 2917647066, 4125231428, 3243916974,
1415501851, 241419995, 3235453602, 1994278041, 3844103697,
1765439693, 2029491646, 1682974475, 2604204134, 3707383297,
4290998805, 374907654, 623680879, 1542076244, 2070598578,
4029913302, 3996324668, 3917812400, 432317798, 3218484448,
308697726, 2422840994, 1970669435, 1776445551, 1255595298,
740584681, 1354746586, 4025594158, 1952629495, 3711828638,
2061284521, 2778309371, 646336854, 3260085641, 3494453847,
1361303196, 4000535316, 3331696895, 1120561261, 2006212216,
2535245246, 1880510924, 597604039, 3291731112, 108967041,
3544154812, 3896910930, 1559345784, 3511193932, 2407531036,
4160133183, 765925613, 2061065800, 2382365510, 1630856348,
3832389507, 1289002313, 1228526345, 2275782516, 1691962348,
3471680587, 2503774260, 1616069078, 2112917413, 1944274249,
3775517156, 3646507747, 65543267, 3037255551, 2554373828,
944672264, 710813750, 3339549401, 1080214693, 1597278508,
1979089837, 2076848381, 3464706227, 381437227, 3757389336,
1054625451, 2480940922, 3153648577, 1946619434, 4075511039,
2132980286, 3264692867, 301992578, 1640959094, 1808330994,
23928160, 538563300, 1935300677, 505323819, 2262004843,
2960932233, 3601154090, 1703543350, 2754627364, 156900840,
2361449198, 1072377668, 3739252943, 808119836, 2481736555,
746236222, 730456161, 2519814155, 3349505741, 2695720888,
19765389, 4115797807, 1722628554, 622751290, 63348976,
4114772656, 4235860086, 4172877107, 948097449, 2650267118,
2338427819, 1631766814, 2672147258, 3345601962], dtype=uint32), 10, 0, 0.0))
In [96]:
samples = sampler.chain[:, 50:, :].reshape((-1, ndim))
In [102]:
%matplotlib inline
import matplotlib
import corner
fig = corner.corner(samples, labels=["$M$", "$Y$", "$Z$", r'$\alpha$', r'$\tau$'])
fig.show()
/scratch/seismo/bellinger/python/anaconda3-5.3.0/lib/python3.7/site-packages/matplotlib/figure.py:457: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
"matplotlib is currently using a non-GUI backend, "
In [100]:
samples[:, 2] = np.exp(samples[:, 2])
list(map(lambda v: (v[1], v[2]-v[1], v[1]-v[0]),
zip(*np.percentile(samples, [16, 50, 84],
axis=0))))
Out[100]:
[(0.7353777174870217, 0.08727497613487079, 0.028190558237808805),
(0.25256167815958314, 0.03790968875300155, 0.020948051206412382),
(15.927722501473882, 0.6985082455182816, 0.5451560921399459),
(1.5400317429711747, 0.07847929804297316, 0.13531837910972477),
(10.705466649181519, 0.045037658095202104, 0.03752895535996714)]
In [ ]:
Content source: earlbellinger/asteroseismology
Similar notebooks: