-
Notifications
You must be signed in to change notification settings - Fork 0
/
Dissimilarities_NeuralData.m
64 lines (54 loc) · 1.72 KB
/
Dissimilarities_NeuralData.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
clear;
%% Regular-Irregular spike matrix and the Distances (euclidean)
load('regular_irregular_SPIKES_yannis.mat');
shapenames = regular_irregular_SPIKES_yannis.Properties.VariableNames;
shapenames = shapenames(3:end);
groupnames = {'R','IC','ISC','ISS','ISCa\_ISSa','ISCb\_ISSb'};%,'ISC\_ISS'};%
REGIRREG = table2array(regular_irregular_SPIKES_yannis);
REGIRREG = 4 .* REGIRREG; % (spikes per second?)
% for i=1:119
% REGIRREG(i,:) = REGIRREG(i,:) / max(REGIRREG(i,:));
% end
R_group = REGIRREG(:,3:18);
IC_group = REGIRREG(:,19:34);
ISC_group = REGIRREG(:,35:50);
ISS_group = REGIRREG(:,51:66);
ISCa_ISSa = REGIRREG(:, [35,51,37,53,39,55,41,57,43,59,45,61,47,63,49,65]);
ISCb_ISSb = REGIRREG(:, [36,52,38,54,40,56,42,58,44,60,46,62,48,64,50,66]);
ISC_ISS = horzcat(ISCa_ISSa, ISCb_ISSb);
c = 1;
for i=1:2:16
DistsR(c) = eupairdist(R_group,i);
DistsIC(c) = eupairdist(IC_group ,i);
DistsISC(c) = eupairdist(ISC_group ,i);
DistsISS(c) = eupairdist(ISS_group ,i);
DistsISCa_ISSa(c) = eupairdist(ISCa_ISSa ,i);
DistsISCb_ISSb(c) = eupairdist(ISCb_ISSb ,i);
c = c+1;
end
% c=1;
% for i=1:2:32
% DistsISC_ISS(c) = eupairdist(ISC_ISS ,i);
% c=c+1;
% end
M(:,1) = DistsR;
M(:,2) = DistsIC;
M(:,3) = DistsISC;
M(:,4) = DistsISS;
% M(:,5) = mean([DistsISCb_ISSb; DistsISCa_ISSa]);
M(:,5) = DistsISCa_ISSa;
M(:,6) = DistsISCb_ISSb;
SE = [];
for i=1:length(groupnames)
SE(i) = std(M(:,i))/sqrt(size(M,1)); % standard errors of the mean
end
y = mean(M);
x = 1:length(y);
h = plot(x,y,'o');
axis([0.5,length(groupnames)+0.5,4.5,8.5]);
set(gca,'XTick',x,'XTickLabel',groupnames,'XTickLabelRotation',60)
title('Original Figure from Paper')
grid on
hold on
errorbar(x,y,SE,'o')
hold off