Homec4science
Diffusion ISFC_Pipeline (master) Browse

Intersubj_pipeline/JOVE_Plotting.m
cd4ccb24d88bmaster

Last Change
Raw File

JOVE_Plotting.m
View Options

%% This script highlights how the figures from the article were generated
n_ROI = 299;
n_conn = 299*298/2;
n_sw_tp = 167;
n_subjects = length(TC_MW);
% Adds paths
addpath(genpath(fullfile(pwd,'Utilities')));
% We select three connections of interest that we will compare across
% parameters
COI = [39,207; 148,207; 207,248];
COI_colors{1} = 1/255*[86,214,114;192,214,116;86,214,192];
COI_colors{2} = 1/255*[255,76,0;255,187,0;251,237,0];
% We will plot the same for the different connections
for con = 1:size(COI,1)
%% 1. Sampling of the regions of interest and associated connection
R1 = COI(con,1);
R2 = COI(con,2);
M = zeros(n_ROI,n_ROI);
M(R1,R2) = 1;
Mvec = logical(jUpperTriMatToVec(M));
%% 2. Sampling of the related ISFC time courses
for s = 1:n_subjects
% W = 10
ROI_ISFC_MW(s,:) = ISFC_MW{s}(Mvec,:);
ROI_ISFC_RS(s,:) = ISFC_RS{s}(Mvec,:);
end
t_MW = 19:2:352;
t_RS = 404:2:678;
%% 3. Plotting of the ISFC time courses for the movie-watching and
% resting-state cases
figure;
set(gcf,'color','w');
% For W = 10
subplot(1,2,1);
hold on;
set(gca,'Box','off');
plot(t_MW,ROI_ISFC_MW','LineWidth',0.5,'color',COI_colors{1}(con,:));
hold on; ylim([-0.5 0.5]); xlim([t_MW(1),t_MW(end)]);
xlabel('Time [s]');
ylabel('ISFC [-]');
subplot(1,2,2);
hold on;
set(gca,'Box','off');
plot(t_RS,ROI_ISFC_RS','LineWidth',0.5,'color',COI_colors{2}(con,:));
hold on; ylim([-0.5 0.5]); xlim([t_RS(1),t_RS(end)]);
xlabel('Time [s]');
ylabel('ISFC [-]');
%% 4. Plotting of the ISFC excursion time courses
Plot_Conn_TimeCourse(SignMat,R1,R2,19,352,2,COI_colors{1}(con,:),COI_colors{1}(con,:));
end
%% 5. Plotting of the spatial patterns at a salient time point
idx_W10 = 19;
Spatial_OI = mean(squeeze(SignMat(:,idx_W10,:)),2);
tmp_MAT = jVecToUpperTriMat(Spatial_OI,n_ROI)+jVecToUpperTriMat(Spatial_OI,n_ROI)';
tmp_DEG = sum(abs(tmp_MAT),1)/2;
Factor_plots = 2*max(tmp_DEG);
%%%%%%%%%%%%%%%%%% FREESURFER NEEDS TO BE ADDED TO THE PATH %%%%%%%%%%%%%%%
PlotBrainGraph_Kmeans_withtitle(tmp_MAT,tmp_DEG,CodeBook,0,...
Factor_plots,Factor_plots,2,1,'hot','jet',1,0.5,[-1 1],...
'W10_Alpha001','W10_Alpha001');
c4science · Help