clear; clc;
warning('off');

%% 
load('ach_non.mat');
load('ach_sti.mat');

ca1 = [128 193 219]./255;    % Blue theme
cl1 = [ 52 148 186]./255;

ca2 = [243 169 114]./255;    % Orange theme
cl2 = [236 112  22]./255;

iosr.statistics.boxPlot([1:4]-0.2,non_data,'theme','colorall',...
            'showScatter',true,'scatterMarker','.','symbolMarker','none',...
            'scatterColor',cl1 ,'scatterSize', 128, 'scatterAlpha', 0.5,...
            'boxWidth',0.4,'boxColor',ca1,'showViolin',false); hold on;
iosr.statistics.boxPlot([1:4]+0.2,sti_data,'theme','colorall',...
            'showScatter',true,'scatterMarker','.','symbolMarker','none',...
            'scatterColor',cl2 ,'scatterSize', 128, 'scatterAlpha', 0.5,...
            'boxWidth',0.4,'boxColor',ca2,'showViolin',false); hold on;
        
xticks([1,2,3,4]); xticklabels({'delay','center','left','right'});
ylabel('dFoverF (z-scored)'); 