sir: Euler scheme for the SIR model

In shill1729/odeSolveR: Implementation of simple Euler schemes and fourth-order Runge-Kutta methods for arbitrary first-order ODE systems

Description

A basic Euler scheme is used to numerically solve the ODE system of the standard SIR model.

Usage

sir(
  beta,
  gamma,
  s0 = 1,
  i0 = 1e-04,
  r0 = 0,
  t0 = 0,
  tn = 1,
  n = 1000,
  verbose = TRUE
)

Arguments

beta	the mean rate of contact, 1/beta the mean time between contact for a given individual (homogeneous population and mixing)
gamma	the mean rate of recovery, 1/gamma the mean time until recovery.
s0	the initial susceptible population
i0	the initial infected population
r0	the initial recovering population, defaults to zer0
t0	initial time point
tn	the terminal time point
n	number of sub-intervals in time-grid
verbose	whether to plot the populations over time and print parameter+reproductive ratio

Value

list

• parameters, a data.frame containing the contact rate, recovery rate, and r0 the basic reproduction ratio.

• sir, a matrix of the three populations in columns over time, susceptible, infected, recovered

Examples

{sir.dat <- sir(beta = 1/2, gamma = 1/3, s0 = 50, i0 = 1, tn = 90, n = 1000)}

shill1729/odeSolveR documentation built on March 31, 2021, 10:52 a.m.