VoltageNMethods2023/Figure 5/SpikeyGi_10Hz.m

load('SpikeyGi_10Hz');
figure
hold on
t = [1:size(temp4,2)];
t = t-1600;
t = t/803.3179;

for i = 1:5
    if i == 1
        rrr = nanmean(temp4(:,1630:1640),2);
        sss = nanmax(temp6(:,1630:1640),[],2);
    elseif i == 2
        rrr = nanmean(temp4(:,1702:1712),2);
        sss = nanmax(temp6(:,1702:1712),[],2);
    elseif i == 3
        rrr = nanmean(temp4(:,1780:1790),2);
        sss = nanmax(temp6(:,1780:1790),[],2);
    elseif i == 4
        rrr = nanmean(temp4(:,1865:1875),2);
        sss = nanmax(temp6(:,1865:1875),[],2);
    elseif i == 5
        rrr = nanmean(temp4(:,1955:1965),2);
        sss = nanmax(temp6(:,1955:1965),[],2);
    end
    ccc=temp4(sss>2&rrr>0.2,1200:1400);
    jjj = temp4(sss>2&rrr>0.2,:);
    confplot(t,nanmean(temp4(sss>2&rrr>0.2,:),1)-nanmean(ccc(:)),nanstd(jjj,1)/sqrt(size(jjj,1)), nanstd(jjj,1)/sqrt(size(jjj,1))) ;
    hold on
    xline(0)
    xline((onepuff(3,2)-onepuff(2,2))/803.3179);
    xline((onepuff(4,2)-onepuff(2,2))/803.3179);
    xline((onepuff(5,2)-onepuff(2,2))/803.3179);
    xline((onepuff(6,2)-onepuff(2,2))/803.3179);
end

xlim([-0.5 1.5])
yline(0)
title('iASAP 10Hz')
ylabel('dFF')
xlabel('seconds')
ylim([-0.05 0.2])