Define the properties of parameter distributions and their correlation structure for probabilistic uncertainty analysis of Markov models.
Usage
define_psa(..., correlation)
define_psa_(.dots = list(), correlation)
Arguments
- ...
Formulas defining parameter distributions.
- correlation
A correlation matrix for parameters or the output of
define_correlation()
.- .dots
Pair/values of expressions coercible to quosures.
Value
An object of class resamp_definition
.
Contains list_qdist
, a list of quantile
functions and correlation
a correlation matrix.
Details
The distributions must be defined within heemod
(see distributions), or defined with
define_distribution()
.
If no correlation matrix is specified parameters are assumed to be independant.
The correlation matrix need only be specified for correlated parameters.
Examples
mc <- define_correlation(
age_init, cost_init, .4
)
define_psa(
age_init ~ normal(60, 10),
cost_init ~ normal(1000, 100),
correlation = mc
)
#> A PSA definition:
#>
#> 2 parameters resampled, 0 multinomial group.
# example with multinomial parameters
define_psa(
rate1 + rate2 + rate3 ~ multinomial(10, 50, 40),
a + b ~ multinomial(15, 30)
)
#> A PSA definition:
#>
#> 5 parameters resampled, 2 multinomial groups.