import matplotlib.pyplot as plt
from javelin.zylc import get_data
from javelin.lcmodel import Cont_Model, Rmap_Model
javdata_con = get_data([
"../data/con_test.txt",
], names=[
"con",
])
javdata_rm = get_data(["../data/con_test.txt", "../data/line_test.txt"],
names=["con", "line"])
cont = Cont_Model(javdata_con)
#cont.do_mcmc(fchain="mychain0.dat",nwalkers=100, nburn=200, nchain=500)
cont.load_chain("mychain0.dat")
cont.show_hist(bins=100)
cont.get_hpd()
conthpd = cont.hpd
print(conthpd[1, :])
cont_best = cont.do_pred(conthpd[1, :])
x = cont_best.jlist[:]
y = cont_best.mlist[:] + cont_best.blist[:]
ey = cont_best.elist[:]
plt.plot(x[0], y[0] + javdata_con.cont_mean)
plt.fill_between(x[0],
y[0] - ey[0] + javdata_con.cont_mean,
y[0] + ey[0] + javdata_con.cont_mean,
where=None,
color='grey')
x = javdata_con.jlist[:]
y = javdata_con.mlist[:] + javdata_con.blist[:]
from javelin.zylc import get_data
from javelin.lcmodel import Cont_Model
from javelin.lcmodel import Rmap_Model
# for i in xrange(50):
for i in xrange(1):
c = get_data(["dat/continuum.dat"])
cmod = Cont_Model(c)
cmod.do_mcmc(nwalkers=200, nburn=100, nchain=200, threads=1)
a = cmod.get_hpd()
javdata4 = get_data(["dat/continuum.dat", "dat/yelm.dat", "dat/zing.dat"],
names=["Continuum", "Yelm", "Zing"])
rmap2 = Rmap_Model(javdata4)
rmap2.do_mcmc(conthpd=a, nwalkers=200, nburn=100, nchain=200, threads=1)
del rmap2
del c
print i
