plotting macro currents

Author: sammykatta

Analysis/SK/Scripts-in-progress/Sim_Distance_Comparison.m

@@ -9,9 +9,10 @@
 simPath = 'C:\Users\Sammy\Dropbox\Goodman Lab\Posters Papers Proposals\2018-Katta-SpatiotemporalDynamics\source-data\sim-channel-distance\';
 simFiles = {'channel-distance-displacement.xlsx','channel-distance-speed.xlsx','channel-distance-frequency.xlsx'};
 simData = cell(0);
-for i = 1:3
-    [~,~,simData{i}] = xlsread(fullfile(simPath,simFiles{i}));
-    simData{i} = 
+simDataMat = cell(0);
+for iPanel = 1:3
+    [~,~,simData{iPanel}] = xlsread(fullfile(simPath,simFiles{iPanel}));
+    simDataMat{iPanel} = cell2mat(simData{iPanel}(2:end,:));
 end
 
 recData = cell(0);
@@ -28,8 +29,8 @@
 % figure('Position',[200 200 800 600]);
 
 for iPanel = 1:3
-    fh(iPanel,1) = figure('Position',[500 100 800 600]);
-    fh(iPanel,2) = figure('Position',[100 100 600 600]);
+    fh(iPanel,1) = figure('Position',[100 100 600 600]); % stim and macro currents
+    fh(iPanel,2) = figure('Position',[500 100 800 600]); % channel currents
 
     switch iPanel
         case 1 % displacement
@@ -45,34 +46,82 @@
 
     % find all subsets of columns now, use overlap in indices later to pull
     % out more specific subsets
-    stimCols = find(~cellfun(@isempty,cellfun(@(x) regexp(x,'^stim'),simHeaders,'un',0)));
-    currTotCols = find(~cellfun(@isempty,cellfun(@(x) regexp(x,'^currTot'),simHeaders,'un',0)));
-    chCurrCols = find(~cellfun(@isempty,cellfun(@(x) regexp(x,'^chcurr'),simHeaders,'un',0)));
-    onCols{1} = find(~cellfun(@isempty,cellfun(@(x) regexp(x,'_on'),simHeaders,'un',0)));
-    onCols{2} = find(~cellfun(@isempty,cellfun(@(x) regexp(x,'_off'),simHeaders,'un',0)));
+    stimCols = ~cellfun(@isempty,cellfun(@(x) regexp(x,'^stim'),simHeaders,'un',0));
+    currTotCols = ~cellfun(@isempty,cellfun(@(x) regexp(x,'^currTot'),simHeaders,'un',0));
+    chCurrCols = ~cellfun(@isempty,cellfun(@(x) regexp(x,'^chcurr'),simHeaders,'un',0));
+    onCols{1} = ~cellfun(@isempty,cellfun(@(x) regexp(x,'_on'),simHeaders,'un',0));
+    onCols{2} = ~cellfun(@isempty,cellfun(@(x) regexp(x,'_off'),simHeaders,'un',0));
 
     intTypes = cellfun(@(x) regexp(x,sprintf('_%s(\\S{1,5})',stType),'tokens'),simHeaders,'un',0);
     intTypes = [intTypes{:}];intTypes = [intTypes{:}];
     intTypes = unique(intTypes,'stable');
     for jIntensity = 1:length(intTypes);
-        intCols{jIntensity} = find(~cellfun(@isempty,cellfun(@(x) regexp(x,sprintf('_%s%s$',stType,intTypes{jIntensity})),simHeaders,'un',0)));
+        intCols{jIntensity} = ~cellfun(@isempty,cellfun(@(x) regexp(x,sprintf('_%s%s$',stType,intTypes{jIntensity})),simHeaders,'un',0));
     end    
 
     distTypes = cellfun(@(x) regexp(x,'_dist(\S{1,4})_','tokens'),simHeaders,'un',0);
     distTypes = [distTypes{:}];distTypes = [distTypes{:}];
     distTypes = unique(distTypes,'stable');
     for kDist = 1:length(distTypes);
-        distCols{kDist} = find(~cellfun(@isempty,cellfun(@(x) regexp(x,sprintf('_dist%s_',distTypes{kDist})),simHeaders,'un',0)));
-    end    
-
-    % NEXT: make subplots and start matching subsets for each subplot
+        distCols{kDist} = ~cellfun(@isempty,cellfun(@(x) regexp(x,sprintf('_dist%s_',distTypes{kDist})),simHeaders,'un',0));
+    end
+    
+    % make subplots, match column subsets for each subplot, and plot
     for jIntensity = 1:length(intTypes)
+        figure(fh(iPanel,1));
 
+        if strcmp(stType,'disp') || strcmp(stType,'speed') % plot on and off
+            % plot stimulus traces
+            axh1(jIntensity,1) = subplot(3,3,1+3*(jIntensity-1));
+            thisCol = sum([stimCols; intCols{jIntensity}; onCols{1}],1)==3;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'b'); % on = blue
+            hold on;
+            thisCol = sum([stimCols; intCols{jIntensity}; onCols{2}],1)==3;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'m'); % off = magenta
+            
+            % plot experimental total current
+            axh1(jIntensity,2) = subplot(3,3,2+3*(jIntensity-1));
+                % separate exp traces into aligned on and off for plotting
+            
+            % plot simulated total current
+            axh1(jIntensity,3) = subplot(3,3,3+3*(jIntensity-1));
+            thisCol = sum([currTotCols; intCols{jIntensity}; onCols{1}],1)==3;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'b');
+            hold on;
+            thisCol = sum([currTotCols; intCols{jIntensity}; onCols{2}],1)==3;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'m');
+
+        elseif strcmp(stType,'freq') % only one current to plot
+            % plot stimulus traces
+            axh1(jIntensity,1) = subplot(3,3,1+3*(jIntensity-1));
+            thisCol = sum([stimCols; intCols{jIntensity}],1)==2;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'b');
+            
+            % plot experimental total current
+            axh1(jIntensity,2) = subplot(3,3,2+3*(jIntensity-1));
+                % separate exp traces into aligned on and off for plotting
+            
+            % plot simulated total current
+            axh1(jIntensity,3) = subplot(3,3,3+3*(jIntensity-1));
+            thisCol = sum([currTotCols; intCols{jIntensity}],1)==2;
+            plot(simTime, simDataMat{iPanel}(:,thisCol),'b');
+        end
+
         for kDist = 1:length(distTypes)
+            figure(fh(iPanel,2));
             plotNo = kDist+4*(jIntensity-1);
             axh(plotNo) = subplot(3,4,plotNo);
+            if strcmp(stType,'disp') || strcmp(stType,'speed')
+                %plot on
+                %plot off
+            elseif strcmp(stType,'freq')
+                %plot current
+            end
+            
         end
 
+        % set xlim and ylim based on iPanel and jIntensity
+        
     end
 
     clear axh;