123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382 |
- %% 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]);
|