def gtCompare(gtSpikes,data,bitmasks): avg,avgB,u = [],[],0 for spk in gtSpikes: for i,mask in enumerate(bitmasks): if u == 0: avg.append(0.0) avgB.append(np.median(roip.find_avg_over_tc(data,mask))) avg[i] += (sum(roip.find_avg_over_tc(data[spk:spk+10],mask))/10.0)-avgB[i] u = 1 return avg
def corCompare(fluor,data,bitmasks): r = [] for mask in bitmasks: tc = roip.find_avg_over_tc(data,mask) r.append(pearsonr(fluor,tc)) return r #def deleteMultiDetection_backup(bitmasks,data): """back up of a different way of dealing with multiple detections of the same physical ROI across z.""" """