# ebvgamma = (fit['gamma'][:,1]-fit['gamma'][:,2])[:,None] * fit['k']
# for s,a,b in zip(snname, numpy.median(ebvdelta,axis=0),numpy.median(ebvgamma,axis=0)):
# print s,a,b
snname = numpy.array(data['snlist'])[use]
for i in xrange(len(sivel)):
print '{0} & ${1:5.1f} \pm {2:3.1f}$ & ${3:5.1f} \pm {4:3.1f}$& ${7:5.0f} \pm {8:3.0f}$ & ${5[0]:6.2f} \pm {6[0]:6.2f}$ & ${5[1]:6.2f} \pm {6[1]:6.2f}$& ${5[2]:6.2f} \pm {6[2]:6.2f}$& ${5[3]:6.2f} \pm {6[3]:6.2f}$& ${5[4]:6.2f} \pm {6[4]:6.2f}$ \\\\'.format(
snname[i], EW_obs[i, 0], numpy.sqrt(EW_cov[i, 0, 0]), EW_obs[i, 1],
numpy.sqrt(EW_cov[i, 1, 1]), mag_obs[i, :],
numpy.sqrt(numpy.diagonal(mag_cov[i, :, :])), sivel[i], sivel_err[i])
import json, codecs
EW_mn = EW_obs.mean(axis=0)
sivel_mn = sivel.mean()
def convert(fit, EW_mn, sivel_mn):
rm = ['c_raw','alpha_raw','beta_raw','eta_raw','gamma01','gamma02','gamma03','gamma04','gamma05','k_unit','lp__', \
'Delta_scale','Delta_unit','R_unit','rho11','rho12','rho13','rho14','rho15','EW','sivel','gamma','rho1','k','R', 'mag_int']
use = dict(fit)
use['EWCa'] = use['EW'][:, :, 0] + EW_mn[0]
use['EWSi'] = use['EW'][:, :, 1] + EW_mn[1]
use['lSi'] = use['sivel'] + sivel_mn
use['gamma0'] = use['gamma']
use['gamma1'] = use['rho1']
use['g0'] = use['k']
use['g1'] = use['R']
for r in rm:
sivel_err = sivel_err[use]
EW_obs=EW_obs[use]
mag_obs=mag_obs[use]
EW_cov= EW_cov[use]
mag_cov=mag_cov[use]
nsne, nmags = mag_obs.shape
# # renormalize the data
EW_mn = EW_obs.mean(axis=0)
EW_renorm = (EW_obs - EW_mn)
mag_mn = mag_obs.mean()
mag_renorm = mag_obs-mag_mn
sivel_mn = sivel.mean()
sivel_renorm = sivel-sivel_mn
data = {'D': nsne, 'N_mags': 5, 'N_EWs': 2, 'mag_obs': mag_renorm, 'EW_obs': EW_renorm, 'EW_cov': EW_cov, 'mag_cov':mag_cov, \
'sivel_obs': sivel_renorm, 'sivel_err': sivel_err, 'a':a}
# pystan.misc.stan_rdump(data, 'data.R')
# wefew
Delta_simplex = numpy.zeros(nsne-1)
# Delta_simplex = numpy.zeros(nsne)+1./nsne
# k_simplex = numpy.zeros(nsne)
R_simplex = ((-1.)**numpy.arange(nsne)*.25 + .5)*2./nsne
R_simplex = R_simplex/R_simplex.sum()
init = [{'EW' : EW_renorm, \
'sivel': sivel_renorm,\
