fp_figure4.m

function fp_figure4
% Plot different pipelines 

% Copyright (c) 2022 Franziska Pellegrini and Stefan Haufe

DIRDATA = './mim_sim5/';
DIRFIG = './figures/mimsim_ana/mim_sim5/Manuscript/';
if ~exist(DIRFIG); mkdir(DIRFIG); end

labs = {'MIM','MIC','Coherence','iCOH','GC-det','GC-dir','TRGC-det','TRGC-dir'};

%%
clearvars -except iname name DIRDATA DIRFIG labs

%default paramenters
nit = 100;
iInt = 2;
iReg=1;
isnr=0.6;
iss = 0.5;
ilag=2;
ifilt='l';
dimred = 'p';
np = 6;

%% load data

for iit= 1:nit
    
        inname = sprintf('pr_iInt%d_iReg%d_snr0%d_iss0%d_lag%d_filt%s_%s_iter%d'...
                ,iInt,iReg,isnr*10,iss*10, ilag,ifilt,dimred, iit);

        load([DIRDATA inname '.mat'])
        
        PR{1}(iit,:) = pr_mim;
        PR{2}(iit,:) = pr_mic;
        PR{3}(iit,:) = pr_aCoh;
        PR{4}(iit,:) = pr_iCoh;
        PR{5}(iit,:) = pr_absgc;
        PR{6}(iit,:) = pr_posgc;
        PR{7}(iit,:) = pr_abstrgc;
        PR{8}(iit,:) = pr_postrgc;  
end

%% plot 
for icon = [1 8] %plot only for MIM and TRGC
    
    figure
    figone(8,37)
    
    if icon < 3 %MIC, MIM
        pips = [1:8 21 22 10 9]; %fixed, sumvox, central,summim, baseline
        npips = 12;
    else %GC
        pips = [1:8 21 22 9];
        npips = 11;
    end
    
    ipip1=1;
    for ipip = pips
        
        %1 to 6: fixed
        %7: 90%
        %8: 99%
        %9: baseline
        %10: sumVox
        %11: 90% corrected
        %12: 99% corrected
        %13 to 20: equal to 1 to 8 but with zscoring
        %21: sum, then MIM
        %22: central voxel
        
        data1 = PR{icon}(:,ipip);
        mean_pr(icon,ipip1) = mean(data1);
        pr_all(icon,ipip1,:) = data1; 
        imlab1 = 'PR';
        
        %select colors 
        if ipip<=np
            cl = [0.8 0.7 0.6];
        elseif ipip <np+3 || ipip == 11 || ipip == 12
            cl = [0.8 0.4 0.5];
        elseif ipip == 21
            cl = [0.4 0.6 0.7];
        elseif ipip == 22
            cl = [0.5 0.7 0.5];
        elseif ipip==9
            cl = [0.8 0.8 0.8];
        elseif ipip == 10
            cl = [0.8 0.8 0.6]  ;
            cl =[0.4 0.5 0.6];
        end
        
        subplot(1,npips,ipip1)
        
        [h, u] = fp_raincloud_plot_a(data1, cl, 1,0.2, 'ks');
        view([-90 -90]);
        set(gca, 'Xdir', 'reverse');
        set(gca, 'XLim', [0 1]);
        ylim([-0.75 2])
        
        if icon < 3
            xtitles = {'1PC','2PCs', '3PCs','4PCs','5PCs','6PCs','90%',...
                '99%','Mean+FC','Central','FC+mean','TrueVox'};
        else
            xtitles = {'1PC','2PCs', '3PCs','4PCs','5PCs','6PCs','90%',...
                '99%','Mean+FC','Central','TrueVox'};
        end
        
        htit = title(xtitles(ipip1));
        htit.Position(1) = -0.12;
        set(gca,'ytick',[])
        box off
        
        if ipip==1
            xlabel([labs{icon} ' ' imlab1])
            
            set(gca,'Clipping','Off')
            xt = xticks;
            for ix = xt
                if icon <3
                    hu = line([ix ix],[2 -40]);
                else
                    hu = line([ix ix],[2 -36]);
                end
                set(hu, 'color',[0.9 0.9 0.9])
                uistack(hu,'bottom')
            end
            if icon <3
                hu1 = line([0 0],[2 -40]);
            else
                hu1 = line([0 0],[2 -36]);
            end
            set(hu1, 'color',[0 0 0])
            
        else
            set(gca,'xticklabel',{[]})
            set(gca,'XColor','none','YColor','none','TickDir','out')
            set(gca,'Clipping','Off')
            for ix = xt
                hu = line([ix ix],[2.2 -0.75]);
                set(hu, 'color',[0.9 0.9 0.9])
                uistack(hu,'bottom')
            end
            hu = line([0 0],[2.2 -0.75]);
            set(hu, 'color',[0 0 0])
        end
        
        ipip1=ipip1+1;        
    end
    

    %% save
    outname = [DIRFIG name{iname} '_' imlab '_' labs{icon} '.eps'];
    print(outname,'-depsc');
    close all    
    
end

%% test performance differences
%mim 
%3 vs 6 PCs
d3 = squeeze(pr_all(1,3,:)); 
d6 = squeeze(pr_all(1,6,:)); 
[p1,~,stats] = signrank(d6,d3,'tail','left')
t1 = sign(stats.zval);

%3 vs 5 PCs
d5 = squeeze(pr_all(1,5,:)); 
[p2,~,stats] = signrank(d5,d3,'tail','left')

%3 PCs vs truevox 
dt = squeeze(pr_all(1,end,:)); 
[p3,~,stats] = signrank(d3,dt,'tail','left')

%VAR90 vs 3 PCs 
d90 = squeeze(pr_all(1,7,:)); 
[p4,~,stats] = signrank(d90,d3,'tail','left')

%VAR99 vs 3 PCs 
d99 = squeeze(pr_all(1,8,:)); 
[p5,~,stats] = signrank(d99,d3,'tail','left')

%MeanFC vs 3 PCs 
dMFC = squeeze(pr_all(1,9,:)); 
[p6,~,stats] = signrank(dMFC,d3,'tail','left')

%Central vs 3 PCs 
dc = squeeze(pr_all(1,10,:)); 
[p7,~,stats] = signrank(dc,d3,'tail','left')

%FCMEAN vs 3 PCs 
dFCM = squeeze(pr_all(1,11,:)); 
[p8,~,stats] = signrank(dFCM,d3,'tail','left')


%% trgc

%3 vs 6 PCs
d3 = squeeze(pr_all(2,3,:)); 
d6 = squeeze(pr_all(8,6,:)); 
[p1,~,stats] = signrank(d6,d3,'tail','left')

%3 vs 5 PCs
d5 = squeeze(pr_all(2,5,:)); 
[p2,~,stats] = signrank(d5,d3,'tail','left')

%3 PCs vs truevox 
dt = squeeze(pr_all(2,11,:)); 
[p3,~,stats] = signrank(d3,dt,'tail','left')

%VAR90 vs 3 PCs 
d90 = squeeze(pr_all(2,7,:)); 
[p4,~,stats] = signrank(d90,d3,'tail','left')

%VAR99 vs 3 PCs 
d99 = squeeze(pr_all(2,8,:)); 
[p5,~,stats] = signrank(d99,d3,'tail','left')

%MeanFC vs 3 PCs 
dMFC = squeeze(pr_all(2,9,:)); 
[p6,~,stats] = signrank(dMFC,d3,'tail','left')

%Central vs 3 PCs 
dc = squeeze(pr_all(2,10,:)); 
[p7,~,stats] = signrank(dc,d3,'tail','left')