Changeset 163 for plugins


Ignore:
Timestamp:
May 18, 2011, 11:32:07 AM (8 years ago)
Author:
tjniskan
Message:
  • Minor changes to partial correlation analysis

M misclib/pcorr.m
M plugins/fmri_plugins/network_analysis.m
M aedes_revision.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • plugins/fmri_plugins/network_analysis.m

    r162 r163  
    6868
    6969% Calculate partial correlations
    70 [C,P] = pcorr(X);
     70[PC,P,CC] = pcorr(X);
    7171
    7272% Linear indexes to unique correlations
    73 ind = find(tril(C,-1)~=0).';
     73ind = find(tril(PC,-1)~=0).';
    7474
    7575% Display results
     
    7878fprintf('***********************************\n');
    7979for ii=ind
    80         [I,J] = ind2sub(size(C),ii);
     80        [I,J] = ind2sub(size(PC),ii);
    8181        if P(ii)<=0.01
    82                 fprintf(2,'%s <-> %s (C=%.4f, P=%.4f)**\n',...
    83                         ROI(J).label,ROI(I).label,C(ii),P(ii));
     82                fprintf(2,'%s <-> %s (PCC=%.4f, CC=%.4f, PCC/CC=%.04f, P=%.4f)**\n',...
     83                        ROI(J).label,ROI(I).label,PC(ii),CC(ii),PC(ii)/CC(ii),P(ii));
    8484        elseif P(ii) <= 0.05
    85                 fprintf(2,'%s <-> %s (C=%.4f, P=%.4f)*\n',...
    86                         ROI(J).label,ROI(I).label,C(ii),P(ii));
     85                fprintf(2,'%s <-> %s (PCC=%.4f, CC=%.4f, PCC/CC=%.04f, P=%.4f)*\n',...
     86                        ROI(J).label,ROI(I).label,PC(ii),CC(ii),PC(ii)/CC(ii),P(ii));
    8787        else
    88                 fprintf(1,'%s <-> %s (C=%.4f, P=%.4f)\n',...
    89                         ROI(J).label,ROI(I).label,C(ii),P(ii));
     88                fprintf(1,'%s <-> %s (PCC=%.4f, CC=%.4f, PCC/CC=%.04f, P=%.4f)\n',...
     89                        ROI(J).label,ROI(I).label,PC(ii),CC(ii),PC(ii)/CC(ii),P(ii));
    9090        end
    9191end
Note: See TracChangeset for help on using the changeset viewer.

Powered by Trac 1.0.9.Copyright © Juha-Pekka Niskanen 2008