Homec4science
Diffusion Research Scripts (Thomas Bolton) (master) Browse

Scripts_Sharing/ISFC/Utils/ISFC_WithinGroup.m
bdfe7aadccccmaster

Last Change
Raw File

ISFC_WithinGroup.m
View Options

%% This function performs within-group ISFC computations between a given
% subject, and all the other subjects from the same group
%
% Inputs:
% - X is the data of all the subjects that we want to consider (cell array with each cell n_ROI x n_TP)
% - W is the window size (in samples)
% - Step is the window Step (in samples)
%
% Outputs:
% - Y is the output (cell array with each cell of size n_pairs x n_swtp)
function [Y,CM_stat] = ISFC_WithinGroup(X,W,Step)
% Number of time points of the data
n_TP = size(X{1},2);
% Number of subjects
n_subjects = length(X);
% Current end of the window; we want to iterate the computations until
% the window end is past the end of the data
WindowEnd = W;
WindowStart = 1;
time_w = 1;
full_index = 1;
% Contains our output
Y = cell(1,n_subjects);
CM_stat = cell(1,n_subjects);
% We will do the same for each sliding window...
while WindowEnd <= n_TP
% For each subject, we want to compute the estimate for a
% given region pair comparing this particular region to all
% the other subject's paired area
for i = 1:n_subjects
% tmp is our data for subject i
tmp = X{i}(:,WindowStart:WindowEnd);
% Takes the indices of all the other subjects
s_idx = 1:n_subjects;
s_idx(i) = [];
tmp_corr = 0;
% We will basically compute all pairs of correlation, and
% average the results
for j = 1:n_subjects-1
to_add = corr(tmp',(X{s_idx(j)}(:,WindowStart:WindowEnd))');
tmp_corr = tmp_corr + ((to_add + to_add')/2);
clear to_add
end
tmp_corr = tmp_corr/(n_subjects - 1);
Y{i}(:,time_w) = jUpperTriMatToVec(tmp_corr);
clear tmp_corr
clear tmp
end
% We shift the window by the step size
WindowStart = WindowStart + Step;
WindowEnd = WindowEnd + Step;
time_w = time_w + 1;
full_index = full_index + 1;
end
% Now we do something similar, but for the stationary part
for i = 1:n_subjects
tmp = X{i};
s_idx = 1:n_subjects;
s_idx(i) = [];
tmp_corr_stat = 0;
for j = 1:n_subjects-1
to_add_stat = corr(tmp',X{s_idx(j)}');
tmp_corr_stat = tmp_corr_stat + ((to_add_stat + to_add_stat')/2);
clear to_add_stat
end
tmp_corr_stat = tmp_corr_stat/(n_subjects - 1);
CM_stat{i} = jUpperTriMatToVec(tmp_corr_stat);
clear tmp_corr_stat
clear tmp
end
end
c4science · Help