-
Notifications
You must be signed in to change notification settings - Fork 2
/
plot_cont.m
37 lines (27 loc) · 1.08 KB
/
plot_cont.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
function plot_cont(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);
file1=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_index_group_continuity.mat');
v1=load(file1);
file2=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_index_indi_continuity.mat');
v2=load(file2);
x=2:MAX_CL_NUM;
m_indi_cont=nanmean(v2.indi_cont);
std_indi_cont=nanstd(v2.indi_cont);
hold on;
plot(x,v1.group_cont(2:end),'-r','Marker','*');
errorbar(x,m_indi_cont(2:end),std_indi_cont(2:end),'-b','Marker','*');
hold off;
set(gca,'XTick',x);
legend('group cont','indi cont','Location','SouthEast');
xlabel('Number of clusters','FontSize',14);ylabel('Indice','FontSize',14);
title(strcat(ROI,'.',LR,' continuity'),'FontSize',14);
set(gcf,'Color','w');
output=strcat(PWD,'/validation_',num2str(sub_num),'_',num2str(VOX_SIZE),'mm/',ROI,'_',LR,'_continuity.jpg');
export_fig(output,'-r300','-painters','-nocrop');
close;