applyBenjaminiHochberg

applyBenjaminiHochberg

applyBenjaminiHochberg(df_fisher, Q = 0.05)

Arguments

df_fisher

The prior Fisher result data with gene & p_value cols

Q

The FDR threshold (default = 0.05)

Value

Returns Fisher results with added cols for adj_p_value, sig_after_bh, Q

Examples

fisher_results <- tibble::tibble(
  gene    = paste0("gene_", 1:5),
  p_value = c(0.001, 0.01, 0.04, 0.2, 0.5)
)
applyBenjaminiHochberg(fisher_results, Q = 0.05)
#> # A tibble: 5 × 5
#>   gene   p_value adj_p_value sig_after_bh     Q
#>   <chr>    <dbl>       <dbl> <lgl>        <dbl>
#> 1 gene_1   0.001      0.005  TRUE          0.05
#> 2 gene_2   0.01       0.025  TRUE          0.05
#> 3 gene_3   0.04       0.0667 FALSE         0.05
#> 4 gene_4   0.2        0.25   FALSE         0.05
#> 5 gene_5   0.5        0.5    FALSE         0.05