A pipe-friendly function to add an adjusted p-value column into a data frame. Supports grouped data.

adjust_pvalue(data, p.col = NULL, output.col = NULL, method = "holm")

Arguments

data

a data frame containing a p-value column

p.col

column name containing p-values

output.col

the output column name to hold the adjusted p-values

method

method for adjusting p values (see p.adjust). Allowed values include "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". If you don't want to adjust the p value (not recommended), use p.adjust.method = "none".

Value

a data frame

Examples

# Perform pairwise comparisons and adjust p-values ToothGrowth %>% t_test(len ~ dose) %>% adjust_pvalue()
#> # A tibble: 3 x 10 #> .y. group1 group2 n1 n2 statistic df p p.adj p.adj.signif #> <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <dbl> <chr> #> 1 len 0.5 1 20 20 -6.48 38.0 1.27e- 7 2.54e- 7 **** #> 2 len 0.5 2 20 20 -11.8 36.9 4.40e-14 1.32e-13 **** #> 3 len 1 2 20 20 -4.90 37.1 1.91e- 5 1.91e- 5 ****