-
Notifications
You must be signed in to change notification settings - Fork 0
/
plotAbrLobes.m
170 lines (136 loc) · 5.67 KB
/
plotAbrLobes.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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
function [statSF,statNSF] = plotAbrLobes(befSurgLobes,aftSurgLobes,metaData)
lobeName = {'Temporal';'SubCortical';'Parietal';'Occipital';'Frontal';'Cingulate'};
%% Compute statistics
for lobe = 1:size(befSurgLobes,1)
[~,statSF.p(lobe,:),CI,stats] = ttest(befSurgLobes(lobe,metaData.ILAE_1==1)',aftSurgLobes(lobe,metaData.ILAE_1==1)');
statSF.meanDiff(lobe,:) = round(mean(befSurgLobes(lobe,metaData.ILAE_1==1))- mean(aftSurgLobes(lobe,metaData.ILAE_1==1)),2);
statSF.CIL(lobe,:) = round(CI(1),2);
statSF.CIU(lobe,:) = round(CI(2),2);
statSF.tstat(lobe,:) = round(stats.tstat,2);
statSF.DF(lobe,:) = stats.df;
[~,statNSF.p(lobe,:),CI,stats] = ttest(befSurgLobes(lobe,metaData.ILAE_1>2)',aftSurgLobes(lobe,metaData.ILAE_1>2)');
statNSF.meanDiff(lobe,:) = round(mean(befSurgLobes(lobe,metaData.ILAE_1>2))- mean(aftSurgLobes(lobe,metaData.ILAE_1>2)),2);
statNSF.CIL(lobe,:) = round(CI(1),2);
statNSF.CIU(lobe,:) = round(CI(2),2);
statNSF.tstat(lobe,:) = round(stats.tstat,2);
statNSF.DF(lobe,:) = stats.df;
end
statSF = struct2table(statSF);
statNSF = struct2table(statNSF);
%% Compute mean
befSurgAbsSF = mean(befSurgLobes(:,metaData.ILAE_1==1),2);
befSurgAbsNSF = mean(befSurgLobes(:,metaData.ILAE_1>2),2);
aftSurgAbsSF = mean(aftSurgLobes(:,metaData.ILAE_1==1),2);
aftSurgAbsNSF = mean(aftSurgLobes(:,metaData.ILAE_1>2),2);
%% Compute sde
befSurgAbsSFsd = std(befSurgLobes(:,metaData.ILAE_1==1),[],2)./numel(metaData.ILAE_1==1);
befSurgAbsNSFsd = std(befSurgLobes(:,metaData.ILAE_1>2),[],2)./numel(metaData.ILAE_1>2);
aftSurgAbsSFsd = std(aftSurgLobes(:,metaData.ILAE_1==1),[],2)./numel(metaData.ILAE_1==1);
aftSurgAbsNSFsd = std(aftSurgLobes(:,metaData.ILAE_1>2),[],2)./numel(metaData.ILAE_1>2);
%% Concatenate sde and mean ILAE 1
meanValSF = [befSurgAbsSF,aftSurgAbsSF]; %befSurgAbsILAE2,befSurgAbsNSF,aftSurgAbsILAE2,aftSurgAbsNSF]; % mean values
sdeValSF = [befSurgAbsSFsd,aftSurgAbsSFsd];%befSurgAbsILAE2sd,befSurgAbsNSFsd,aftSurgAbsSFsd,aftSurgAbsILAE2sd,aftSurgAbsNSFsd]; % standard error
meanValSF_Ipsi = meanValSF(1:6,:);
meanValSF_Contra = meanValSF(7:12,:);
sdeValSF_Ipsi = sdeValSF(1:6,:);
sdeValSF_Contra = sdeValSF(7:12,:);
% Plot Ipsilateral for SF patients
figure
hold on
hb = bar(1:6,meanValSF_Ipsi);
hb(1).FaceColor = hex2rgb('DB984D');
hb(2).FaceColor = hex2rgb('75A5ED');
% For each set of bars, find the centers of the bars, and write error bars
pause(0.1); %pause allows the figure to be created
for ib = 1:numel(hb)
%XData property is the tick labels/group centers; XOffset is the offset
%of each distinct group
xData = hb(ib).XData+hb(ib).XOffset;
errorbar(xData,meanValSF_Ipsi(:,ib),sdeValSF_Ipsi(:,ib),'k.')
end
set(gca,'xTickLabel',lobeName);
set(gca,'xTick',1:6);
ylabel('Proportion abnormal nodes')
set(gca,'FontSize',12);
set(gca,'XTickLabelRotation',45);
legend({'Before Surgery'; 'After Surgery'})
ylim([0 0.3])
legend off
title('Figure 4B: ILAE 1 Ipsilateral')
% Plot Contralateral for SF patients
figure
hold on
hb = bar(1:6,meanValSF_Contra);
hb(1).FaceColor = hex2rgb('DB984D');
hb(2).FaceColor = hex2rgb('75A5ED');
% For each set of bars, find the centers of the bars, and write error bars
pause(0.1); %pause allows the figure to be created
for ib = 1:numel(hb)
%XData property is the tick labels/group centers; XOffset is the offset
%of each distinct group
xData = hb(ib).XData+hb(ib).XOffset;
errorbar(xData,meanValSF_Contra(:,ib),sdeValSF_Contra(:,ib),'k.')
end
set(gca,'xTickLabel',lobeName);
set(gca,'xTick',1:6);
ylabel('Proportion abnormal nodes')
set(gca,'FontSize',12);
set(gca,'XTickLabelRotation',45);
legend({'Before Surgery'; 'After Surgery'})
ylim([0 0.3])
legend off
title('Figure 4B: ILAE 1 Contralateral')
%% Concatenate sde and mean ILAE 3+
meanValNSF = [befSurgAbsNSF,aftSurgAbsNSF]; %befSurgAbsILAE2,befSurgAbsNSF,aftSurgAbsILAE2,aftSurgAbsNSF]; % mean values
sdeValNSF = [befSurgAbsNSFsd,aftSurgAbsNSFsd];%befSurgAbsILAE2sd,befSurgAbsNSFsd,aftSurgAbsSFsd,aftSurgAbsILAE2sd,aftSurgAbsNSFsd]; % standard error
meanValNSF_Ipsi = meanValNSF(1:6,:);
meanValNSF_Contra = meanValNSF(7:12,:);
sdeValNSF_Ipsi = sdeValNSF(1:6,:);
sdeValNSF_Contra = sdeValNSF(7:12,:);
% Plot Ipsilateral for NSF patients
figure
hold on
hb = bar(1:6,meanValNSF_Ipsi);
hb(1).FaceColor = hex2rgb('DB984D');
hb(2).FaceColor = hex2rgb('75A5ED');
% For each set of bars, find the centers of the bars, and write error bars
pause(0.1); %pause allows the figure to be created
for ib = 1:numel(hb)
%XData property is the tick labels/group centers; XOffset is the offset
%of each distinct group
xData = hb(ib).XData+hb(ib).XOffset;
errorbar(xData,meanValNSF_Ipsi(:,ib),sdeValNSF_Ipsi(:,ib),'k.')
end
set(gca,'xTickLabel',lobeName);
set(gca,'xTick',1:6);
ylabel('Proportion abnormal nodes')
set(gca,'FontSize',12);
set(gca,'XTickLabelRotation',45);
legend({'Before Surgery'; 'After Surgery'})
ylim([0 0.3])
legend off
title('Figure 4D: ILAE 3+ Ipsilateral')
% Plot Contralateral for NSF patients
figure
hold on
hb = bar(1:6,meanValNSF_Contra);
hb(1).FaceColor = hex2rgb('DB984D');
hb(2).FaceColor = hex2rgb('75A5ED');
% For each set of bars, find the centers of the bars, and write error bars
pause(0.1); %pause allows the figure to be created
for ib = 1:numel(hb)
%XData property is the tick labels/group centers; XOffset is the offset
%of each distinct group
xData = hb(ib).XData+hb(ib).XOffset;
errorbar(xData,meanValNSF_Contra(:,ib),sdeValNSF_Contra(:,ib),'k.')
end
set(gca,'xTickLabel',lobeName);
set(gca,'xTick',1:6);
ylabel('Proportion abnormal nodes')
set(gca,'FontSize',12);
set(gca,'XTickLabelRotation',45);
legend({'Before Surgery'; 'After Surgery'})
ylim([0 0.3])
legend off
title('Figure 4D: ILAE 3+ Contralateral')
end