Homec4science
Diffusion Motion HCP analysis (master) Browse

MOT_ANA/Matlab_Scripts/Functions/MOT_ANA_SummariseMotion.m
09f7a984c7ffmaster

Last Change
Raw File

MOT_ANA_SummariseMotion.m
View Options

% MOT has size n_subjects x 36
function [] = MOT_ANA_SummariseMotion(MOT,IDX,Time_Labels,Space_Labels)
n_TB = size(Time_Labels,1);
n_clusters = length(unique(IDX));
n_subjects = size(MOT,1);
% We now want to compute, for all subject/session pairs, representative values
for s = 1:n_subjects
% Data of the considered subject, and associated cluster index
tmp_subject = MOT(s,:);
% Computing the time
for t = 1:n_TB
Summary_time(s,t) = mean(tmp_subject(Time_Labels(t,:)));
end
for p = 1:6
Summary_position(s,p) = mean(tmp_subject(Space_Labels(p,:)));
end
end
figure;
set(gcf,'color','w');
% Along time
to_plot = [];
to_plot2 = [];
% Colors denoting all the 6 time points
colmap = [240,240,240; 220, 220, 220; 190, 190, 190; 160, 160, 160; 130, 130, 130; 100, 100, 100]/255;
for k = 1:n_clusters
for t = 1:n_TB
tmp = Summary_time(IDX == k,t);
%%%%% UNCOMMENT FOR VIOLIN PLOTS GENERATION
% tmpf = figure;
% COL{1} = colmap(t,:);
% COL{2} = [0 0 0];
% MakeViolin(tmp',gca,{''},'',COL,1,1);
% ylim([-2 6.5]);
% print(tmpf,'-depsc',['UCLA_VIOLIN_CLUSTER',num2str(k),'_TIME',num2str(t)]);
% close(tmpf);
to_plot = [to_plot,mean(tmp)];
to_plot2 = [to_plot2,std(tmp)];
end
end
% Computes SEM for display
to_plot2 = to_plot2/sqrt(n_subjects);
subplot(2,1,1);
bar(1:n_TB*n_clusters,to_plot);
hold on;
errorbar(1:n_TB*n_clusters,to_plot,to_plot2);
xlabel('Time x Clusters');
% Along space
to_plot = [];
to_plot2 = [];
colmap = [238,116,97; 149, 197, 106; 129, 149, 204; 169, 219, 241; 190, 135, 187; 254, 244, 161]/255;
for k = 1:n_clusters
for p = 1:6
tmp = Summary_position(IDX == k,p);
%%%%%% UNCOMMENT FOR VIOLIN DISPLAYS
% tmpf = figure;
% COL{1} = colmap(p,:);
% COL{2} = [0 0 0];
% MakeViolin(tmp',gca,{''},'',COL,1,1);
% ylim([-3.5 9]);
% print(tmpf,'-depsc',['UCLA_VIOLIN_CLUSTER',num2str(k),'_POS',num2str(p)]);
% close(tmpf);
to_plot = [to_plot,mean(tmp)];
to_plot2 = [to_plot2,std(tmp)];
end
end
to_plot2 = to_plot2/sqrt(n_subjects);
subplot(2,1,2);
bar(1:6*n_clusters,to_plot);
hold on;
errorbar(1:6*n_clusters,to_plot,to_plot2);
xlabel('Space x Clusters');
to_plot = [];
to_plot2 = [];
end
c4science · Help