Update 'code/neuralAnalysesDelayShiftAligned.m'

code/neuralAnalysesDelayShiftAligned.m

@@ -9,10 +9,9 @@ clearvars -except subjsData poolstack* pooltw*;
 reload=0; recomp_shLRalgn=1;
 recomputelm=1;
 Nlag=10; Nspan=200; %Nspan=200;
-Nshf=10;
+Nshf=100;
 
-
-selflag=1; %LookL during offer1 is most common
+%selflag=1; %LookL during offer1 is most common
 %selflag=4; %LookR during offer2 is most common
 %selflag=2; %LookR during offer1 is less common
 %selflag=3; %LookL during offer2 is less common
@@ -23,7 +22,7 @@ selflag=1; %LookL during offer1 is most common
 %caseflag=4; % same as caseflag=1 but forcing offerLL/LR
 %caseflag=5; % same as caseflag=2 but forcing offerLL/LR
 %caseflag=6; % same as caseflag=3 but forcing offerLL/LR
-for selflag=4%[1 4]
+for selflag=[1 4]
 for caseflag=2
 clearvars -except subjsData reload recomp_shLRalgn recomputelm Nlag Nspan Nshf selflag caseflag; close all; clc;
 
@@ -35,43 +34,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_apr24/evLevR_off1L_del1LR_' casesstrs{caseflag} '/'];
+        most_recent_folder=['../../Figures/evLevR_off1L_del1LR_' 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(3,1)-46):ntce_ds(3,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 
+        most_recent_folder=['../../Figures/evLevR_off1R_del1LR_' casesstrs{caseflag} '/'];
+        selectedTimes='((ntce_ds(3,1)-46):ntce_ds(3,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
+        most_recent_folder=['../../Figures/evLevR_off2L_del2LR_' casesstrs{caseflag} '/'];
+        selectedTimes='((ntce_ds(5,1)-46):ntce_ds(5,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_apr24/evLevR_off2R_del2LR_' casesstrs{caseflag} '/'];
+        most_recent_folder=['../../Figures/evLevR_off2R_del2LR_' 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(5,1)-46):ntce_ds(5,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);
@@ -120,7 +105,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)]];
@@ -146,7 +130,6 @@ 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);
 ntcv_ds=[fcgh_ds(diff([fcgh_ds Ntps_ds])>1)' fchl_ds([diff(fchl_ds) Ntps_ds]>1)'];
 vtps_ds=ones(1,Ntps_ds);  for jj=1:7; vtps_ds(ntcv_ds(jj,1):ntcv_ds(jj,2))=nan; end
@@ -185,24 +168,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));
-% 
-% fulldel2time=ntce_ds(5,1):ntce_ds(6,1)-20-2;
-% 
-% fullfig();
-% subplot(1,2,1);
-% spy(vectorisen(~isnan(pooleposxmat_ds_bln(1,cumsum(nccs),fulldel2time)),1:2));
-% subplot(1,2,2);
-% spy(vectorisen(~isnan(poolspksqmat_ds_bln(1,cumsum(nccs),fulldel2time)),1:2));
-
-
 %%
 off1time=ntce_ds(2,2)-(19:-1:0);
 del1time=ntce_ds(3,2)-(19:-1:0);
@@ -211,7 +179,7 @@ del2time=ntce_ds(5,2)-(19:-1:0);
 
 %fulldel1time=(ntce_ds(3,1)):ntce_ds(3,2);
 fulldel1time=[(ntce_ds(3,1)-46):(ntce_ds(3,1)-7) ntce_ds(3,1):ntce_ds(3,2)];
-%fulldel2time=(ntce_ds(5,1)):ntce_ds(6,1)-20-2; %due to spikecount/subsampling setup
+%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];
 
 num1stshiftsL=zeros(248,1);