Returns a list of options defining the secondary model used in
estimate_secondary()
. This model is a combination of a convolution of
previously observed primary reports combined with current primary reports
(either additive or subtractive). It can optionally be cumulative. See the
documentation of type
for sensible options to cover most use cases and the
returned values of secondary_opts()
for all currently supported options.
Usage
secondary_opts(type = c("incidence", "prevalence"), ...)
Arguments
- type
A character string indicating the type of observation the secondary reports are. Options include:
"incidence": Assumes that secondary reports equal a convolution of previously observed primary reported cases. An example application is deaths from an infectious disease predicted by reported cases of that disease (or estimated infections).
"prevalence": Assumes that secondary reports are cumulative and are defined by currently observed primary reports minus a convolution of secondary reports. An example application is hospital bed usage predicted by hospital admissions.
- ...
Overwrite options defined by type. See the returned values for all options that can be passed.
Value
A <secondary_opts>
object of binary options summarising secondary
model used in estimate_secondary()
. Options returned are cumulative
(should the secondary report be cumulative), historic
(should a
convolution of primary reported cases be used to predict secondary reported
cases), primary_hist_additive
(should the historic convolution of primary
reported cases be additive or subtractive), current
(should currently
observed primary reported cases contribute to current secondary reported
cases), primary_current_additive
(should current primary reported cases be
additive or subtractive).
Examples
# incidence model
secondary_opts("incidence")
#> $cumulative
#> [1] 0
#>
#> $historic
#> [1] 1
#>
#> $primary_hist_additive
#> [1] 1
#>
#> $current
#> [1] 0
#>
#> $primary_current_additive
#> [1] 0
#>
#> attr(,"class")
#> [1] "secondary_opts" "list"
# prevalence model
secondary_opts("prevalence")
#> $cumulative
#> [1] 1
#>
#> $historic
#> [1] 1
#>
#> $primary_hist_additive
#> [1] 0
#>
#> $current
#> [1] 1
#>
#> $primary_current_additive
#> [1] 1
#>
#> attr(,"class")
#> [1] "secondary_opts" "list"