% PLOT_SEIR.m % % analyze and plot the model results LW=2;FS=16;FW='demi'; % Figure 1: 2 subplots S, I figure subplot(2,2,1) plot(t,S,'b-','LineWidth',LW) title('Susceptible','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Individuals m^{-2}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) subplot(2,2,2) plot(t,E,'g-','LineWidth',LW) title('Exposed','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Individuals m^{-2}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) subplot(2,2,3) plot(t,I,'r-','LineWidth',LW) title('Infected','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Individuals m^{-2}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) subplot(2,2,4) plot(t,R,'k-','LineWidth',LW) title('Recovered','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Individuals m^{-2}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) %xlim([0 100]) % print and sent plot to a folder % print('-dpng',['/Path/CAse1.png']) % Figure 2: Infection rate by contact figure N=2199711; infI=PAR.Iinfect .* I .* S; plot(t,infI,'r-','LineWidth',LW) title('Infection rate','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Infection rate (individuals m^{-2} day^{-1})','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) % print('-dpng',['Figures/Ab' PAR.case 'InfRate.png']) % Figure 3 for Susceptible, Infected, Dead together figure plot(t,S,'b-',t,E,'g-',t,I,'r-',t,R,'k-','LineWidth',LW) %title(' Susceptible, Infectd','FontSize',FS,'FontWeight',FW) legend({'Susceptible','Exposed','Infected','Recovered'},'Location','northeast','FontSize',16) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Individuals m^{-2}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) % print('-dpng',['Figures/Ab' PAR.case 'InfRate.png']) %xlim([0 100])