sim_exponential {epifitter}

Simulate an epidemic using the Exponential model

Description

Simulate a stochastic epidemic curve using the Exponential model.

Usage

sim_exponential(N = 10,dt = 1, y0 = 0.01, r, n,  alpha = 0.2)

Arguments

N

Total time course of the epidemic

dt

Time step

y0

Initial inoculum or initial disease intensity

r

Infection rate

n

Number or replicates or sample size for each time step

alpha

Variation parameter. stands for the variation for the replicates for each time step. The standard deviation is calculated as sd = alpha * y * (1 - y), being y the disease intensity for each time step.

Value

rep

Replicates

time

Time after epidemic start

y

Disease intensity

random_y

Disease intensity after applying the random alpha error

Examples

sim_exponential(N = 30, y0 = 0.01,dt = 5, r = 0.1, alpha = 0.5, n = 4)