benda/nonlinearbaseline2025
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

improved results

Browse Source
This commit is contained in:
Jan Benda
2025-05-27 09:59:33 +02:00
parent fcf9cfa431
commit 742fa901c7
4 changed files with 26 additions and 16 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
dataoverview.py
View File

@@ -141,8 +141,6 @@ def si_stats(title, data, sicol, si_thresh, nsegscol):
print(f' high SI cells: n={len(hcells):3d}, {100*len(hcells)/ncells:4.1f}%')
print(f' high SI recordings: n={np.sum(sidata > si_thresh):3d}, '
f'{100*np.sum(sidata > si_thresh)/nrecs:4.1f}%')
nsegs = data[nsegscol]
print(f' number of segments: {np.min(nsegs):4.0f} - {np.max(nsegs):4.0f}, median={np.median(nsegs):4.0f}, mean={np.mean(nsegs):4.0f}, std={np.std(nsegs):4.0f}')
nrecs = []
for cell in cells:
nrecs.append(len(data[data["cell"] == cell, :]))
@@ -154,6 +152,16 @@ def si_stats(title, data, sicol, si_thresh, nsegscol):
contrasts = 100*data['contrast']
print(' contrasts: ', ' '.join([f'{c:.2g}%' for c in np.unique(contrasts)]))
print(f' contrasts: {np.min(contrasts):.2g}% - {np.max(contrasts):.2g}%, median={np.median(contrasts):.2g}%, mean={np.mean(contrasts):.2g}%, std={np.std(contrasts):.2g}%')
nsegs = data[nsegscol]
print(f' number of segments: {np.min(nsegs):4.0f} - {np.max(nsegs):4.0f}, median={np.median(nsegs):4.0f}, mean={np.mean(nsegs):4.0f}, std={np.std(nsegs):4.0f}')
nsegs = data['nsegs']
print(f' available segments: {np.min(nsegs):4.0f} - {np.max(nsegs):4.0f}, median={np.median(nsegs):4.0f}, mean={np.mean(nsegs):4.0f}, std={np.std(nsegs):4.0f}')
trials = data['trials']
print(f' trials: {np.min(trials):.0f} - {np.max(trials):.0f}, median={np.median(trials):.0f}, mean={np.mean(trials):.0f}, std={np.std(trials):.0f}')
duration = data['duration']
print(f' duration: {np.min(duration):.1f}s - {np.max(duration):.1f}s, median={np.median(duration):.1f}s, mean={np.mean(duration):.1f}s, std={np.std(duration):.1f}s')
duration *= trials
print(f' total duration: {np.min(duration):.1f}s - {np.max(duration):.1f}s, median={np.median(duration):.1f}s, mean={np.mean(duration):.1f}s, std={np.std(duration):.1f}s')
cols = ['cvbase', 'respmod2', 'ratebase', 'vsbase', 'serialcorr1', 'burstfrac', 'ratestim', 'cvstim']
for i in range(len(cols)):
for j in range(i + 1, len(cols)):