dust_example: Access dust's built-in examples

In mrc-ide/dust: Iterate Multiple Realisations of Stochastic Models

Access dust's built-in examples

Description

Access dust's built-in examples. These are compiled into the package so that examples and tests can be run more quickly without having to compile code directly via dust(). These examples are all "toy" examples, being small and fast to run.

Usage

dust_example(name)

Arguments

name	The name of the example to use. There are five examples: sir, sirs, variable, volatility, walk and logistic (see Details).

Details

• sir: a basic SIR (Susceptible, Infected, Resistant) epidemiological model. Draws from the binomial distribution to update the population between each time step.

• sirs: an SIRS model, the SIR model with an added R->S transition. This has a non-zero steady state, so can be run indefinitely for testing.

• volatility: A volatility model that might be applied to currency fluctuations etc.

• walk: A 1D random walk, following a Gaussian distribution each time step.

• logistic: Logistic growth in continuous time

Value

A dust_generator object that can be used to create a model. See examples for usage.

Examples

# A SIR (Susceptible, Infected, Resistant) epidemiological model
sir <- dust::dust_example("sir")
sir

# Initialise the model at time step 0 with 50 independent trajectories
mod <- sir$new(list(), 0, 50)

# Run the model for 400 steps, collecting "infected" every 4th time step
times <- seq(0, 400, by = 4)
mod$set_index(2L)
y <- mod$simulate(times)

# A plot of our epidemic
matplot(times, t(drop(y)), type = "l", lty = 1, col = "#00000044",
        las = 1, xlab = "Time", ylab = "Number infected")

mrc-ide/dust documentation built on Sept. 25, 2023, 10:09 p.m.