%% Import training data
clear;

load('InstTrainingAnalysis.mat');
load('PavTrainingAnalysis.mat');

%% Figure 1 - Averages across training and boxplots of final performance
figure; 
% Instrumental conditioning: PE latency and PE Probability vs. Day
subplot(3,2,1);
e = errorbar(PAInst.AverageCSPlusCPLatMean,PAInst.Error_AverageCSPlusCPLatMean);
e.Marker = 'o';
e.MarkerSize = 4;
e.MarkerFaceColor = '[0.4940, 0.1840, 0.5560]';
e.Color = '[0.4940, 0.1840, 0.5560]';
ylabel('Port Entry Latency (sec)');
title('Instrumental Conditioning');
hold on
errorbar(PAInst.AverageCSMinCPLatMean,PAInst.Error_AverageCSMinCPLatMean,'-ok','MarkerFaceColor','k','MarkerSize',4);
legend('DS','NS');
axis([0 25 0 25]);

subplot(3,2,3);
e = errorbar(PAInst.AverageCSPlusRatio,PAInst.Error_AverageCSPlusRatio);
hold on;
e.Marker = 'o';
e.MarkerSize = 4;
e.MarkerFaceColor = '[0.4940, 0.1840, 0.5560]';
e.Color = '[0.4940, 0.1840, 0.5560]';
xlabel('Session');
ylabel('Port Entry Probability');
errorbar(PAInst.AverageCSMinRatio,PAInst.Error_AverageCSMinRatio,'-ok','MarkerFaceColor','k','MarkerSize',4);
axis([0 25 0 1]);

% Pavlovian conditioning: PE latency and PE Probability vs. Training Day
subplot(3,2,2);
e = errorbar(PAPav.AverageCSPlusCPLatMean,PAPav.Error_AverageCSPlusCPLatMean);
e.Marker = 'o';
e.MarkerSize = 4;
e.MarkerFaceColor = '[0.9290, 0.6940, 0.1250]';
e.Color = '[0.9290, 0.6940, 0.1250]';
title('Pavlovian Conditioning');
hold on
errorbar(PAPav.AverageCSMinCPLatMean,PAPav.Error_AverageCSMinCPLatMean,'-ok','MarkerFaceColor','k','MarkerSize',4);
legend({'CS+','CS-'},'Location','Northwest');
axis([0 15 0 25]);

subplot(3,2,4);
e = errorbar(PAPav.AverageCSPlusRatio,PAPav.Error_AverageCSPlusRatio);
e.Marker = 'o';
e.MarkerSize = 4;
e.MarkerFaceColor = '[0.9290, 0.6940, 0.1250]';
e.Color = '[0.9290, 0.6940, 0.1250]';
hold on
errorbar(PAPav.AverageCSMinRatio,PAPav.Error_AverageCSMinRatio,'-ok','MarkerFaceColor','k','MarkerSize',4);
axis([0 15 0 1]);
xlabel('Session');

% Calculation of Days to criteria
%PAInst.CriteriaDayRat=cell2mat(PAInst.CriteriaDayRat)
InstSelection=ismember(PAInst.CriteriaDayRat,[105 107 139 151]);
InstCriteria = table(PAInst.CriteriaDay(InstSelection));
nrow = size(InstCriteria, 1);
InstCriteria.Task = ones(nrow, 1);
InstCriteria.Properties.VariableNames = {'Days_to_Criteria','Task'};

%PAPav.CriteriaDayRat=cell2mat(PAPav.CriteriaDayRat)
PavSelection=ismember(PAPav.CriteriaDayRat,[36 122 137 140 154]);
PavCriteria = table(PAPav.CriteriaDay(PavSelection));
nrow = size(PavCriteria, 1);
PavCriteria.Task = repmat(2, nrow, 1);
PavCriteria.Properties.VariableNames = {'Days_to_Criteria','Task'};

Criteria = vertcat(InstCriteria,PavCriteria);

subplot(3,3,7);
names = {'Instrumental'; 'Pavlovian'};
boxplot(Criteria.Days_to_Criteria,Criteria.Task,'colorgroup',[1 2],'colors',[0.4940, 0.1840, 0.5560;0.9290, 0.6940, 0.1250]);
hold on
scatter(InstCriteria.Task,InstCriteria.Days_to_Criteria,15,'c','MarkerEdgeColor',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560]);
hold on
scatter(PavCriteria.Task,PavCriteria.Days_to_Criteria,15,'c','MarkerEdgeColor',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250]);
hold off
title('Days to Criteria')
ylim([0 25]);
set(gca,'xtick',1:2,'xticklabel',names);

% Final port entry probability
PAInst.CriteriaDayCSPlusRatio = cell2mat(PAInst.CriteriaDayCSPlusRatio);
PAInst.CriteriaDayCSMinRatio = cell2mat(PAInst.CriteriaDayCSMinRatio);
PAPav.CriteriaDayCSPlusRatio = cell2mat(PAPav.CriteriaDayCSPlusRatio);
PAPav.CriteriaDayCSMinRatio = cell2mat(PAPav.CriteriaDayCSMinRatio);


%InstSelection=ismember(PAInst.RatNum,[105 106 107 139 151]);
%PavSelection=ismember(PAPav.RatNum,[36 122 137 140 154]);

InstDSPEProb = table(PAInst.CriteriaDayCSPlusRatio(InstSelection));
InstNSPEProb = table(PAInst.CriteriaDayCSMinRatio(InstSelection));
nrowDS = size(InstDSPEProb, 1);
InstDSPEProb.Cue = ones(nrowDS, 1);
nrowNS = size(InstNSPEProb, 1);
InstNSPEProb.Cue = repmat(2,nrowNS, 1);

PavCSPlusPEProb = table(PAPav.CriteriaDayCSPlusRatio(PavSelection));
PavCSMinPEProb = table(PAPav.CriteriaDayCSMinRatio(PavSelection));
nrowCSPlus = size(PavCSPlusPEProb, 1);
PavCSPlusPEProb.Cue = repmat(3, nrowCSPlus, 1);
nrowCSMin = size(PavCSMinPEProb, 1);
PavCSMinPEProb.Cue = repmat(4, nrowCSMin, 1);

FinalPEProb = vertcat(InstDSPEProb,InstNSPEProb,PavCSPlusPEProb,PavCSMinPEProb);
FinalPEProb.Properties.VariableNames = {'Probability','Cue'};

subplot(3,3,8);
names = {'DS'; 'NS'; 'CS+'; 'CS-'};
boxplot(FinalPEProb.Probability,FinalPEProb.Cue,'colorgroup',[1 2 3 4],'colors',[0.4940, 0.1840, 0.5560;0,0,0;0.9290, 0.6940, 0.1250;0,0,0]);
hold on
scatter(InstDSPEProb.Cue,InstDSPEProb.Var1,15,'c','MarkerEdgeColor',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560]);
hold on
scatter(InstNSPEProb.Cue,InstNSPEProb.Var1,15,'c','MarkerEdgeColor',[0,0,0],'MarkerFaceColor',[0,0,0]);
hold on
scatter(PavCSPlusPEProb.Cue,PavCSPlusPEProb.Var1,15,'c','MarkerEdgeColor',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250]);
hold on
scatter(PavCSMinPEProb.Cue,PavCSMinPEProb.Var1,15,'c','MarkerEdgeColor',[0,0,0],'MarkerFaceColor',[0,0,0]);
hold off
title('Port Entry Probability')
ylim([0 1]);
set(gca,'xtick',1:4,'xticklabel',names);

% Final port entry latency (sec)
InstDSLat = table(PAInst.CriteriaDayCSPlusLat(InstSelection));
InstNSLat = table(PAInst.CriteriaDayCSMinLat(InstSelection));
nrowDS = size(InstDSLat, 1);
InstDSLat.Cue = ones(nrowDS, 1);
nrowNS = size(InstNSLat, 1);
InstNSLat.Cue = repmat(2,nrowNS, 1);

PavCSPlusLat = table(PAPav.CriteriaDayCSPlusLat(PavSelection));
PavCSMinLat = table(PAPav.CriteriaDayCSMinLat(PavSelection));
nrowCSPlus = size(PavCSPlusLat, 1);
PavCSPlusLat.Cue = repmat(3,nrowCSPlus, 1);
nrowCSMin = size(PavCSMinLat, 1);
PavCSMinLat.Cue = repmat(4,nrowCSMin, 1);

FinalLat = vertcat(InstDSLat,InstNSLat,PavCSPlusLat,PavCSMinLat);
FinalLat.Properties.VariableNames = {'Latency','Cue'};

subplot(3,3,9);
names = {'DS'; 'NS'; 'CS+'; 'CS-'};
boxplot(FinalLat.Latency,FinalLat.Cue,'colorgroup',[1 2 3 4],'colors',[0.4940, 0.1840, 0.5560;0,0,0;0.9290, 0.6940, 0.1250;0,0,0]);
hold on
scatter(InstDSLat.Cue,InstDSLat.Var1,15,'c','MarkerEdgeColor',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560]);
hold on
scatter(InstNSLat.Cue,InstNSLat.Var1,15,'c','MarkerEdgeColor',[0,0,0],'MarkerFaceColor',[0,0,0]);
hold on
scatter(PavCSPlusLat.Cue,PavCSPlusLat.Var1,15,'c','MarkerEdgeColor',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250]);
hold on
scatter(PavCSMinLat.Cue,PavCSMinLat.Var1,15,'c','MarkerEdgeColor',[0,0,0],'MarkerFaceColor',[0,0,0]);
hold off
title('Port Entry Latency (sec)')
ylim([0 9]);
set(gca,'xtick',1:4,'xticklabel',names);


%% Figure of PE latency and probabilities of individual rats for instrumental conditioning

Rat105 = (PAInst.RatNum == 105 & PAInst.Criteria == 0);
Rat106 = (PAInst.RatNum == 106 & PAInst.Criteria == 0);
Rat107 = (PAInst.RatNum == 107 & PAInst.Criteria == 0);
Rat108 = (PAInst.RatNum == 108 & PAInst.Criteria == 0);
Rat139 = (PAInst.RatNum == 139 & PAInst.Criteria == 0);
Rat151 = (PAInst.RatNum == 151 & PAInst.Criteria == 0);


figure;
subplot(4,2,1);
plot(PAInst.Day(Rat105),PAInst.CSPlusCPLatMean(Rat105),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat106),PAInst.CSPlusCPLatMean(Rat106),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
plot(PAInst.Day(Rat107),PAInst.CSPlusCPLatMean(Rat107),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
plot(PAInst.Day(Rat108),PAInst.CSPlusCPLatMean(Rat108),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
plot(PAInst.Day(Rat139),PAInst.CSPlusCPLatMean(Rat139),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
plot(PAInst.Day(Rat151),PAInst.CSPlusCPLatMean(Rat151),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
title('Instrumental Conditioning')
ylabel('Port Entry Latency (sec)')
legend('DS')
xlim([0 25]);
ylim([0 40]);

subplot(4,2,3);
plot(PAInst.Day(Rat105),PAInst.CSPlusRatio(Rat105),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat106),PAInst.CSPlusRatio(Rat106),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat107),PAInst.CSPlusRatio(Rat107),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat108),PAInst.CSPlusRatio(Rat108),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat139),PAInst.CSPlusRatio(Rat139),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold on
plot(PAInst.Day(Rat151),PAInst.CSPlusRatio(Rat151),'-o','Color',[0.4940, 0.1840, 0.5560],'MarkerFaceColor',[0.4940, 0.1840, 0.5560],'MarkerSize',4);
hold off
ylabel('Port Entry Probability')
xlim([0 25]);
ylim([0 1]);

subplot(4,2,5);
plot(PAInst.Day(Rat105),PAInst.CSMinCPLatMean(Rat105),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat106),PAInst.CSMinCPLatMean(Rat106),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat107),PAInst.CSMinCPLatMean(Rat107),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat108),PAInst.CSMinCPLatMean(Rat108),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat139),PAInst.CSMinCPLatMean(Rat139),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat151),PAInst.CSMinCPLatMean(Rat151),'-ok','MarkerFaceColor','k','MarkerSize',4)
legend('NS','Location','Northeast');
ylabel('Port Entry Latency (sec)');
xlim([0 25]);
ylim([0 10]);

subplot(4,2,7);
plot(PAInst.Day(Rat105),PAInst.CSMinRatio(Rat105),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat106),PAInst.CSMinRatio(Rat106),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat107),PAInst.CSMinRatio(Rat107),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat108),PAInst.CSMinRatio(Rat108),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat139),PAInst.CSMinRatio(Rat139),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAInst.Day(Rat151),PAInst.CSMinRatio(Rat151),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold off
ylabel('Port Entry Probability');
xlabel('Session');
xlim([0 25])
ylim([0 1]);

% PE latency and probabilities of individual rats for pavlovian conditioning

Rat122 = (PAPav.RatNum == 122 & PAPav.Criteria == 0);
Rat123 = (PAPav.RatNum == 123 & PAPav.Criteria == 0);
Rat137 = (PAPav.RatNum == 137 & PAPav.Criteria == 0);
Rat138 = (PAPav.RatNum == 138 & PAPav.Criteria == 0);
Rat140 = (PAPav.RatNum == 140 & PAPav.Criteria == 0);
Rat154 = (PAPav.RatNum == 154 & PAPav.Criteria == 0);
Rat29 = (PAPav.RatNum == 29 & PAPav.Criteria == 0);
Rat30 = (PAPav.RatNum == 30 & PAPav.Criteria == 0);
Rat31 = (PAPav.RatNum == 31 & PAPav.Criteria == 0);
Rat32 = (PAPav.RatNum == 32 & PAPav.Criteria == 0);
Rat33 = (PAPav.RatNum == 33 & PAPav.Criteria == 0);
Rat34 = (PAPav.RatNum == 34 & PAPav.Criteria == 0);
Rat35 = (PAPav.RatNum == 35 & PAPav.Criteria == 0);
Rat36 = (PAPav.RatNum == 36 & PAPav.Criteria == 0);

subplot(4,2,2);
plot(PAPav.Day(Rat122),PAPav.CSPlusCPLatMean(Rat122),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat123),PAPav.CSPlusCPLatMean(Rat123),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat137),PAPav.CSPlusCPLatMean(Rat137),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat138),PAPav.CSPlusCPLatMean(Rat138),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat140),PAPav.CSPlusCPLatMean(Rat140),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat154),PAPav.CSPlusCPLatMean(Rat154),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat29),PAPav.CSPlusCPLatMean(Rat29),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat30),PAPav.CSPlusCPLatMean(Rat30),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat31),PAPav.CSPlusCPLatMean(Rat31),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat32),PAPav.CSPlusCPLatMean(Rat32),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat33),PAPav.CSPlusCPLatMean(Rat33),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat34),PAPav.CSPlusCPLatMean(Rat34),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat35),PAPav.CSPlusCPLatMean(Rat35),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat36),PAPav.CSPlusCPLatMean(Rat36),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold off
xlim([0 25]);
ylim([0 40]);
legend('CS+','Location','Northeast');
title('Pavlovian Conditioning')

subplot(4,2,4);
plot(PAPav.Day(Rat122),PAPav.CSPlusRatio(Rat122),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat123),PAPav.CSPlusRatio(Rat123),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat137),PAPav.CSPlusRatio(Rat137),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat138),PAPav.CSPlusRatio(Rat138),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat140),PAPav.CSPlusRatio(Rat140),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat154),PAPav.CSPlusRatio(Rat154),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat29),PAPav.CSPlusRatio(Rat29),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat30),PAPav.CSPlusRatio(Rat30),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat31),PAPav.CSPlusRatio(Rat31),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat32),PAPav.CSPlusRatio(Rat32),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat33),PAPav.CSPlusRatio(Rat33),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat34),PAPav.CSPlusRatio(Rat34),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat35),PAPav.CSPlusRatio(Rat35),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold on
plot(PAPav.Day(Rat36),PAPav.CSPlusRatio(Rat36),'-o','Color',[0.9290, 0.6940, 0.1250],'MarkerFaceColor',[0.9290, 0.6940, 0.1250],'MarkerSize',4);
hold off
xlim([0 25]);
ylim([0 1]);

subplot(4,2,6);
plot(PAPav.Day(Rat122),PAPav.CSMinCPLatMean(Rat122),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat123),PAPav.CSMinCPLatMean(Rat123),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat137),PAPav.CSMinCPLatMean(Rat137),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat138),PAPav.CSMinCPLatMean(Rat138),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat140),PAPav.CSMinCPLatMean(Rat140),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat154),PAPav.CSMinCPLatMean(Rat154),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat29),PAPav.CSMinCPLatMean(Rat29),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat30),PAPav.CSMinCPLatMean(Rat30),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat31),PAPav.CSMinCPLatMean(Rat31),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat32),PAPav.CSMinCPLatMean(Rat32),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat33),PAPav.CSMinCPLatMean(Rat33),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat34),PAPav.CSMinCPLatMean(Rat34),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat35),PAPav.CSMinCPLatMean(Rat35),'-ok','MarkerFaceColor','k','MarkerSize',4);
plot(PAPav.Day(Rat36),PAPav.CSMinCPLatMean(Rat36),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold off
legend({'CS-'},'Fontsize',4.1,'Position',[0.81878 0.332 0.083 0.02]);
xlim([0 25]);

subplot(4,2,8);
plot(PAPav.Day(Rat122),PAPav.CSMinRatio(Rat122),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat123),PAPav.CSMinRatio(Rat123),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat137),PAPav.CSMinRatio(Rat137),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat138),PAPav.CSMinRatio(Rat138),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat140),PAPav.CSMinRatio(Rat140),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat154),PAPav.CSMinRatio(Rat154),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat29),PAPav.CSMinRatio(Rat29),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat30),PAPav.CSMinRatio(Rat30),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat31),PAPav.CSMinRatio(Rat31),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat32),PAPav.CSMinRatio(Rat32),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat33),PAPav.CSMinRatio(Rat33),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat34),PAPav.CSMinRatio(Rat34),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat35),PAPav.CSMinRatio(Rat35),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold on
plot(PAPav.Day(Rat36),PAPav.CSMinRatio(Rat36),'-ok','MarkerFaceColor','k','MarkerSize',4);
hold off
xlabel('Session');
xlim([0 25]);
ylim([0 1]);