Error in plotting quality metrics of data sorted with kilosort4 about bombcell

[JohannaMankel]

Hi Julie, this works perfectly! Thank you so much,
Johanna

[Julie-Fabre]

Hi Johanna,

Thanks for your message! Bombcell hasn't been fully tested with Kilosort4 yet. Would you be able to send me your qMetric and param structures that bombcell computes (the inputs to bc_plotGlobalQualityMetric) so I can troubleshoot?

[JohannaMankel]

Hi Julie, sure. Here are the files:
param_Johanna.csv
qMetric_Johanna.csv

[Julie-Fabre]

Hi Johanna,

Thank you very much!

I had a quick look at the error was because in this dataset no units were being classified as Good or MUA (so there was nothing to plot on the second upSet plot, and that caused on error) - you can see in the first upSet plot (below)
. This is now fixed - bombcell will display a text warning you there are no units in a catagery but won't error anymore.

This classification was due to most of your units being classified as noise with a waveform baseline flatness above the threshold, and this is because Kilosort4 extracts templates differently than previous Kilosort versions. It used to be 82 samples and it is now 61, which means the baseline window we define with param.waveformBaselineWindowStart, param.waveformBaselineWindowStop awas now wrong, and usually included the units' main trough/peak. I have now added an input kilosortVersion in bc_qualityParamValues to get values that make more sense for kilosort4.

Let me know if this works better for you!