/
plot_pairwise.m
69 lines (54 loc) · 1.94 KB
/
plot_pairwise.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
function plot_pairwise(PWD,ROI,SUB_LIST,VOX_SIZE,MAX_CL_NUM,LorR)
if LorR == 1
LR='L';
elseif LorR == 0
LR='R';
end
sub=textread(SUB_LIST,'%s');
sub_num=length(sub);
file=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_index_pairwise.mat');
v=load(file);
x=2:MAX_CL_NUM;
mat_cv=[];
mat_dice=[];
mat_nmi=[];
mat_vi=[];
mask=triu(ones(sub_num),1); % upper triangular part
for kc=2:MAX_CL_NUM
col_cv=v.cv(:,:,kc);
mat_cv(:,kc)=col_cv(find(mask));
col_dice=v.dice(:,:,kc);
mat_dice(:,kc)=col_dice(find(mask));
col_nmi=v.nminfo(:,:,kc);
mat_nmi(:,kc)=col_nmi(find(mask));
col_vi=v.vi(:,:,kc);
mat_vi(:,kc)=col_vi(find(mask));
end
hold on;
errorbar(x,mean(mat_dice(:,2:kc)),std(mat_dice(:,2:kc)),'-r','Marker','*');
errorbar(x,mean(mat_nmi(:,2:kc)),std(mat_nmi(:,2:kc)),'-b','Marker','*');
errorbar(x,mean(mat_cv(:,2:kc)),std(mat_cv(:,2:kc)),'-g','Marker','*');
hold off;
set(gca,'XTick',x);
legend('Dice','NMI','CV','Location','SouthEast');
xlabel('Number of clusters','FontSize',14);ylabel('Indice','FontSize',14);
title(strcat(ROI,'.',LR,' pairwise'),'FontSize',14);
set(gcf,'Color','w');
output=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_pairwise.jpg');
export_fig(output,'-r300','-painters','-nocrop');
close;
% VI with non-significant label
errorbar(x,mean(mat_vi(:,2:kc)),std(mat_vi(:,2:kc)),'-r','Marker','.');
for k=2:MAX_CL_NUM-1
h=ttest2(v.vi(:,k),v.vi(:,k+1),0.05,'left');
if h==0
sigstar({[k,k+1]},[nan]);
end
end
set(gca,'XTick',x);
xlabel('Number of clusters','FontSize',14);ylabel('VI','FontSize',14);
title(strcat(ROI,'.',LR,' pairwise VI'),'FontSize',14);
set(gcf,'Color','w');
output=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_pairwise_vi.jpg');
export_fig(output,'-r300','-painters','-nocrop');
close;