Update 'code/neuralAnalysesOfferShiftAligned.m'

code/neuralAnalysesOfferShiftAligned.m

@@ -35,45 +35,29 @@ disp(['running selflag=' num2str(selflag) '; caseflag=' num2str(caseflag) '; ('
 
 switch selflag
     case 1 
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off1L_del1LR_sacc/';
-%        most_recent_folder=['../../Figures/Temp/NoBaselineNorm/evLevR_off_del_focus_minus200ms/evLevR_off1L_del1LR_' casesstrs{caseflag} '/'];
-        most_recent_folder=['../../Figures/Temp/NoBaselineNorm/evLevR_off_only_focus_minus200ms/evLevR_off1LR_' casesstrs{caseflag} '/'];
+        most_recent_folder=['../../Figures/evLevR_off1LR_' casesstrs{caseflag} '/'];
         if ~exist(most_recent_folder,'dir'); mkdir(most_recent_folder); mkdir([most_recent_folder 'span' num2str(Nspan) 'ms_lag' num2str(Nlag) 'ms']); end
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off1L_del1LR_NOshift_NOoff1L/';
         selectedTimes='((ntce_ds(2,1)-40):ntce_ds(2,2))';
         %selectedCases{1}='find((cciso1LL) & (ccepx(:,tt)<0))';
         %selectedCases{2}='find((cciso1LL) & (ccepx(:,tt)>0))';
-    %case 2
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off1R_del1LR_sacc/';
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off1R_del1LR_sacc/';
-        %selectedTimes='(ntce_ds(3,1):ntce_ds(3,2))';
+    %case 2 % DEPRECATED
+        %most_recent_folder='../../Figures/';
+        %selectedTimes='(ntce_ds(2,1)-40:ntce_ds(2,2))';
         %selectedCases{1}='find((~cciso1LL) & (ccepx(:,tt)<0))';
         %selectedCases{2}='find((~cciso1LL) & (ccepx(:,tt)>0))';
-    %case 3
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off2L_del2LR_sacc/';
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off2L_del2LR_shift/';
-        %selectedTimes='(ntce_ds(5,1):ntce_ds(5,2))';
-        %selectedTimes='(ntce_ds(5,1):ntce_ds(6,1)-2)'; % FULL DELAY2
+    %case 3 % DEPRECATED
+        %most_recent_folder='../../Figures/';
+        %selectedTimes='(ntce_ds(4,1)-40:ntce_ds(4,2))';
         %selectedCases{1}='find((cciso2LL) & (ccepx(:,tt)<0))';
         %selectedCases{2}='find((cciso2LL) & (ccepx(:,tt)>0))';
     case 4
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off2R_del2LR_sacc/';
-        %most_recent_folder=['../../Figures/Temp/NoBaselineNorm/evLevR_off_del_focus_minus200ms/evLevR_off2R_del2LR_' casesstrs{caseflag} '/'];
-        most_recent_folder=['../../Figures/Temp/NoBaselineNorm/evLevR_off_only_focus_minus200ms/evLevR_off2LR_' casesstrs{caseflag} '/'];
+        most_recent_folder=['../../Figures/evLevR_off2LR_' casesstrs{caseflag} '/'];
         if ~exist(most_recent_folder,'dir'); mkdir(most_recent_folder); mkdir([most_recent_folder 'span' num2str(Nspan) 'ms_lag' num2str(Nlag) 'ms']); end
-        %most_recent_folder='../../Figures/Temp/NoBaselineNorm/evLevR_off2R_del2LR_NOshift_NOoff2R/';
         selectedTimes='((ntce_ds(4,1)-40):ntce_ds(4,2))';
-        %selectedTimes='(ntce_ds(5,1):ntce_ds(6,1)-2)'; % FULL DELAY2
         %selectedCases{1}='find((~cciso2LL) & (ccepx(:,tt)<0))';
         %selectedCases{2}='find((~cciso2LL) & (ccepx(:,tt)>0))';
 end
 
-% for off1L del1L/R
-% for off1R del1L/R
-% for off2L del2L/R
-% for off2R del2L/R
-
-
 if reload
 % if ~exist('subjsData','var')
 %     subjsData(1).sbdata=getSubjectData(1);
@@ -82,12 +66,15 @@ if reload
 
 if ~exist('poolstackvars','var')
     load latestpoolstackvars_14feb22.mat
+	% or getPoolStackVars.m 
 end
 if ~exist('pooltweyepos','var')
     load latestpooltweyepos_14feb22.mat
+    % or getTimewinsPoolEyePos.m 
 end
 if ~exist('pooltwspksq','var')
     load latestpooltwspks_14feb22.mat
+    % or getTimewinsPoolSpikeSeqs.m
 end
 
 twLabels={'preoffer1','offer1','delay1','offer2','delay2','refixate','choice-go','ch-hold'};
@@ -122,7 +109,6 @@ poolcovhist=cat(2,pooltwspksq(1).covhist,pooltwspksq(2).covhist,pooltwspksq(3).c
     pooltwspksq(5).covhist,pooltwspksq(6).covhist,pooltwspksq(7).covhist,pooltwspksq(8).covhist);
 %clearvars pooltwspksq pooltweyepos
 
-%fchg=find(poolcovhist>=.8); fchl=find(poolcovhist<.8);
 fchg=find(poolcovhist>=.999); fchl=find(poolcovhist<.999);
 ntcv=[fchg(diff([fchg sum(Ntps12)])>1)' fchl([diff(fchl) sum(Ntps12)]>1)'];
 ntce=[[1; 401; ntcv(1:end-1,2)] [400; ntcv(:,1)]];
@@ -145,10 +131,8 @@ end
 
 clearvars poolspksqmat pooleposxmat
 
-
 %%
 
-
 poolcovhist_ds=nanmovmean(poolcovhist,Nlag,Nlag);
 %fcgh_ds=find(poolcovhist_ds>=.8); fchl_ds=find(poolcovhist_ds<.8);
 fcgh_ds=find(poolcovhist_ds>=.999); fchl_ds=find(poolcovhist_ds<.999);
@@ -189,24 +173,9 @@ rtmes_ds=arrayfun(@(jj) length(ntce_ds(jj,1):ntce_ds(jj,2)), 1:8);
 poolspksqmat_ds_bln=poolspksqmat_ds;
 else
     load pool_all_light.mat
+    % save time at run-time
 end
 
-%%
-% fullfig();
-% subplot(1,2,1);
-% spy(vectorisen(~isnan(pooleposxmat_ds_bln(1,cumsum(nccs),:)),1:2));
-% subplot(1,2,2);
-% spy(vectorisen(~isnan(poolspksqmat_ds_bln(1,cumsum(nccs),:)),1:2));
-% 
-% fulloff2time=ntce_ds(5,1):ntce_ds(6,1)-20-2;
-% 
-% fullfig();
-% subplot(1,2,1);
-% spy(vectorisen(~isnan(pooleposxmat_ds_bln(1,cumsum(nccs),fulloff2time)),1:2));
-% subplot(1,2,2);
-% spy(vectorisen(~isnan(poolspksqmat_ds_bln(1,cumsum(nccs),fulloff2time)),1:2));
-
-
 %%
 off1time=ntce_ds(2,2)-(19:-1:0);
 %del1time=ntce_ds(3,2)-(19:-1:0);
@@ -217,8 +186,8 @@ off2time=ntce_ds(4,2)-(19:-1:0);
 %fulldel1time=[(ntce_ds(3,1)-46):(ntce_ds(3,1)-7) ntce_ds(3,1):ntce_ds(3,2)];
 %fulloff1time=[(ntce_ds(2,1)-46):(ntce_ds(2,1)-7) ntce_ds(2,1):ntce_ds(2,2)];
 fulloff1time=[(ntce_ds(2,1)-40):ntce_ds(2,2)];
-fulloff1time=fulloff1time(fulloff1time>0);
-%fulldel2time=(ntce_ds(5,1)):ntce_ds(6,1)-20-2; %due to spikecount/subsampling setup
+fulloff1time=fulloff1time(fulloff1time>0); % ensure it is not negative
+%fulldel2time=(ntce_ds(5,1)):ntce_ds(6,1)-20-2; % due to spikecount/subsampling setup
 %fulldel2time=[(ntce_ds(5,1)-46):(ntce_ds(5,1)-7) ntce_ds(5,1):ntce_ds(6,1)-20-2];
 fulloff2time=[(ntce_ds(4,1)-46):(ntce_ds(4,1)-7) ntce_ds(4,1):ntce_ds(4,2)];
 %fulloff2time=[(ntce_ds(4,1)-40):ntce_ds(4,2)];
@@ -337,7 +306,7 @@ else
     end
 end
 
-%%
+%% EXTRA CHECKS: plotting percentiles of eye tracks
 pooleposxmat_L=nan(Ntrs12,Nccs,Ntps_ds); pooleposxmat_R=nan(Ntrs12,Nccs,Ntps_ds);
 if any(caseflag==[1 2 4 5])
     pooleposxmat_L=pooleposxmat_shLalgn(:,:,:);