sim_monomolecular {epifitter}

Simulate an epidemic using the Monomolecular model

Description

Simulate a stochastic epidemic curve using the Monomolecular model.

Usage

sim_monomolecular(N = 10,dt = 1, y0 = 0.01, r,K = 1, 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

K

Maximum asymptote

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_monomolecular(N = 30, y0 = 0.01,dt = 5, r = 0.3, K = 1, alpha = 0.5, n = 4)