def find_very_good_beta_examples():
minsnr = 7
mincount = 1000
minrate = 9
minqual = 3
preobject = 6,-1000,0 # pre-object
prego = 8,-1000,0 # pre-go
for s,a in sessions_areas():
print s,a
NUNITS = len(metaloadvariable(s,a,'unitIds')[0])
print 'No. Units = ',NUNITS
for u in range(1,1+NUNITS):
if not is_unit_on_good_lfp_channel(s,a,u): continue
if get_unit_quality(s,a,u)<minqual: continue
if get_total_spikes(s,a,u,preobject)<mincount: continue
if get_total_spikes(s,a,u,prego )<mincount: continue
if get_average_rate(s,a,u,preobject)<minrate: continue
if get_average_rate(s,a,u,prego )<minrate: continue
if get_unit_snr(s,a,u)<minsnr: continue
mf1 = isimodefreq(s,a,u,-1000,0000,FS=1000)
if mf1<15 or mf1>30: continue
mf2 = isimodefreq(s,a,u, 2000,3000,FS=1000)
if mf2<15 or mf2>30: continue
printstats(s,a,u)
unitsumfig(s,a,u,savein='./unit_info_figures_exemplar/')
def dumpallunits():
for s,a in sessions_areas():
print s,a
NUNITS = len(metaloadvariable(s,a,'unitIds')[0])
print 'No. Units = ',NUNITS
for i in range(NUNITS):
print '\t',s,a,i
unitsumfig(s,a,i+1,savein='./unit_info_figures/')
for (s,a),us in allunitsbysession.iteritems():
for u in us:
unitsumfig(s,a,u,savein='./unit_info_figures/')
def surveyallunits():
preobject = 6,-1000,0 # pre-object
prego = 8,-1000,0 # pre-go
allrates = []
allcounts = []
allsnr = []
for s,a in sessions_areas():
print s,a
NUNITS = len(metaloadvariable(s,a,'unitIds')[0])
print 'No. Units = ',NUNITS
for u in range(1,1+NUNITS):
printstats(s,a,u)
allrates .append(get_average_rate(s,a,u,preobject))
allrates .append(get_average_rate(s,a,u,prego))
allcounts.append(get_total_spikes(s,a,u,preobject))
allcounts.append(get_total_spikes(s,a,u,prego))
allsnr .append(get_unit_snr(s,a,u))
def filterunits():
minsnr = 5
mincount = 500
minrate = 1
minqual = 2
preobject = 6,-1000,0 # pre-object
prego = 8,-1000,0 # pre-go
for s,a in sessions_areas():
print s,a
NUNITS = len(metaloadvariable(s,a,'unitIds')[0])
print 'No. Units = ',NUNITS
for u in range(1,1+NUNITS):
if not is_unit_on_good_lfp_channel(s,a,u): continue
if get_unit_quality(s,a,u)<minqual: continue
if get_total_spikes(s,a,u,preobject)<mincount: continue
if get_total_spikes(s,a,u,prego )<mincount: continue
if get_average_rate(s,a,u,preobject)<minrate: continue
if get_average_rate(s,a,u,prego )<minrate: continue
if get_unit_snr(s,a,u)<minsnr: continue
printstats(s,a,u)
unitsumfig(s,a,u,savein='./unit_info_figures_good/')
