Legacy function to compute superseded borrowing metrics — hbl

Calculate legacy/superseded borrowing metrics using summary output from a fitted borrowing model and analogous summaries from the benchmark models. hbl_ess() is preferred over hbl_metrics().

hbl_metrics(borrow, pool, independent)

Arguments

borrow: A data frame returned by hbl_summary() for the hierarchical model.
pool: A data frame returned by hbl_summary() for the pooled model.
independent: A data frame returned by hbl_summary() for the independent model.

Value

A data frame with borrowing metrics.

Examples

if (!identical(Sys.getenv("HBL_TEST", unset = ""), "")) {
set.seed(0)
data <- hbl_sim_independent(
  n_study = 2,
  n_group = 2,
  n_patient = 5,
  n_rep = 3
)$data
tmp <- utils::capture.output(
  suppressWarnings(
    mcmc_borrow <- hbl_mcmc_hierarchical(
      data,
      chains = 1,
      warmup = 10,
      iter = 20,
      seed = 0
    )
  )
)
tmp <- utils::capture.output(
  suppressWarnings(
    mcmc_pool <- hbl_mcmc_pool(
      data,
      chains = 1,
      warmup = 10,
      iter = 20,
      seed = 0
    )
  )
)
tmp <- utils::capture.output(
  suppressWarnings(
    mcmc_independent <- hbl_mcmc_independent(
      data,
      chains = 1,
      warmup = 10,
      iter = 20,
      seed = 0
    )
  )
)
borrow <- hbl_summary(mcmc_borrow, data)
pool <- hbl_summary(mcmc_pool, data)
independent <- hbl_summary(mcmc_independent, data)
hbl_metrics(
  borrow = borrow,
  pool = pool,
  independent = independent
)
}