SIRD model

## 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

`ModelSIRD`

function returns a model of class epiworld_model.

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)
```