% PLOT_SIPDB.m % % Analyze and plot the model results % Fonts, selected options LW=2;FS=16;FW='demi'; % Figure 1: 4 subplots S, I, D, 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,D,'k-','LineWidth',LW) title('Dead','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,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) %xlim([0 100]) % print and sent plot to a folder % print('-dpng',['/Path/CAse1.png']) % Figure 2: Infection rate by contact with P, D and I figure N=2199711; infP=PAR.Pinfect .* P .* S; infI=PAR.Iinfect .* P .* S; infD=PAR.Dinfect .* P .* S; plot(t,infI,'r-',t,infD,'b',t,infP,'k-.','LineWidth',LW) title('Infection rate','FontSize',FS,'FontWeight',FW) legend({'Contact with I','Contact with D' 'Contact with P'},'Location','northeast','FontSize',FS) 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 Susceptibel, Infected, Dead together figure plot(t,S,'b-',t,I,'r-',t,D,'k-','LineWidth',LW) %title(' Susceptible,Infected, Dead','FontSize',FS,'FontWeight',FW) legend({'Susceptible','Infected', 'Dead'},'Location','northeast','FontSize',FS) 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])