Skip to contents

Simulate daily double censored PMF. From epinowcast: https://package.epinowcast.org/dev/reference/simulate_double_censored_pmf.html #nolint

Source: R/delay_distribs.R
simulate_double_censored_pmf.Rd

This function simulates the probability mass function of a daily double-censored process. The process involves two distributions: a primary distribution which represents the censoring process for the primary event and another distribution (which is offset by the primary).

Usage

simulate_double_censored_pmf(
  max,
  fun_primary = stats::runif,
  primary_args = list(),
  fun_dist = stats::rlnorm,
  dist_args = list(...),
  n = 1e+06,
  ...
)

Arguments

max

Maximum value for the computed CDF. If not specified, the maximum value is the maximum simulated delay.

fun_primary

Primary distribution function (default is runif).

primary_args

List of additional arguments to be passed to the primary distribution function.

fun_dist

Distribution function to be added to the primary (default is rlnorm).

dist_args

List of additional arguments to be passed to the distribution function.

n

Number of simulations (default is 1e6).

...

Additional arguments to be passed to the distribution function. This is an alternative to dist_args.

Value

A numeric vector representing the PMF.

Details

Based off of: https://www.medrxiv.org/content/10.1101/2024.01.12.24301247v1