log_prior: Evaluate log-prior 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

Evaluate log-prior of calibrated parameters

Description

log_prior computes a log-prior value for one (or multiple) parameter set(s) based on their prior distributions.

Usage

log_prior(
  v_params,
  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

v_params	Vector (or matrix) of model parameters.
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 scalar (or vector) with log-prior values.

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
log_prior(v_params = sample.prior(n_samp = 5))

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