Skip to contents

Simulate daily double censored PMF

Source: R/model-module-helpers.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.

See also

Helper functions for model modules add_max_observed_delay(), add_pmfs(), convolution_matrix(), enw_reference_by_report(), enw_reps_with_complete_refs(), extract_obs_metadata(), extract_sparse_matrix(), latest_obs_as_matrix()

Examples

simulate_double_censored_pmf(10, meanlog = 0, sdlog = 1)
#> Warning: `simulate_double_censored_pmf()` was deprecated in epinowcast 0.5.0.
#> ! Users should instead use the primarycensored package for simulating double
#>   censored processes.
#> ℹ See <https://github.com/epinowcast/primarycensored> for more information.
#>  [1] 0.238044 0.409753 0.168789 0.076520 0.039802 0.022420 0.013734 0.008822
#>  [9] 0.005843 0.004064 0.002851