Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/epiverse-trace/epichains
Methods for simulating and analysing the sizes and lengths of infectious disease transmission chains from branching process models
https://github.com/epiverse-trace/epichains
branching-processes epidemic-dynamics epidemic-modelling epidemic-simulations epidemiology epidemiology-models outbreak-simulator r-package r-stats transmission-chain transmission-chain-reconstruction
Last synced: about 1 month ago
JSON representation
Methods for simulating and analysing the sizes and lengths of infectious disease transmission chains from branching process models
- Host: GitHub
- URL: https://github.com/epiverse-trace/epichains
- Owner: epiverse-trace
- License: other
- Created: 2023-05-02T21:12:43.000Z (over 1 year ago)
- Default Branch: main
- Last Pushed: 2024-11-19T16:33:46.000Z (about 1 month ago)
- Last Synced: 2024-12-01T00:50:11.281Z (about 1 month ago)
- Topics: branching-processes, epidemic-dynamics, epidemic-modelling, epidemic-simulations, epidemiology, epidemiology-models, outbreak-simulator, r-package, r-stats, transmission-chain, transmission-chain-reconstruction
- Language: R
- Homepage: https://epiverse-trace.github.io/epichains/
- Size: 9.91 MB
- Stars: 7
- Watchers: 5
- Forks: 2
- Open Issues: 27
-
Metadata Files:
- Readme: README.Rmd
- License: LICENSE
- Citation: CITATION.cff
Awesome Lists containing this project
README
---
output: github_document
bibliography: vignettes/references.json
link-citations: true
---
```{r setup, include = FALSE}
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.path = file.path("man", "figures", "README-"),
out.width = "100%",
echo = TRUE
)
```
# _{{ packagename }}_: Methods for simulating and analysing the size and length of transmission chains from branching process models 
[](https://opensource.org/licenses/MIT)
[](https://github.com/epiverse-trace/epichains/actions/workflows/R-CMD-check.yaml)
[](https://app.codecov.io/gh/{{ gh_repo }}?branch=main)
[](https://lifecycle.r-lib.org/articles/stages.html#experimental)
[](https://CRAN.R-project.org/package=epichains)
_{{ packagename }}_ is an R package to simulate, analyse, and visualize the size
and length of branching processes with a given offspring distribution. These
models are often used in infectious disease epidemiology, where the chains represent chains of
transmission, and the offspring distribution represents the distribution of
secondary infections caused by an infected individual.
_{{ packagename }}_ re-implements [bpmodels](https://github.com/epiforecasts/bpmodels/)
by providing bespoke functions and data structures that allow easy
manipulation and interoperability with other Epiverse-TRACE packages, for example, [superspreading](https://github.com/epiverse-trace/superspreading/) and [epiparameter](https://github.com/epiverse-trace/epiparameter/), and potentially some existing packages for handling transmission chains, for example, [epicontacts](https://github.com/reconhub/epicontacts).
_{{ packagename }}_ is developed at the [Centre for the Mathematical Modelling of Infectious Diseases](https://www.lshtm.ac.uk/research/centres/centre-mathematical-modelling-infectious-diseases) at the London School of Hygiene and Tropical Medicine as part of the [Epiverse Initiative](https://data.org/initiatives/epiverse/).
## Installation
Install the released version of the package:
```{r install_cran, include=TRUE,eval=FALSE}
install.packages("{{ packagename }}")
```
The latest development version of the _{{ packagename }}_ package can be installed via
```{r install_with_remotes, include=TRUE,eval=FALSE}
# check whether {remotes} is installed
if (!require("remotes")) install.packages("remotes")
remotes::install_github("{{ gh_repo }}")
```
If this fails, try using the `pak` R package via
```{r install_with_pak, include=TRUE,eval=FALSE}
# check whether {pak} is installed
if (!require("pak")) install.packages("pak")
pak::pak("{{ gh_repo }}")
```
If both of these options fail, please [file an issue](https://github.com/epiverse-trace/epichains/issues) with a full log of the error messages. Here is an [example of an issue reporting an installation failure](https://github.com/epiverse-trace/epichains/issues/262). This will help us to improve the installation process.
To load the package, use
```{r load_pkg, eval=TRUE}
library("{{ packagename }}")
```
## Quick start
_{{ packagename }}_ provides three main functions:
* `simulate_chains()`: simulates transmission chains using a simple
branching process model that accepts an index number of cases that seed
the outbreak, a distribution of offspring per case, and a chain statistic
to track (size or length/duration). It optionally accepts other population
related inputs such as the population size (defaults to Inf) and percentage
of the population initially immune (defaults to 0). This function returns
an object with columns that track information on who infected whom, the
generation of infection and, if a generation time function is specified, the
time of infection.
* `simulate_chain_stats()`: provides a performant version of `simulate_chains()`
that only tracks and return a vector of realized chain sizes or
lengths/durations for each index case without details of the infection tree.
* `likelihood()`: calculates the loglikelihood (or likelihood, depending
on the value of `log`) of observing a vector of transmission chain sizes or
lengths.
The objects returned by the `simulate_*()` functions can be summarised with
`summary()`. Running `summary()` on the output of `simulate_chains()` will
return the same output as `simulate_chain_stats()` using the same inputs.
Objects returned from `simulate_chains()` can be aggregated into a
`` of cases per time or generation
with the function `aggregate()`.
The simulated `` object can be plotted in various ways using
`plot()`. See the plotting section in `vignette("epichains")` for two
use cases.
### Simulation
For the simulation functionality, let's look at a simple example where we
simulate a transmission chain with $20$ index cases, a constant generation time
of $3$, and a poisson offspring distribution with mean $1$. We are tracking the
chain "size" statistic and will cap all chain sizes at $25$ cases. We will then
look at the summary of the simulation, and aggregate it into cases per
generation.
```{r simulate_chains, eval=TRUE}
set.seed(32)
# Simulate chains
sim_chains <- simulate_chains(
n_chains = 20,
statistic = "size",
offspring_dist = rpois,
stat_threshold = 25,
generation_time = function(n) {
rep(3, n)
}, # constant generation time of 3
lambda = 1 # mean of the Poisson distribution
)
# View the head of the simulation
head(sim_chains)
# Summarise the simulation
summary(sim_chains)
# Aggregate the simulation into cases per generation
chains_agregegated <- aggregate(sim_chains, by = "generation")
# view the time series of cases per generation
chains_agregegated
```
### Inference
Let's look at the following example where we estimate the log-likelihood of
observing a hypothetical `chain_lengths` dataset.
```{r likelihood_example, eval=TRUE}
set.seed(32)
# randomly generate 20 chain lengths between 1 to 40
chain_lengths <- sample(1:40, 20, replace = TRUE)
chain_lengths
# estimate loglikelihood of the observed chain sizes
likelihood_eg <- likelihood(
chains = chain_lengths,
statistic = "length",
offspring_dist = rpois,
lambda = 0.99
)
# Print the estimate
likelihood_eg
```
Each of the listed functionalities is demonstrated in detail
in the ["Getting Started" vignette](https://epiverse-trace.github.io/epichains/articles/epichains.html).
## Package websites
The package has two websites: one for [the stable release version on CRAN](https://epiverse-trace.github.io/epichains/), and another for [the version in development](https://epiverse-trace.github.io/epichains/dev/).
## Package vignettes
The theory behind the models provided here can be
found in the [theory vignette](https://epiverse-trace.github.io/epichains/articles/theoretical_background.html).
We have also collated a bibliography of branching process applications in
epidemiology. These can be found in the [literature vignette](https://epiverse-trace.github.io/epichains/articles/branching_process_literature.html).
Specific use cases of _{{ packagename }}_ can be found in
the [online documentation as package vignettes](https://epiverse-trace.github.io/epichains/), under "Articles".
## Related R packages
As far as we know, below are the existing R packages for simulating branching
processes and transmission chains.
Click to expand
* [bpmodels](https://github.com/epiforecasts/bpmodels/): provides methods
for analysing the size and length of transmission chains from branching
process models. `{epichains}` supersedes `{bpmodels}`, which has been retired.
* [ringbp](https://github.com/epiforecasts/ringbp): a branching process
model, parameterised to the 2019-nCoV outbreak, and used to quantify
the potential effectiveness of contact tracing and isolation of cases.
* [covidhm](https://github.com/biouea/covidhm): code for simulating
COVID-19 dynamics in a range of scenarios across a real-world social
network. The model is conceptually based on `{ringbp}`.
* [epicontacts](https://github.com/reconhub/epicontacts): provides methods
for handling, analysing, and visualizing transmission chains and contact-tracing
data/linelists.
* [simulist](https://epiverse-trace.github.io/simulist/): uses a branching
process model to simulate individual-level infectious disease outbreak data,
including line lists and contact tracing data. This package is part of the
Epiverse-TRACE Initiative.
* [superspreading](https://epiverse-trace.github.io/superspreading/): provides
a set of functions to estimate and understand individual-level variation in
transmission of infectious diseases from data on secondary cases. These are
useful for understanding the role of superspreading in the spread of infectious
diseases and for informing public health interventions.
* [earlyR](https://github.com/reconhub/earlyR): estimates the reproduction
number (R), in the early stages of an outbreak. The model requires a
specified serial interval distribution, characterised by the mean and
standard deviation of the (Gamma) distribution, and data on daily disease
incidence, including only confirmed and probable cases.
* [projections](https://github.com/reconhub/projections): uses data on daily
incidence, the serial interval (time between onsets of infectors and
infectees) and the reproduction number to simulate plausible epidemic
trajectories and project future incidence. It relies on a branching process
where daily incidence follows a Poisson or a Negative Binomial distribution
governed by a force of infection.
* [simulacr](https://github.com/reconhub/simulacr): simulates outbreaks
for specified values of reproduction number, incubation period, duration
of infectiousness, and optionally reporting delays. Outputs a linelist
stored as a `data.frame` with the class `outbreak`, including information
on transmission chains; the output can be converted to ``
objects for visualisation.
* [outbreakr2](https://github.com/reconhub/outbreaker2): a Bayesian
framework for integrating epidemiological and genetic data to reconstruct
transmission trees of densely sampled outbreaks. It re-implements,
generalises and replaces the model of outbreaker, and uses a modular
approach which enables fine customisation of priors, likelihoods
and parameter movements.
* [o2geosocial](https://github.com/alxsrobert/o2geosocial): integrates
geographical and social contact data to reconstruct transmission chains.
It combines the age group, location, onset date and genotype of cases
to infer their import status, and their likely infector.
* [nosoi](https://github.com/slequime/nosoi): simulates agent-based
transmission chains by taking into account the influence of multiple
variables on the transmission process (e.g. dual-host systems
(such as arboviruses), within-host viral dynamics, transportation,
population structure), alone or taken together, to create complex
but relatively intuitive epidemiological simulations.
* [TransPhylo](https://xavierdidelot.github.io/TransPhylo/index.html):
reconstructs infectious disease transmission using genomic data.
## Reporting bugs
To report a bug please open an [issue](https://github.com/epiverse-trace/epichains/issues/new/choose).
## Contribute
Contributions to {epichains} are welcomed. Please follow the [package contributing guide](https://github.com/epiverse-trace/.github/blob/main/CONTRIBUTING.md).
## Code of conduct
Please note that the _{{ packagename }}_ project is released with a [Contributor Code of Conduct](https://github.com/epiverse-trace/.github/blob/main/CODE_OF_CONDUCT.md).
By contributing to this project, you agree to abide by its terms.
## Citing this package
```{r citation, message=FALSE, warning=FALSE}
citation("epichains")
```