Hi Arianna,
I understand that you are trying to find correlation coefficients of three 3D matrices ‘A’, ‘B’, ‘C’ with a 4th 3D matrix ‘D’. And want to know which one of three has maximum correlation with the matrix ‘D’.
Few assumptions that I have made:
- All the 4 matrices A, B, C, D are of same size.
- Correlation coefficients are to be determined along the 2nd dimension i.e., along columns.
Based on the above assumptions, here is a sample script which you can refer to:
%create correlation coefficients cell array where each element of array i.e cca{i,j} is
%a square matrix which represents correlation of ith column of matrix P to
%jth column of matrix R
cca = cell(3,10);
ccb = cell(3,10);
ccc = cell(3,10);
% matrixSize is an array that holds the dimensions of the matrices
% Calculate correlation coefficients over time (2nd dimension)
for i=1:matrixSize(3)
for t = 1:matrixSize(2)
cca{i,t} = corrcoef(reshape(A(:, t, :), [], matrixSize(3)), reshape(D(:, t, :), [], matrixSize(3)));
ccb{i,t} = corrcoef(reshape(B(:, t, :), [], matrixSize(3)), reshape(D(:, t, :), [], matrixSize(3)));
ccc{i,t} = corrcoef(reshape(C(:, t, :), [], matrixSize(3)), reshape(D(:, t, :), [], matrixSize(3)));
end
end
%Now just iterate over each of the cell arrays and store all the
%coefficients other than 1 into a vector lets say iterate over cca to
%get a vector coeffA
%Once all 3 vectors are obtained coeffA,coeffB,coeffC just find max for each
%of three and use them according to use-case.
For further details, refer the following documentations:
