Robust to arbitrary correlation structure. No clean CRAN
equivalent (the ACAT package is GitHub-only); kept as an in-house
implementation.
Usage
cauchy_combination(p_values, weights = NULL)
Arguments
- p_values
Numeric vector of raw p-values.
- weights
Optional non-negative weights summing to 1.
Value
An object of class "morie_multiple_testing_result".
Examples
res <- cauchy_combination(c(0.01, 0.3, 0.5))
res$p_value
#> [1] 0.0292573
resw <- cauchy_combination(c(0.01, 0.3, 0.5), weights = c(0.5, 0.25, 0.25))
resw$p_value
#> [1] 0.01975536
