Phase 1.g delegates to rbounds when installed and the
wilcoxon or sign method is requested; otherwise
falls back to the base-R normal-approximation implementation
originally shipped with rmorie. The mcnemar path is always
served by the inline binomial formula (rbounds does not expose a
McNemar entry point on CRAN).
Usage
rosenbaum_bounds(
treated_outcomes,
control_outcomes,
gamma_range = NULL,
method = "wilcoxon"
)