如何将所有变量相互交叉并在 R 中收集卡方测试值？

Question

我想将数据文件中的所有 40 个分类变量相互设置（= 160 个交叉表）并收集所有卡方检验的 p 值，最好在一个列表中，以便查看哪些变量最密切相关.

是否有 R 代码以简单的方式执行此请求？

Answer 1

您可以使用comb函数来查找所有组合并相互运行任意数量的变量。

作为一个简单的解决方案，如果您有一个名为dt的data.table ，并且自变量是result ，则使用以下代码。

library(data.table)
library(magrittr)
library(dplyr)

chi_dt <- dt %>%
  map(~chisq.test(.x, dt$result)) %>%
  tibble(names = names(.), data = .) %>%
  mutate(stats = map(data, broom::tidy)) %>%
  unnest(stats)  %>% select(-data) %>%
  arrange(p.value, desc(statistic))