load('Fig5e.mat')

ASAP3firing = permute(ASAP3firing, [2,3,1]);
SpikeyGi2firing = permute(SpikeyGi2firing, [2,3,1]);
SpikeyGifiring = permute(SpikeyGifiring, [2,3,1]);
x=[1:50]/10;

figure
temp = nanmean(ASAP3firing(1,:,:),3);
stemp = std(ASAP3firing(1,:,:),[],3)/sqrt(size(ASAP3firing,3)/3);
confplot(x,temp,stemp,stemp);
hold on
temp = nanmean(SpikeyGifiring(1,:,:),3);
stemp = std(SpikeyGifiring(1,:,:),[],3)/sqrt(size(SpikeyGifiring,3)/3);
confplot(x,temp,stemp,stemp);
hold on
temp = nanmean(SpikeyGi2firing(1,:,:),3);
stemp = std(SpikeyGi2firing(1,:,:),[],3)/sqrt(size(SpikeyGi2firing,3)/3);
confplot(x,temp,stemp,stemp);
xlim([3.5 5])
title('spontaneous activity')
xlabel('SNR')
ylabel('firing rate (Hz)')

figure; 
temp = nanmean(ASAP3firing(2,:,:),3);
stemp = std(ASAP3firing(2,:,:),[],3)/sqrt(size(ASAP3firing,3)/3);
confplot(x,temp,stemp,stemp);
hold on
temp = nanmean(SpikeyGifiring(2,:,:),3);
stemp = std(SpikeyGifiring(2,:,:),[],3)/sqrt(size(SpikeyGifiring,3)/3);
confplot(x,temp,stemp,stemp);
hold on
temp = nanmean(SpikeyGi2firing(2,:,:),3);
stemp = std(SpikeyGi2firing(2,:,:),[],3)/sqrt(size(SpikeyGi2firing,3)/3);
confplot(x,temp,stemp,stemp);
xlim([3.5 5])
title('evoked activity')
xlabel('SNR')
ylabel('firing rate (Hz)')