SinePowerSpec now plots representative traces with calibrated photodiode trace, and fat ISI analysis outputs habituation for both on and off responses

Author: sammykatta

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

@@ -339,7 +339,32 @@
 stimTraces = [stimTraces{:}];
 tVec = ((1:length(stimTraces))/sf)';
 
-yyh = plotyy(tVec,respTraces,tVec,stimTraces);
+% Plot calibrated photodiode trace as alternative
+a = ephysData.FAT170.data{3,52}; % probe steps
+a = a-mean(a(1:1499));
+b = ephysData.FAT170.data{3,53}-ephysData.FAT170.data{3,54}; %pd steps minus worm only pd steps
+b = b-mean(b(1:1499));
+c = [0 2.5 5 7.5 10 12.5];
+
+% calibCells = {'FAT170'};
+% calibData = cell(0);
+% calibData{1} = b;
+% 
+% clear handles;
+% plotData = calibData{:,1};
+% handles = selectCalibSteps(plotData,calibCells{1});
+% stepIdx = handles.cursorPoints(:,1);
+% stepIdx = reshape(stepIdx,2,[])';
+% for j = 1:size(stepIdx,1)
+%     stepValues(1,j) = mean(plotData(stepIdx(j,1):stepIdx(j,2)))';
+% end
+% close;      
+% d = stepValues(1:6);
+d = [0.000719859832312041,-0.196562355775423,-0.460558335247858,-0.755470591227155,-1.11678911522346,-1.51308976348689];
+e = stimTraces;
+pdTraces = interp1(-d,c,-e,'pchip');
+
+yyh = plotyy(tVec,respTraces,tVec,pdTraces);
 hline(respMean(1),'b'); hline(respMean(2),'r'); hline(respMean(3),'m');
 
 %% Plot zoomed traces for 5 cycles
@@ -352,10 +377,10 @@
 yyh = cell(0);
 for i = 1:3
     figure();
-    yyh{i} = plotyy(tVec,respTraces(:,i),tVec,stimTraces(:,i));
+    yyh{i} = plotyy(tVec,respTraces(:,i),tVec,pdTraces(:,i));
     hline(respMean(i),'k:');
     set(yyh{i},'XLim',cycleLims(i,:));
-    set(yyh{i}(2),'YLim',[-1 0]);
+    set(yyh{i}(2),'YLim',[0 10]);
     set(yyh{i}(1),'Ylim',[respMean(i)-15e-12 respMean(i)+15e-12]);
 end
 

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

@@ -121,14 +121,28 @@
 
 %% Plot them
 
-figure(); hold on; ylim(gca,[0 1.75]);
+fh(1) = figure(); hold on; ylim(gca,[0 1.75]);
 plot(wtOnMean(1:40),'b','LineWidth',2);
 plot(fatOnMean(1:40),'Color',[0 1 0],'LineWidth',2);
+axh(1) = gca;
+fh(2) = figure(); hold on; ylim(gca,[0 1.75]);
+plot(wtOffMean(1:40),'b','LineWidth',2);
+plot(fatOffMean(1:40),'Color',[0 1 0],'LineWidth',2);
+axh(2) = gca;
 plotfixer;
-plot(wtOnNorm(1:40,:),'Color',[.7 .7 1])
-plot(fatOnNorm(1:40,:),'Color',[.2 .7 .2])
-ylabel('Normalized on current');
-xlabel('Stimulus #');
-chH = get(gca,'children');
-set(gca,'children',flipud(chH));
+
+plot(axh(1),wtOnNorm(1:40,:),'Color',[.7 .7 1])
+plot(axh(1),fatOnNorm(1:40,:),'Color',[.2 .7 .2])
+ylabel(axh(1),'Normalized on current');
+xlabel(axh(1),'Stimulus #');
+chH = get(axh(1),'children');
+set(axh(1),'children',flipud(chH));
+
+plot(axh(2),wtOnNorm(1:40,:),'Color',[.7 .7 1])
+plot(axh(2),fatOnNorm(1:40,:),'Color',[.2 .7 .2])
+ylabel(axh(2),'Normalized off current');
+xlabel(axh(2),'Stimulus #');
+chH = get(axh(2),'children');
+set(axh(2),'children',flipud(chH));
+