pval_correct {sparrpowR} | R Documentation |
Calculate p-value corrections for multiple testing
Description
Internal function to calculate various p-value corrections for use within the spatial_power
and jitter_power
functions.
Usage
pval_correct(
input,
type = c("FDR", "Sidak", "Bonferroni"),
alpha = 0.05,
nbc = NULL
)
Arguments
input |
An object of class 'rrs' from the |
type |
Character string specifying which correction for multiple comparisons. Options include a False Discovery Rate |
alpha |
Numeric. The alpha level for significance threshold (default in |
Details
This function provides functionality for multiple testing correction in five ways:
Computes a False Discovery Rate by Benjamini and Hochberg doi:10.1111/j.2517-6161.1995.tb02031.x (
p_correct = "FDR"
) by: 1) sorting the p-values (p_i) of each knot in ascending order (p_1 <= p_2 <= ... <= p_m), 2) starting from p_m find the first p_i for which p_i <= (i/m) * alpha.Computes a Sidak correction doi:10.2307/2283989 (
p_correct = "Sidak"
) by 1 - (1 -alpha
) ^ (1 / total number of gridded knots across the estimated surface). The default in therisk
function is a resolution of 128 x 128 or n = 16,384 knots and a custom resolution can be specified using theresolution
argument within therisk
function.Computes a Bonferroni correction (
p_correct = "Bonferroni"
) byalpha
/ total number of gridded knots across the estimated surface. The default in therisk
function is a resolution of 128 x 128 or n = 16,384 knots and a custom resolution can be specified using theresolution
argument within therisk
function.
Value
An object of class 'numeric' with the corrected alpha level.