-
Notifications
You must be signed in to change notification settings - Fork 225
Add SLAy auto-merge preset #4190
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: main
Are you sure you want to change the base?
Conversation
| return sig | ||
|
|
||
|
|
||
| def _sliding_RP_viol_pair( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Could we try to unify with this code @chrishalcrow ?
https://github.com/SpikeInterface/spikeinterface/blob/main/src/spikeinterface/qualitymetrics/misc_metrics.py#L1493-L1511
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Probably...
| bin_rate_max = np.max(smoothed_acg) | ||
| max_conts_max = np.array(bTest) / bin_size_ms * 1000 * (bin_rate_max * acceptThresh) | ||
| # compute confidence of less than acceptThresh contamination at each refractory period | ||
| confs = 1 - poisson.cdf(sum_res, max_conts_max) | ||
| rp_viol = 1 - confs.max() |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
One very minor detail!
Add SLAy as a method for
compute_merge_unit_groupshttps://www.biorxiv.org/content/10.1101/2025.06.20.660590v1SLAy combines a template similarity measure, a cross-correlation significance measure and a sliding refractory period violation measure to make a merge decision. You need an analyzer with templates, correlograms and template_similarity. If you do, you can try it using:
or, for more control, e.g.
This is an approximate implementation. Some details
_sliding_RP_violand_compute_xcorrfunctions from https://github.com/saikoukunt/SLAy . I think the sliding RP computation is fairly different between our implementation in quality_metrics. Will take a closer look...To do:
With help and permission from @saikoukunt