% update for calculating spatial correlations between different channel
% residuals (v3_3)
ResSS_nirs = zeros(nCh, nCh);
for aa = 1:nCh
for bb = 1:nCh
if aa == bb
ResSS_nirs(aa,aa) = ResSS(aa);
elseif aa ~= bb
ResSS_nirs(aa,bb) = (KWY(:,aa)'*KWY(:,bb) - KWY(:,i)'*xX.xKXs.X*xX.pKX* KWY(:,bb)); [This line index i is not defined]
end
end
end
% end of updates


In the file "nirs_spm_indiv.m", the index i is not defined. I think it needs to be aa