我使用R中的pairwise.t.test函数进行了一系列T检验,为所有比较创建了一个P值的时髦矩阵,如
date p.value.2016.04.18 p.value.2016.04.20 p.value.2016.04.22
2016-04-20 1.00E+00 NA NA
2016-04-22 1.00E+00 1.22E-01 NA
2016-04-24 1.00E+00 1.00E+00 1.00E+00
2016-04-26 1.00E+00 5.01E-01 1.00E+00
2016-04-28 1.28E-02 1.03E-01 6.20E-06
我想在组的平均值之间建立匹配的差异矩阵。我的google-fu让我失望了。你可以帮助一个简陋的统计数据吗?
答案 0 :(得分:1)
以下是此famous listserv的示例:
Table1 <- matrix(10:6, ncol = 1)
rownames(Table1) <- letters[1:5]
Table1
t(outer(Table1[,1], Table1[,1], `-`))