'BRT034f-02-1', 'BRT034f-07-1', 'BRT034f-34-1', 'BRT036b-06-1',
'BRT036b-07-1', 'BRT036b-10-1', 'BRT036b-12-1', 'BRT036b-20-1',
'BRT036b-28-1', 'BRT036b-31-1', 'BRT036b-36-1', 'BRT036b-38-1',
'BRT036b-40-1', 'BRT036b-45-1', 'BRT036b-46-1', 'BRT037b-03-1',
'BRT037b-13-1', 'BRT037b-30-1', 'BRT037b-33-1', 'BRT037b-36-1',
'BRT037b-63-1', 'TAR010c-01-1', 'TAR010c-06-1', 'TAR010c-12-1',
'TAR010c-15-1', 'TAR010c-18-1', 'TAR010c-19-1', 'TAR010c-21-1',
'TAR010c-27-1', 'TAR010c-27-2', 'TAR010c-30-1', 'TAR010c-33-1',
'TAR010c-58-1'
]
path = '/auto/users/daniela/docs/AC_IC_project/307_per_file/'
for cellid in cellids:
model_per_time_wrapper(cellid,
batch,
basemodel=basemodel,
state_list=state_list)
fig = plt.gcf()
fig.savefig(path + cellid + '.pdf')
plt.close(fig)
state_list_beh = [
'st.pup0.beh0', 'st.pup0.beh', 'st.pup.beh0', 'st.pup.beh'
]
for cellid in cellids:
model_act_pas_DS(cellid,
batch,
basemodel=basemodel,
state_list=state_list_beh)
fig = plt.gcf()
from nems_lbhb.stateplots import model_per_time_wrapper, beta_comp
from nems_lbhb.stateplots import pb_model_plot, ppas_model_plot
state_list = ['st.pup0.hlf0', 'st.pup0.hlf', 'st.pup.hlf0', 'st.pup.hlf']
cellids = ["BRT009c-a1", "ley046f-44-1"] # ["ley046f-41-1", "ley046g-04-1"]
batch = 309
for cellid in cellids:
fh, stats = ppas_model_plot(cellid=cellid,
batch=batch,
loader="psth.fs20.pup",
basemodel="psthfr_stategain.S",
fitter="jk.nf20-basic")
for cellid in cellids:
f = model_per_time_wrapper(cellid,
batch=batch,
loader="psth.fs20.pup-ld-",
fitter="_jk.nf20-basic",
basemodel="-ref-psthfr_stategain.S",
state_list=None,
plot_halves=True)
"""
BRT009c: r2taskonly = 0.60, r2puponly=0.62, r2full=0.62
ley046g: r2taskonly = 0.19, r2puponly=0.22, r2full=0.22
"""