Skip to main content
Revision bc6bc5271ca8c2cbf8107bdf243352e543f64b8e authored by Brian Angeles (Desktop) on 20 March 2023, 22:34:45 UTC, committed by Brian Angeles (Desktop) on 20 March 2023, 22:34:45 UTC
Added .gitignore file and removed RigChoice.s2s
As a first step to generalizing the code, we removed the need to
comment/uncomment the appropriate rig-specific parameters file and
replace it with a local rig_config.s2s file that is ignored by Git. This
way, changes can be made to this simple file without affecting the
Github repository at all.
1 parent 3d78fa7
Raw File Download
Boris_spikeAnalysis.m 
function Boris_spikeAnalysis(ephysdata, struct, thresh)

%temp
%thresh = 120;

% Prep & params
window = [20 40];

% manual spike time extraction
putitive_spikes = ephysdata > thresh;
rc = size(putitive_spikes);
RC = rc(1) * rc(2);
spikeVec = reshape(putitive_spikes', 1,RC);

% create spike window
for j = 1:length(spikeVec)
    if spikeVec(j) == 1
        % block out the rest of spike. Only keep initial detection point
        spikeVec(j+1:j+window(2)-1) = 0;
    end
end

putitive_spikes_mat = vec2mat(spikeVec, struct.segLen * 100);

% x1 bins (raw)
putitive_spikes_mat_x1 = sum(putitive_spikes_mat);
% x2 bins
putitive_spikes_mat_x2 = putitive_spikes_mat_x1(1:2:end) + putitive_spikes_mat_x1(2:2:end);
% x4 bins
putitive_spikes_mat_x4 = putitive_spikes_mat_x2(1:2:end) + putitive_spikes_mat_x2(2:2:end);
% x8 bins
putitive_spikes_mat_x8 = putitive_spikes_mat_x4(1:2:end) + putitive_spikes_mat_x4(2:2:end);


%% plotting
figure()
% subplot(2,1,1)
bar([1:4:struct.segLen * 100]./100, putitive_spikes_mat_x4); hold on
xlim([0 struct.segLen * 100]./100)
% subplot(2,1,2)
%plot(ephysdata'); hold on
% plot(nanmean(ephysdata), 'k');
%hline(thresh, ':k')
xlabel('Seconds')
Boris_figureCosmetics



end