#modelname1 = "env.fs100-ld-sev_dlog.f-fir.2x15-lvl.1-dexp.1_init-mt.shr-basic"
modelname1 = "env.fs100-ld-sev_dlog.f-fir.2x15-lvl.1-dexp.1_init-basic"
modelname2 = "env.fs100-ld-sev_dlog.f-wc.2x3.c-stp.3-fir.3x15-lvl.1-dexp.1_init-basic"
#modelname2 = "env.fs100-ld-sev_dlog.f-wc.2x3.c.n-stp.3-fir.3x15-lvl.1-dexp.1_init-basic"
#modelname2 = "env.fs100-ld-sev_dlog.f-wc.2x2.c.n-stp.2-fir.2x15-lvl.1-dexp.1_init-basic"
save_figs = False
outpath = "/auto/users/svd/docs/current/two_band_spn/eps_rev/"
#if save_figs:
plt.close('all')
#cellid="por077a-c1"
#cellid = "chn003c-a1"
cellid = "eno009d-a1"
cellid = "eno027d-c1"
cellid = "eno029d-c1"
cellid = "por074b-d2"
cellid = "por074b-c2"
cellid = "por020a-c1"
fh1, ctx1, ctx2 = lplt.compare_model_preds(cellid, batch, modelname1,
modelname2)
# xf2, ctx2 = lplt.get_model_preds(cellid, batch, modelname2)
fh2 = nplt.diagnostic(ctx2, pre_dur=0.25, dur=1.5)
if save_figs:
print("saving figures...")
fh1.savefig(outpath + "fig1.detailed_response_" + cellid + ".pdf")
fh2.savefig(outpath + "fig4.model_steps_" + cellid + ".pdf")
cellid = "TAR010c-24-2"
cellid = "TAR010c-38-2"
cellid = "AMT003c-33-1"
cellid = "bbl104h-12-1"
xf1, ctx1 = load_model_xform(batch=batch,
modelname=modelnames[0],
cellid=cellid)
xf2, ctx2 = load_model_xform(batch=batch,
modelname=modelnames[1],
cellid=cellid)
fh, ctx1, ctx2 = lplt.compare_model_preds(cellid,
batch,
modelnames[0],
modelnames[1],
max_pre=0.5,
max_dur=4,
stim_ids=[0, 1],
ctx1=ctx1,
ctx2=ctx2)
#xf,ctx1=load_model_xform(batch=batch, modelname=modelnames[0], cellid=cellid)
#xf,ctx2=load_model_xform(batch=batch, modelname=modelnames[1], cellid=cellid)
#fh1=nplt.quickplot(ctx1)
#fh2=nplt.quickplot(ctx2)
if save_fig:
fh.savefig(outpath + fileprefix + "_example_" + cellid + "_" +
"modelcomp_batch_" + str(batch) + ".pdf")