- Author:
- schdemp <sdem348@aucklanduni.ac.nz>
- Date:
- 2022-06-14 17:43:11-04:00
- Desc:
- Added Matlab codes and renamed input param filenames for ease of use.
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/8a2/rawfile/6905429673d8289e3908cff6db2149f51894c44d/MatlabFiles/Figure3.m
f=figure(3);
tiledlayout(5,1);
DataSets={'Fig1A','Fig1B','Fig1C','Fig1D','Fig1E'};
TitleName={'Sensory stimulation (awake)','GABAergic interneurons OG (awake)',...
'Sensory Stimulation (anesthesia)','GABAergic interneurons OG (anesthesia)',...
'Pyramidal neurons OG (anesthesia)'};
PlotFlags=[2000,2000,2500,4500,4500];
for i=1:length(DataSets)
nexttile%subplot(5,1,i)
T=readtable(strcat('sten2020_NVC_main_ISO_data_',DataSets{i},'.csv'));
%Time is Column 1, Signal is Column 2, Dilation variable is column 13
PlotData=table2array(T(:,[1,2,13]));
plot(PlotData(1:PlotFlags(i),1),PlotData(1:PlotFlags(i),3),'r-')
hold on
plot(PlotData(1:PlotFlags(i),1),PlotData(1:PlotFlags(i),2)-10,'k-')
ylabel('\Delta Dilation (%)','FontSize',10)
title(TitleName{i},'FontSize',10,'FontWeight','Bold')
ylim([-11 20]);
xlim([0 45])
yticks(-5:5:15)
xticks(0:10:40)
if i==1
legset=legend('Best Model Simulation','Neuronal Stim');
end
if i==5
xlabel('Time (Seconds)','FontSize',10)
end
end
set(gcf,'Position',[100 100 400 700])
exportgraphics(f,'Figure3.png','Resolution',600)
