wilcoxon test in loop with rstatix

I have data frame like

library(rstatix)
grp1 <- runif(10, min = 0, max = 100)  # Exemple de données numériques
grp2 <- runif(10, min = 0, max = 100)
grp3 <- runif(10, min = 0, max = 100)
grp4 <- runif(10, min = 0, max = 100)
status <- sample(c("alive", "death", "sick"), 10, replace = TRUE)  # Exemple de données caractère

my_data_frame <- data.frame(grp1, grp2, grp3, grp4, status)

and I want to do for each group.

  wilcox_test(grp1~ status, p.adjust.method = "fdr")

I try

for (i in 4) {
 
  var <- names(my_data_frame)
  test <- my_data_frame %>% wilcox_test(var[i] ~ status, p.adjust.method = "fdr")}

but Column var[i] doesn't exist.

with other test like glm this syntax works

Solution

rstatix is designed to fit nicely into Tidyverse-based workflows and pipelines, so instead of a loop, you could first pivot your dataset to longer, group and pass grouped tibble to wilcox_test():

library(rstatix)
library(dplyr)
library(tidyr)

grp1 <- runif(10, min = 0, max = 100)  # Exemple de données numériques
grp2 <- runif(10, min = 0, max = 100)
grp3 <- runif(10, min = 0, max = 100)
grp4 <- runif(10, min = 0, max = 100)
status <- sample(c("alive", "death", "sick"), 10, replace = TRUE)  # Exemple de données caractère
my_data_frame <- data.frame(grp1, grp2, grp3, grp4, status)

my_data_frame |>
  pivot_longer(-status, names_to = "grp") |>
  group_by(grp) |>
  wilcox_test(value ~ status, p.adjust.method = "fdr")
#> # A tibble: 12 × 10
#>    grp   .y.   group1 group2    n1    n2 statistic     p p.adj p.adj.signif
#>  * <chr> <chr> <chr>  <chr>  <int> <int>     <dbl> <dbl> <dbl> <chr>       
#>  1 grp1  value alive  death      3     3         3 0.7   0.857 ns          
#>  2 grp1  value alive  sick       3     4         3 0.4   0.857 ns          
#>  3 grp1  value death  sick       3     4         5 0.857 0.857 ns          
#>  4 grp2  value alive  death      3     3         3 0.7   0.7   ns          
#>  5 grp2  value alive  sick       3     4         8 0.629 0.7   ns          
#>  6 grp2  value death  sick       3     4        12 0.057 0.171 ns          
#>  7 grp3  value alive  death      3     3         7 0.4   0.629 ns          
#>  8 grp3  value alive  sick       3     4         8 0.629 0.629 ns          
#>  9 grp3  value death  sick       3     4         4 0.629 0.629 ns          
#> 10 grp4  value alive  death      3     3         3 0.7   0.7   ns          
#> 11 grp4  value alive  sick       3     4         4 0.629 0.7   ns          
#> 12 grp4  value death  sick       3     4         8 0.629 0.7   ns

^{Created on 2024-02-09 with reprex v2.1.0}