close all;
clear;
% plot classification in response to natural sounds
int=load('all_data.mat'); int=int.all_di;
chi=load('chi_data.mat'); chi=chi.all_di;
pt=load('pt_data.mat');   pt=pt.all_di;
pair=load('pair_data.mat');   pair=pair.all_di;

awake_=int(:,3)==0;  
int_1=int(:,3)==1;
int_2=pt(:,3)==1;
chi_1=chi(:,3)==0;
chi_2=pt(:,3)==2;
pt_=pt(:,3)==3;
pair_=pair(:,3)==0;

% classification lists
awake=int(awake_,2);
anes_intact=[int(int_1,2); pt(int_2,2)];
chi_alone=chi(chi_1,2);
pair_chi=pt(chi_2,2);
pair_inj=pair(pair_,2);


%% plotting
grey_light=[0.6 0.6 0.6];
grey=[0.3 0.3 0.3];

figure(1); set(gcf, 'Position',  [200, 400, 700, 180])
subplot(151)
y=histc(awake,0:4);
% y=y/max(y)*100;
b=bar(y([2,3,4,5,1],1),'FaceColor', grey_light);
b.FaceColor = 'flat';
b.CData(5,:) = grey;
get_y = get(gca, 'YTick');
y_pos=max(get_y)*0.75;
y_pos2=max(get_y)*0.60;
text(1,y_pos,['n=' num2str(y(1))],'FontSize',8,'Color',grey);
text(1,y_pos2,['n=' num2str(sum(y))],'FontSize',8,'Color',grey_light);

xticklabels({'ech','com','ech>com','com>ech','ech=com'})
ylabel('# of units')
title('\fontsize{8}awake','FontWeight', 'Normal','FontSize',8)

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

subplot(152)
y=histc(anes_intact,0:4);
% y=y/max(y)*100;
b=bar(y([2,3,4,5,1],1),'FaceColor',grey_light);
b.FaceColor = 'flat';
b.CData(5,:) =grey;
get_y = get(gca, 'YTick');
y_pos=max(get_y)*0.75;
y_pos2=max(get_y)*0.60;
text(1,y_pos,['n=' num2str(y(1))],'FontSize',8,'Color',grey);
text(1,y_pos2,['n=' num2str(sum(y))],'FontSize',8,'Color',grey_light);

xticklabels({'ech','com','ech>com','com>ech','ech=com'})
title('\fontsize{8}anesthetized','FontWeight', 'Normal')

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


subplot(153)
y=histc(chi_alone,0:4);
% y=y/max(y)*100;
b=bar(y([2,3,4,5,1],1),'FaceColor',grey_light);
b.FaceColor = 'flat';
b.CData(5,:) = grey;
get_y = get(gca, 'YTick');
y_pos=max(get_y)*0.75;
y_pos2=max(get_y)*0.60;
text(1,y_pos,['n=' num2str(y(1))],'FontSize',8,'Color',grey);
text(1,y_pos2,['n=' num2str(sum(y))],'FontSize',8,'Color',grey_light);

xticklabels({'ech','com','ech>com','com>ech','ech=com'})
title('\fontsize{8}chimeras alone','FontWeight', 'Normal')

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

subplot(154)
y=histc(pair_chi,0:4);
b=bar(y([2,3,4,5,1],1),'FaceColor',grey_light);
b.FaceColor = 'flat';
b.CData(5,:) = grey;
get_y = get(gca, 'YTick');
y_pos=max(get_y)*0.75;
y_pos2=max(get_y)*0.60;
text(1,y_pos,['n=' num2str(y(1))],'FontSize',8,'Color',grey);
text(1,y_pos2,['n=' num2str(sum(y))],'FontSize',8,'Color',grey_light);

xticklabels({'ech','com','ech>com','com>ech','ech=com'})
title('\fontsize{8}paired chimeras','FontWeight', 'Normal')

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


subplot(155)
y=histc(pair_inj,0:4);
% y=y/max(y)*100;
b=bar(y([2,3,4,5,1],1),'FaceColor',grey_light);
b.FaceColor = 'flat';
b.CData(5,:) = grey;
get_y = get(gca, 'YTick');
y_pos=max(get_y)*0.75;
y_pos2=max(get_y)*0.60;
text(1,y_pos,['n=' num2str(y(1))],'FontSize',8,'Color',grey);
text(1,y_pos2,['n=' num2str(sum(y))],'FontSize',8,'Color',grey_light);

xticklabels({'ech','com','ech>com','com>ech','ech=com'});
title('\fontsize{8}paired injection','FontWeight', 'Normal')

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

% exportgraphics(gcf,'E:\Users\User\Desktop\delay paper\anesthetized\figures_paper\class_sets.pdf',...
%     'Resolution',300','ContentType','vector','BackgroundColor','none')