-
Notifications
You must be signed in to change notification settings - Fork 0
/
tylers_similarity.m
31 lines (30 loc) · 1.25 KB
/
tylers_similarity.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
function outString = tyler_similarity(fixMaps, removedAt, t)
% Initialization
outString = '';
fixMapsSize = size(fixMaps);
% For each participant...
for i=1:fixMapsSize(3)
% If it was one we skipped due to missing data, just write N/A.
if (ismember(i,removedAt))
outString = [outString 'N/A,'];
else
holder = load(['young_IVT_MovingMedian_avg/young_avg_' num2str(t) '.mat']);
meanMatrix = holder.meanMatrix;
% Reshape 2D data to 1D samples and meld together
% Done with magic matlab syntax
linearMeanMatrix = meanMatrix(:);
linearMeanMatrix = linearMeanMatrix';
% linearMeanMatrix = linearMeanMatrix'; % Just flipping again - trust me
linearIndivMatrix = fixMaps(:,:,i);
linearIndivMatrix = linearIndivMatrix(:);
linearIndivMatrix = linearIndivMatrix';
statsMatrix = cat(2,linearIndivMatrix',linearMeanMatrix');
% Do the stats!
cMat=corr(statsMatrix);
cMat=cMat(find(triu(ones(size(cMat)),1)));
r=tanh(nanmean(atanh(cMat)));
% Save to output string
outString = [outString num2str(r) ','];
end
end
end