anesthesia_bats/Ephys data/Fig7D-E.m

%% recovery plots: comparing 60 vs 416 ms 
% col 1 --> unit ID 
% col 2 --> neuronal classification
% col 3 --> condition 0: awake / 1: anesthetized
% col 4 --> discriminability index isolation - 60 ms
% col 5 --> discriminability index navigation - 60 ms
% col 6 --> discriminability index communication - 60 ms
% col 7 --> discriminability index isolation - 416 ms
% col 8 --> discriminability index navigation - 416 ms
% col 9 --> discriminability index communication - 416 ms
% col 10 --> context effect nav-exp - 60 ms
% col 11 --> context effect nav-unexp - 60 ms
% col 12 --> context effect com-exp - 60 ms
% col 13 --> context effect com-unexp - 60 ms
% col 14 --> context effect nav-exp - 416 ms
% col 15 --> context effect nav-unexp - 416 ms
% col 16 --> context effect com-exp - 416 ms
% col 17 --> context effect com-unexp - 416 ms
% col 18 --> mean spikecount isolation nav 
% col 19 --> mean spikecount isolation com 
%%
load('all_data.mat')
% add data from injection data set 
inj=load('inj_data.mat');   inj=inj.all_di;
all_di=[all_di; inj];

% add data from chimeras data set 
all_chi=load('chi_data.mat'); all_chi=all_chi.all_di;
all_chi(1:33,3)=2;
all_di=[all_di; all_chi(1:33,:)];
%%
c_=all_di(:,2)==0;   % unit classification
aw_=all_di(:,3)==0; aw_(89:end)=0;  % awake data 
an_=all_di(:,3)==1;                 % anesthetized data 
chi_=all_di(:,3)==2;                % chimeras data 
p1=logical(c_.*aw_);
p2=logical(c_.*an_);
p3=logical(c_.*chi_);

% navigation
% ce exp-60-416 unexp-60-416 awake
aw1=all_di(p1,10);
aw2=all_di(p1,14);
aw3=all_di(p1,11);
aw4=all_di(p1,15);
awN=[aw1 aw2 aw3 aw4];
% ce exp-60-416 unexp-60-416 anesthetized
an1=all_di(p2,10);
an2=all_di(p2,14);
an3=all_di(p2,11);
an4=all_di(p2,15);
anN=[an1 an2 an3 an4];
% ce exp-60-416 unexp-60-416 chimeras
chi1=all_di(p3,10);
chi2=all_di(p3,14);
chi3=all_di(p3,11);
chi4=all_di(p3,15);
chiN=[chi1 chi2 chi3 chi4];


% communication
% ce exp-60-416 unexp-60-416 awake
aw1=all_di(p1,12);
aw2=all_di(p1,16);
aw3=all_di(p1,13);
aw4=all_di(p1,17);
awC=[aw1 aw2 aw3 aw4];
% ce exp-60-416 unexp-60-416 anesthetized
an1=all_di(p2,12);
an2=all_di(p2,16);
an3=all_di(p2,13);
an4=all_di(p2,17);
anC=[an1 an2 an3 an4];
% ce exp-60-416 unexp-60-416 chimeras
chi1=all_di(p3,12);
chi2=all_di(p3,16);
chi3=all_di(p3,13);
chi4=all_di(p3,17);
chiC=[chi1 chi2 chi3 chi4];

% compare paired samples 
pv1=signrank(awN(:,1),awN(:,2));
pv2=signrank(awN(:,3),awN(:,4));
pv3=signrank(awC(:,1),awC(:,2));
pv4=signrank(awC(:,3),awC(:,4));
pv13=signrank(awN(:,2),awN(:,4));
pv14=signrank(awC(:,2),awC(:,4));

disp(CliffDelta(awN(:,3),awN(:,4)))
disp(CliffDelta(awC(:,3),awC(:,4)))
disp(CliffDelta(anN(:,3),anN(:,4)))
disp(CliffDelta(anC(:,3),anC(:,4)))

pv5=signrank(anN(:,1),anN(:,2));
pv6=signrank(anN(:,3),anN(:,4));
pv7=signrank(anC(:,1),anC(:,2));
pv8=signrank(anC(:,3),anC(:,4));
pv15=signrank(anN(:,2),anN(:,4));
pv16=signrank(anC(:,2),anC(:,4));

pv9=signrank(chiN(:,1),chiN(:,2));
pv10=signrank(chiN(:,3),chiN(:,4));
pv11=signrank(chiC(:,1),chiC(:,2));
pv12=signrank(chiC(:,3),chiC(:,4));
pv17=signrank(chiN(:,2),chiN(:,4));
pv18=signrank(chiC(:,2),chiC(:,4));

%% plots for awake data - intact - both gaps 
red_light =[0.8500,0.3250, 0.0980];
blue_light=[0,0.4470,0.7410];
red=[1,0.5,0.5].*red_light;
blue=[1,0.5,0.5].*blue_light;
green_light=[0.4863,0.3059,0.7098];    % green: [0.3922,0.8314,0.0745];
green=[0.2980,0.1294,0.5098];     % [0.4667,0.6745,0.1882];

colors=[blue_light;blue_light;blue;blue];
% colors=[green_light;green_light;green;green];

x=1:4;

figure(1); set(gcf,'Position',[400 200 200 250])
t=tiledlayout(2,1);
ax1 = nexttile;

hold(ax1);
for i = 1:4
    h=boxchart(x(i)*ones(size(awN(:,i))), awN(:,i), 'BoxFaceColor', colors(i,:),'Notch','on');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on
plot(1:2,median(awN(:,1:2)),'-o','LineWidth',1,'Color',colors(1,:),'MarkerSize',3,'MarkerFaceColor',colors(1,:))
plot(3:4,median(awN(:,3:4)),'-o','LineWidth',1,'Color',colors(3,:),'MarkerSize',3,'MarkerFaceColor',colors(3,:))
hold off

xlim([0.5 4.5])
ylabel({'echolocation'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
line1=0.6;
line2=1;

sl=0.75; 

hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);

ylim([-1.2 1.2])

% plot level of significance
if pv1 <= 0.05 && pv1 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv1 <=0.01 && pv1 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv1 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

if pv2 <= 0.05 && pv2 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv2 <=0.01 && pv2 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv2 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

if pv13 <= 0.05 && pv13 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv13 <=0.01 && pv13 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv13 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

ax2 = nexttile;
hold(ax2);
for i = 1:4
    h=boxchart(x(i)*ones(size(awC(:,i))), awC(:,i), 'BoxFaceColor', colors(i,:),'Notch','on');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on 
plot(1:2,median(awC(:,1:2)),'-o','LineWidth',1,'Color',colors(1,:),'MarkerSize',3,'MarkerFaceColor',colors(1,:))
plot(3:4,median(awC(:,3:4)),'-o','LineWidth',1,'Color',colors(3,:),'MarkerSize',3,'MarkerFaceColor',colors(3,:))
hold off
xlim([0.5 4.5])
ylabel({'communication'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);


ylim([-1.2 1.2])

% plot level of significance
if pv3 <= 0.05 && pv3 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv3 <=0.01 && pv3 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv3 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

if pv4 <= 0.05 && pv4 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv4 <=0.01 && pv4 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv4 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

if pv14 <= 0.05 && pv14 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv14 <=0.01 && pv14 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv14 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',8,'HorizontalAlignment','center')
end

linkaxes([ax1,ax2],'x');
xlabel(t,'gap','FontSize',8)
% ylabel(t,'Context effect','FontSize',8)
xticklabels(ax1,{});
xticklabels(ax2,{'60','416','60','416'});
t.TileSpacing = 'compact';



%% plots for anesthetized data - intact - both gaps 
colors=[red_light;red_light;red;red];

x=1:4;

figure(2); set(gcf,'Position',[600 200 200 250])
t=tiledlayout(2,1);
ax1 = nexttile;
hold(ax1);

for i = 1:4
    h=boxchart(x(i)*ones(size(anN(:,i))), anN(:,i), 'BoxFaceColor', colors(i,:),'Notch','on');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on
plot(1:2,median(anN(:,1:2)),'-o','LineWidth',1,'Color',red_light,'MarkerSize',3,'MarkerFaceColor',red_light)
plot(3:4,median(anN(:,3:4)),'-o','LineWidth',1,'Color',red,'MarkerSize',3,'MarkerFaceColor',red)
hold off

xlim([0.5 4.5])
ylabel({'echolocation'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
line1=1;
line2=1.4;

sl=0.75; 

hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);

ylim([-1.2 1.5])

% plot level of significance
if pv5 <= 0.05 && pv5 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv5 <=0.01 && pv5 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv5 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv6 <= 0.05 && pv6 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv6 <=0.01 && pv6 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv6 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv15 <= 0.05 && pv15 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv15 <=0.01 && pv15 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv15 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end


ax2 = nexttile;
hold(ax2);

for i =1:4
    h=boxchart(x(i)*ones(size(anC(:,i))), anC(:,i), 'BoxFaceColor', colors(i,:),'Notch','on');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on 
plot(1:2,median(anC(:,1:2)),'-o','LineWidth',1,'Color',red_light,'MarkerSize',3,'MarkerFaceColor',red_light)
plot(3:4,median(anC(:,3:4)),'-o','LineWidth',1,'Color',red_light,'MarkerSize',3,'MarkerFaceColor',red_light)
hold off

xlim([0.5 4.5])
ylabel({'communication'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
line1=0.8;
line2=1.2;

sl=0.75; 

hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);

ylim([-1.2 1.4])

% plot level of significance
if pv7 <= 0.05 && pv7 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv7 <=0.01 && pv7 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv7 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv8 <= 0.05 && pv8 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv8 <=0.01 && pv8 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv8 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv16 <= 0.05 && pv16 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv16 <=0.01 && pv16 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv16 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end


linkaxes([ax1,ax2],'x');
xlabel(t,'gap','FontSize',8)
% ylabel(t,'Context effect','FontSize',8)
xticklabels(ax1,{});
xticklabels(ax2,{'60','416','60','416'});
t.TileSpacing = 'compact';




%% plot of chimeras data - anesthetized

colors=[green_light;green_light;green;green];

x=1:4;

figure(3); set(gcf,'Position',[800 200 200 250])
t=tiledlayout(2,1);
ax1 = nexttile;

hold(ax1);
for i = 1:4
    h=boxchart(x(i)*ones(size(chiN(:,i))), chiN(:,i), 'BoxFaceColor', colors(i,:),'Notch','off');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on
plot(1:2,median(chiN(:,1:2)),'-o','LineWidth',1,'Color',colors(1,:),'MarkerSize',3,'MarkerFaceColor',colors(1,:))
plot(3:4,median(chiN(:,3:4)),'-o','LineWidth',1,'Color',colors(3,:),'MarkerSize',3,'MarkerFaceColor',colors(3,:))
hold off

xlim([0.5 4.5])
ylabel({'fast echo'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
line1=0.25;
line2=0.65;

sl=0.75; 

hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);

ylim([-1.2 0.8])

% plot level of significance
if pv9 <= 0.05 && pv9 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv9 <=0.01 && pv9 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv9 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv10 <= 0.05 && pv10 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv10 <=0.01 && pv10 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv10 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv17 <= 0.05 && pv17 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv17 <=0.01 && pv17 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv17 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end


ax2 = nexttile;
hold(ax2);
for i = 1:4
    h=boxchart(x(i)*ones(size(chiC(:,i))), chiC(:,i), 'BoxFaceColor', colors(i,:),'Notch','off');
    h.MarkerStyle='.';
    h.MarkerColor=colors(i,:);
end
hold on 
plot(1:2,median(chiC(:,1:2)),'-o','LineWidth',1,'Color',colors(1,:),'MarkerSize',3,'MarkerFaceColor',colors(1,:))
plot(3:4,median(chiC(:,3:4)),'-o','LineWidth',1,'Color',colors(3,:),'MarkerSize',3,'MarkerFaceColor',colors(3,:))
hold off
xlim([0.5 4.5])
ylabel({'slow com'; 'effect'})

set(gca, 'box', 'off')
set(gca, 'Color','none')
set(gca,'linewidth',1);set(gca,'fontsize',8);

% plot lines for comparison 
hold on 
plot([1 2], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*1,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*2,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([3 4], [1 1]*line1, '-k','LineWidth',sl) 
plot([1;1]*3,[line1,line1-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line1,line1-0.1],'-k', 'LineWidth',sl);

plot([2 4], [1 1]*line2, '-k','LineWidth',sl) 
plot([1;1]*2,[line2,line2-0.1],'-k', 'LineWidth',sl);
plot([1;1]*4,[line2,line2-0.1],'-k', 'LineWidth',sl);

ylim([-1.2 0.8])

% plot level of significance
if pv11 <= 0.05 && pv11 > 0.01
    text(mean(1:2)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv11 <=0.01 && pv11 > 0.001
    text(mean(1:2)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv11 <= 0.001
    text(mean(1:2)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(1:2)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv12 <= 0.05 && pv12 > 0.01
    text(mean(3:4)-0.05, line1+0.05, '*','HorizontalAlignment','center')
elseif pv12 <=0.01 && pv12 > 0.001
    text(mean(3:4)-0.05, line1+0.05, '**','HorizontalAlignment','center')
elseif pv12 <= 0.001
    text(mean(3:4)-0.05, line1+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(3:4)-0.05, line1+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

if pv18 <= 0.05 && pv18 > 0.01
    text(mean(2:4)-0.05, line2+0.05, '*','HorizontalAlignment','center')
elseif pv18 <=0.01 && pv18 > 0.001
    text(mean(2:4)-0.05, line2+0.05, '**','HorizontalAlignment','center')
elseif pv18 <= 0.001
    text(mean(2:4)-0.05, line2+0.05, '***','HorizontalAlignment','center')
else 
    text(mean(2:4)-0.05, line2+0.2, 'ns','FontSize',7,'HorizontalAlignment','center')
end

linkaxes([ax1,ax2],'x');
xlabel(t,'gap','FontSize',8)
% ylabel(t,'Context effect','FontSize',8)
xticklabels(ax1,{});
xticklabels(ax2,{'60','416','60','416'});
t.TileSpacing = 'compact';