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
それを行う簡単な方法はありますか?
乾杯!