Hi Andrew,
I understand that you want to run a t-test on each each entry of the correlation matrices in the 1x8 CTRL and 1x8 TREATED cell.
The following code illustrates how to achieve the same on randomly generated dummy data.
clear;
clc;
% Creating sample correlation matrices for 8 subjects in CTRL and TREATED conditions
CTRL = cell(1, 8);
TREATED = cell(1, 8);
for subject = 1:8
% Generating random 40x40 correlation matrix for each subject
correlationMatrix1 = rand(40);
correlationMatrix2 = rand(40);
% Storing the correlation matrix in the respective cell arrays
CTRL{subject} = correlationMatrix1;
TREATED{subject} = correlationMatrix2;
end
% Initialize the output matrix
outputMatrix = zeros(40);
% Iterate over each entry of the correlation matrices
for i = 1:40
for j = 1:40
% Generate arrays with (i, j)th entry from each subject's correlation matrix
ctrlValues = cellfun(@(x) x(i, j), CTRL);
treatedValues = cellfun(@(x) x(i, j), TREATED);
% Perform the t-test between the two arrays
[~, pValue] = ttest2(ctrlValues, treatedValues);
% Check if the p-value is significant (e.g., pValue < 0.05)
if pValue < 0.05
outputMatrix(i, j) = 1; % Mark as significant
end
end
end
% Resulting matrix
disp(size(outputMatrix));
Please refer to the follow documentation links for information on the functions used in the code above.
- cellfun - https://www.mathworks.com/help/matlab/ref/cellfun.html
- ttest2 - https://www.mathworks.com/help/stats/ttest2.html
Hope this resolves your query.
