trial_totals: Trial Probability and Expected Counts

In mcmodule: Modular Monte Carlo Risk Analysis

Trial Probability and Expected Counts

Description

Calculates probabilities and expected counts across hierarchical levels (trial, subset, set) in a structured population. Uses trial probabilities and handles nested sampling with conditional probabilities.

Usage

trial_totals(
  mcmodule,
  mc_names,
  trials_n,
  subsets_n = NULL,
  subsets_p = NULL,
  name = NULL,
  prefix = NULL,
  combine_prob = TRUE,
  all_suffix = NULL,
  level_suffix = c(trial = "trial", subset = "subset", set = "set"),
  mctable = set_mctable(),
  agg_keys = NULL,
  agg_suffix = NULL,
  keep_variates = FALSE,
  summary = TRUE,
  data_name = NULL
)

Arguments

mcmodule	mcmodule object containing input data and node structure
mc_names	Vector of node names to process
trials_n	Trial count column name
subsets_n	Subset count column name (optional)
subsets_p	Subset prevalence column name (optional)
name	Custom name for output nodes (optional)
prefix	Prefix for output node names (optional)
combine_prob	Combine probability of all nodes assuming independence (default: TRUE)
all_suffix	Suffix for combined node name (default: "all")
level_suffix	A list of suffixes for each hierarchical level (default: c(trial="trial",subset="subset",set="set"))
mctable	Data frame containing Monte Carlo nodes definitions (default: set_mctable())
agg_keys	Column names for aggregation (optional)
agg_suffix	Suffix for aggregated node names (default: "hag")
keep_variates	whether to preserve individual values (default: FALSE)
summary	Include summary statistics if TRUE (default: TRUE)
data_name	Data name used to create trials_n, subsets_n and subsets_p nodes if they don't exist in mcmodule (optional)

Value

Updated mcmodule object containing:

• Combined node probabilities

• Probabilities and counts at trial level

• Probabilities and counts at subset level

• Probabilities and counts at set level

Examples

imports_mcmodule <- trial_totals(
  mcmodule = imports_mcmodule,
  mc_names = "no_detect_a",
  trials_n = "animals_n",
  subsets_n = "farms_n",
  subsets_p = "h_prev",
  mctable = imports_mctable
)
print(imports_mcmodule$node_list$no_detect_a_set$summary)

mcmodule documentation built on Nov. 25, 2025, 5:10 p.m.