sample.prior: Sample from prior distributions of calibrated parameters

In feralaes/cdx2cea: A cost-efectiveness analysis of testing stage II colon cancer patients for the absence of CDX2 biomarker followed by adjuvant chemotherapy

Sample from prior distributions of calibrated parameters

Description

sample.prior generates a sample of parameter sets from their prior distribution.

Usage

sample.prior(
  n_samp,
  v_param_names = c("r_DieMets", "r_RecurCDX2pos", "hr_RecurCDX2neg", "p_Mets"),
  v_lb = c(r_DieMets = 0.037, r_RecurCDX2pos = 0.001, hr_RecurCDX2neg = 1.58, p_Mets =
    0.9),
  v_ub = c(r_DieMets = -log(1 - (1 - 0.03))/60, r_RecurCDX2pos = 0.03, hr_RecurCDX2neg
    = 4.72, p_Mets = 0.99)
)

Arguments

n_samp	Number of samples.
v_param_names	Vector with parameter names.
v_lb	Vector with upper bounds for each parameter.
v_ub	Vector with lower bounds for each parameter.

Value

A matrix with 3 rows and n_samp rows. Each row corresponds to a parameter set sampled from their prior distributions

Examples

v_param_names  <- c("r_DieMets",
                    "r_RecurCDX2pos",
                    "hr_RecurCDX2neg",
                    "p_Mets")
n_param        <- length(v_param_names)
v_lb <- c(r_DieMets       = 0.037, 
          r_RecurCDX2pos  = 0.001,
          hr_RecurCDX2neg = 1.58, 
          p_Mets          = 0.9))  # lower bound
v_ub <- c(r_DieMets       = -log(1-(1-0.03))/60, 
          r_RecurCDX2pos  = 0.03,
          hr_RecurCDX2neg = 4.72, 
          p_Mets          = 0.99) # upper bound
sample.prior(2)

feralaes/cdx2cea documentation built on April 7, 2024, 10:12 a.m.