Question regarding Unit IDs about bombcell

[Julie-Fabre]

Hi Simon,
Thank you for question and I am glad you're enjoying bombcell!

Bombcell stores the unit ID for each unit in qMetric.clusterID, this corresponds to the kilosort/phy label + 1 (the +1 makes it 1-indexed, and easier to use in MATLAB). I've just pushed a change so it's clearer by adding a field qMetric.phy_clusterID that directly gives you the kilosort/phy label.
The value in qMetric.clusterID value does not necessarily correspond to the row in the qMetric structure the unit is in (ie for unit say 5, it might be in row 3 for instance), because kilosort drops units in the last stages of the algorithm but does not rename them, leading to some empty units.

I pushed changes yesterday to make this clearer. There's now a little blurb and example at the end of bc_qualityMetrics_pipeline, and the GUI now displays all three values: the clusterID, the phy clusterID and the row the unit is in (see image below for what it should look like now)

[Julie-Fabre]

To answer you question about spikes, spike times per unit are extracted with the function bc_loadEphysData. If you use the bc_qualityMetrics_pipeline script, they will be stored in spikeTimes_samples. To get them into seconds, divide spikeTimes_samples by your ephys sampling rate (usually 30000 samples/second).

[simon37robledo]

Thanks for the answers! It is clear now.

[Julie-Fabre]

Great! Closing the issue now. Don't hesitate to open it back up or create a new one if you have any other questions or issues.