sir-model

About This Project

Test SIR model based on Ottar's book book code from ch 1 births and deaths included

this uses fractional populations

$$\begin{align} \frac{dS}{dt} &= \mu (N - S) -\beta S \frac{I}{N} \\\ \frac{dI}{dt} &= \beta S \frac{I}{N} - \gamma I - \mu I \\\ \frac{dR}{dt} &= \gamma I - \mu R \end{align}$$ $$\begin{align} \mu &= \frac{1}{50*52} \\\ \beta &= 2 \\\ \gamma &= \frac{1}{2} \\\\\ N &= 1000 \\\ S_0 &= 999.0 \\\ I_0 &= 1.0 \\\ R_0 &= 0.0 \end{align}$$

Repository Structure

.
├── data/
├── figs/
├── funs/
├── out/
└── src/

• data/ contains ...
• figs/ contains ...
• funs/ contains ...
• out/ contains ...
• src/ contains ...

Built With

• R version 4.2.1 (2022-06-23)
• {renv} for package management
• {deSolve} for solving differential equations

Getting Started

I've used the {renv} package to manage the R environment for this project. For more details on how to use {renv}, see this article, but in brief, it creates a snapshot of the installed packages and their versions.

To get started, you will need to install {renv} as usual (i.e., install.packages("renv")), and then run renv::restore() to install the packages that are used in this project (the record in the renv.lock file).

Usage

License

This project is licensed under the MIT License - see the LICENSE file for details

Contact

nita at psu dot edu

Acknowledgements

Ottar Bjornstadt book