Intro_Biocomp_ND_317_Tutorial10 Code used to generate different parameters of beta and gamma modelSim=ode(y=initial,times=times,func=SIR,parms=c(beta, gamma))
Sim = data.frame(modelSim)
ggplot(Sim, aes(Sim$Time)) +
geom_line(aes(y=Sim1$Susceptible, color="mediumpurple3")) +
geom_line(aes(y=Sim1$Infected, color="mediumvioletred")) +
geom_line(aes(y=Sim1$Recovered, color="mediumseagreen")) +
theme_classic() +
xlab("Time (days)") +
ylab("Number of Individuals") +
scale_color_manual(labels=c("Susceptible", "Infected", "Recovered"),
values = c("mediumpurple3", "mediumvioletred", "mediumseagreen")) +
ggtitle("Simulation Plot") +
theme(plot.title = element_text(hjust = 0.5))
Increase in beta leads to faster increase to threshold capacity Decrease in beta leads to slower rate of increase to threshold capacity Beta is responsible for rate of change of infected and recovered Examples seen in simulation comparisons of parameters 1 and 6 Increase in gamma leads to lower threshold capacity of infected and recovered Decrease in gamma leads to lower threshold capacity of infected and recovered Gamma is responsible for carrying capacity of they system Examples seen in simulation comparisons of parameters 5 and 7