chi-square-tests-for-variable-pairs.R

library(dplyr)
library(broom)
library(tidyr)

# Make example data
  my_data <- data.frame(x = sample(1:5, size = 10, replace = TRUE),
                        y = sample(1:5, size = 10, replace = TRUE))

# Generate every pair of variables
  var_list <- colnames(my_data)
  var_combinations <- expand.grid(Var_1 = var_list, Var_2 = var_list,
                                  stringsAsFactors = FALSE)
  
  var_combinations |>
    # Group data by rows
    rowwise() |>
    # Do a Chi-square test and tidy the output
    mutate(
      chisq_test = chisq.test(x = my_data[[Var_1]],
                              y = my_data[[Var_2]]) |>
        broom::tidy()
    ) |>
    unnest(cols = chisq_test)