%% Network Comparisons
% Calculate covariance measures & visualise networks
% Use covariance measures to derive eigenvectors for PCA
% Perform DR & compare PCs

% Get covariances matrices for all & groups:
g1_cov = cov(g1_coherence); % Group 1 coherence
g2_cov = cov(g2_coherence); % Group 2 coherence

all_pca = pca(all_coherence(:,3:end)');  % Get PCA coefficients for all subjects

g1_pca = all_pca(group1,:); % Isolate group 1
g2_pca = all_pca(group2,:); % Isolate group 2

% Compare networks based on first principal component:
p = ranksum(g1_pca(:,1), g2_pca(:,1)); % Compare first principal components
disp(['Group 1 vs. Group 2: P = ' num2str(p) ', Wilcoxon rank sum test']);