% PLOT_SIPF.m % % analyze and plot the model results LW=2;FS=16;FW='demi'; % Figure 1: 3 subplots S, I, P 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,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,3) plot(t,P,'g-','LineWidth',LW) title('Infectious Particles','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Pathogens m^{-3}','FontSize',FS,'FontWeight',FW) set(gca,'LineWidth',LW,'FontSize',FS,'FontWeight',FW) subplot(2,2,4) plot(t,F./S,'g-','LineWidth',LW) title('Infectious Particles','FontSize',FS,'FontWeight',FW) xlabel('Time (days)','FontSize',FS,'FontWeight',FW) ylabel('Pathogens S^{-1}','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.Finfect .* F; 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']) 3 % Figure 3 for Susceptibel, Infected, Dead together figure plot(t,S,'b-',t,I,'r-','LineWidth',LW) %title(' Susceptible, Infected',FontSize',FS,'FontWeight',FW) legend({'Susceptible','Infected'},'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])