Expand spatial-temporal netmat interaction tests

Add within subject (i.e. to complement across subject).

Scoring/calculateDecompositionAccuracy.m

@@ -224,6 +224,14 @@ scores.pcA_xs.range  = [-1.0, 1.0];
 % Repeat spatial correlations to get per-run metrics
 cPz = repelem(cPz, 1, params.R);
 inf_cPz = repelem(inf_cPz, 1, params.R);
+
+% Each subject, across elements (i.e. are netmats similar)
+% Cosine sim as zero is meaningful for netmats
+scores.cP_pcA.data ...
+    = diag(cosine_sim(inf_cPz, inf_pcAz)) - diag(cosine_sim(cPz, pcAz));
+scores.cP_pcA.name   = 'Spatial-temporal netmat interactions';
+scores.cP_pcA.metric = 'Cosine similarity (difference from GT)';
+scores.cP_pcA.range  = [-2.0, 2.0];
 % Each element, across subjects (i.e. behavioural prediction)
 scores.cP_pcA_xs.data ...
     = diag(corr(inf_cPz', inf_pcAz')) - diag(corr(cPz', pcAz'));