set(0, 'DefaultLineLineWidth', 2);
set(0,'defaultTextFontSize', 12);
set(0,'defaultAxesFontSize',12);
set(0,'DefaultFigureColormap',jet)
set(0,'defaultfigurecolor',[1 1 1])
set(0,'DefaultAxesTitleFontWeight', 'normal')

%%
Ach_change = readtable('ACh_OXT_behavior change.xlsx','Sheet','ACh_behavior');
Oxt_change = readtable('ACh_OXT_behavior change.xlsx','Sheet','OXT_behavior');

%%
G1 = Ach_change.REM*1e2;
%G1 = G1(~isnan(G1));
G2 = Ach_change.NREM*1e2;
%G2 = G2(~isnan(G2));
G3 = Ach_change.WAKE*1e2;
%G3 = G3(~isnan(G3));
figure
set(gcf,'Position',[100 100 400 600])

subplot(2,1,1)
boxplot([G1,G2,G3],'positions',[0,1,2],'notch','on','color',[1 0 0],'symbol','');
box off


G1 = Oxt_change.REM*1e2;
%G1 = G1(~isnan(G1));
G2 = Oxt_change.NREM*1e2;
%G2 = G2(~isnan(G2));
G3 = Oxt_change.WAKE*1e2;
%G3 = G3(~isnan(G3));

subplot(2,1,2)
boxplot([G1,G2,G3],'positions',[0,1,2],'notch','on','color',[0 0.6 0],'symbol','');
box off
