Function to run two-sample Wilcoxon tests on multiple variables across multiple grouping variables.

Source: R/grouped_wilcox.R

Function to run two-sample Wilcoxon tests on multiple variables across multiple grouping variables.

Running Wilcox test across multiple grouping variables.

grouped_wilcox(
  data,
  dep.vars,
  indep.vars,
  grouping.vars,
  paired = FALSE,
  correct = TRUE
)

Arguments

data

Dataframe from which variables are to be taken.
dep.vars

List dependent variables for a two-sample Wilcoxon tests (y in y ~ x).
indep.vars

List independent variables for a two-sample Wilcoxon tests (x in y ~ x).
grouping.vars

List of grouping variables (if NULL, the entire dataframe will be used).
paired

A logical indicating whether you want a paired two-sample Wilcoxon tests (Default: paired = FALSE).
correct

A logical indicating whether to apply continuity correction in the normal approximation for the p-value (Default: correct = TRUE).

Value

A tibble dataframe with tidy results from two-sample Wilcoxon tests analyses.

See also

grouped_tidy

Examples

# for reproducibility
set.seed(123)

# only with one grouping variable
groupedstats::grouped_wilcox(
  data = dplyr::filter(.data = ggplot2::diamonds, color == "E" | color == "J"),
  dep.vars = depth:table,
  indep.vars = color,
  grouping.vars = clarity,
  paired = FALSE
)

#> # A tibble: 16 x 9
#>    clarity formula  method        statistic estimate conf.low conf.high  p.value
#>    <ord>   <chr>    <chr>             <dbl>    <dbl>    <dbl>     <dbl>    <dbl>
#>  1 SI2     depth ~… Wilcoxon ran…   380600  -2.00e-1 -3.00e-1  -3.64e-5 1.54e- 2
#>  2 SI1     depth ~… Wilcoxon ran…   918892   3.93e-5 -1.00e-1   1.00e-1 6.77e- 1
#>  3 VS1     depth ~… Wilcoxon ran…   276743  -4.00e-1 -5.00e-1  -3.00e-1 7.03e-12
#>  4 VVS2    depth ~… Wilcoxon ran…    52618  -4.00e-1 -6.00e-1  -2.00e-1 4.18e- 4
#>  5 VS2     depth ~… Wilcoxon ran…   816784. -2.00e-1 -3.00e-1  -1.00e-1 8.86e- 5
#>  6 I1      depth ~… Wilcoxon ran…     1570. -2.00e+0 -3.00e+0  -1.00e+0 1.21e- 4
#>  7 VVS1    depth ~… Wilcoxon ran…    19455  -3.00e-1 -5.00e-1  -1.00e-1 5.07e- 3
#>  8 IF      depth ~… Wilcoxon ran…     3286  -3.00e-1 -7.00e-1  -2.51e-5 4.79e- 2
#>  9 SI2     table ~… Wilcoxon ran…   388138. -2.08e-5 -3.91e-5   4.98e-6 6.85e- 2
#> 10 SI1     table ~… Wilcoxon ran…   880016.  0       -1.30e-5   3.62e-5 1.71e- 1
#> 11 VS1     table ~… Wilcoxon ran…   306400  -4.82e-6 -1.00e+0  -2.87e-5 6.16e- 5
#> 12 VVS2    table ~… Wilcoxon ran…    54524  -1.00e+0 -1.00e+0  -5.87e-5 2.60e- 3
#> 13 VS2     table ~… Wilcoxon ran…   763942. -1.00e+0 -1.00e+0  -4.61e-5 1.65e-10
#> 14 I1      table ~… Wilcoxon ran…     2900   1.00e+0 -2.80e-5   2.00e+0 1.68e- 1
#> 15 VVS1    table ~… Wilcoxon ran…    22659  -7.72e-6 -1.00e+0   4.05e-6 3.42e- 1
#> 16 IF      table ~… Wilcoxon ran…     3987  -3.35e-5 -1.00e+0   1.00e+0 9.11e- 1
#> # … with 1 more variable: significance <chr>