# 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,\