% run RSA analysis on whole video set
%% paths and files
clear
addpath(fullfile(pwd,'functions'))
addpath(genpath(fullfile(fileparts(pwd),'utils')))
%where to load data and store results
analysispath = '/Users/dianadima/OneDrive - Johns Hopkins/Desktop/MomentsInTime/mot_action/results/video_ratings';
loadfile = fullfile(analysispath,'videoset_full.mat');
savefile = fullfile(analysispath,'rsamodels_full.mat');
%% get models
[models, modelnames] = get_rsa_models(loadfile,savefile);
%% plot models
load(loadfile,'categories','categories_idx')
plot_rsa_models(models,modelnames,'Model RDMs',categories, categories_idx)
%% model corrs
nmod = numel(modelnames);
mcorrS = corr(models,'type','Spearman','rows','pairwise');
mcorrK = nan(nmod,nmod);
for i = 1:nmod
mcorrK(i,i) = 1;
for j = i+1:nmod
mcorrK(i,j) = rankCorr_Kendall_taua(models(:,i),models(:,j));
mcorrK(j,i) = mcorrK(i,j);
end
end
figure;plot_rdm(mcorrS,modelnames,[],0,1)
title('Spearman`s model correlations')
figure;plot_rdm(mcorrK,modelnames,[],0,1)
title('Kendall`s tau-A model correlations')
save(savefile,'-append','mcorr*')