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/')