将pairwise.wilcox.test测试的输出转换为列表

Question

我想转换来自pairwise.wilcox.test的输出矩阵：

> df <- data.frame(obs = c(rep("a", 10), rep("b", 10), rep("c", 10), rep("d", 10)), 

                   vars = c(runif(10, 1, 3), runif(10, 1, 6), runif(10, 3, 7), runif(10, 2, 9)))

> pairwise.wilcox.test(df$vars, df$obs)$p.value
             a          b         c
b 6.299985e-01         NA        NA
c 6.495053e-05 0.02052437        NA
d 6.495053e-05 0.01554483 0.6299985

添加到一个列表中，例如：

a  b  6.299985e-01
a  c  6.495053e-05
a  d  6.495053e-05
b  c  0.02052437
b  d  0.01554483
c  d  0.6299985

有什么简单的方法做到这一点吗？

干杯!

Answer 1

你可以使用这个

temp <-pairwise.wilcox.test(df$vars, df$obs)$p.value
df <- data.frame(expand.grid(dimnames(temp)),array(temp))
na.omit(df)