ayalab1/neurocode

Problematic use of UID to get unique units in RankOrder.m

ryanharvey1 opened this issue · 0 comments

If restricting data prior to call to RankOrder.m, rank order will be incorrect.

spikes = importSpikes('brainRegion',"CA1",'cellType',"Pyramidal Cell");

spikes

UID: [1×240 double]
times: {1×240 cell}

Only 240 units. However,

[spkEventTimes] = getRipSpikes('spikes',spikes,'events',eventInt);

[rankStats] = RankOrder('spkEventTimes',spkEventTimes)

size(rankUnits)
ans =
275 4803

units = unique(evtTimes{event}(2,:),'stable');