Skip to contents

SIRD model

Source: R/ModelSIRD.R
ModelSIRD.Rd

SIRD model

Usage

ModelSIRD(name, prevalence, transmission_rate, recovery_rate, death_rate)

# S3 method for epiworld_sird
plot(x, main = get_name(x), ...)

Arguments

name

String. Name of the virus

prevalence

Double. Initial proportion of individuals with the virus.

transmission_rate

Numeric scalar between 0 and 1. Virus's rate of infection.

recovery_rate

Numeric scalar between 0 and 1. Rate of recovery_rate from virus.

death_rate

Numeric scalar between 0 and 1. Rate of death from virus.

x

Object of class SIR.

main

Title of the plot

...

Additional arguments passed to graphics::plot.

Value

  • The plot function returns a plot of the SIRD model of class epiworld_model.

Details

The initial_states function allows the user to set the initial state of the model. The user must provide a vector of proportions indicating the following values: (1) proportion of non-infected agents already removed, and (2) proportion of non-ifected agents already deceased.

See also

epiworld-methods

Other Models: ModelDiffNet(), ModelSEIR(), ModelSEIRCONN(), ModelSEIRD(), ModelSEIRDCONN(), ModelSIR(), ModelSIRCONN(), ModelSIRDCONN(), ModelSIRLogit(), ModelSIS(), ModelSISD(), ModelSURV(), epiworld-data

Examples

model_sird <- ModelSIRD(
 name              = "COVID-19",
 prevalence        = 0.01, 
 transmission_rate = 0.9,
 recovery_rate     = 0.1,
 death_rate        = 0.01
)

# Adding a small world population
agents_smallworld(
  model_sird,
  n = 1000,
  k = 5,
  d = FALSE,
  p = .01
  )
  
# Running and printing
run(model_sird, ndays = 100, seed = 1912)
#> _________________________________________________________________________
#> Running the model...
#> ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| done.
#>  done.
model_sird
#> ________________________________________________________________________________
#> Susceptible-Infected-Recovered-Deceased (SIRD)
#> It features 1000 agents, 1 virus(es), and 0 tool(s).
#> The model has 4 states.
#> The final distribution is: 0 Susceptible, 5 Infected, 912 Recovered, and 83 Deceased.

# Plotting
plot(model_sird)