sim_gompertz: Simulate an epidemic using the Gompertz model

In epifitter: Analysis and Simulation of Plant Disease Progress Curves

Description

Simulate a stochastic epidemic curve using the Gompertz model.

Usage

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

epifitter documentation built on June 14, 2021, 5:08 p.m.