我可以应用以下代码在某个数据帧 A 的每一列上针对依赖参数 DEP 应用 pairwiset.t.test 函数。
apply(A, 2, function(x) pairwise.t.test(DEP, x, p.adj = 'bonferroni')$p.value)
此函数的输出产生一个列表类型的数据。如下所示:
$A
1 2 3 4
1 0.7236751 NA NA NA
2 0.7236751 0.7236751 NA NA
3 0.7236751 0.7236751 0.7236751 NA
4 0.7236751 0.7236751 0.7236751 0.7236751
$B
0
1 0.7236751
我想问一下是否有一种方法可以组合这两个列表模式数据,例如:
1 2 3 4 5
1 0.7236751 NA NA NA 0.7236751
2 0.7236751 0.7236751 NA NA NA
3 0.7236751 0.7236751 0.7236751 NA NA
4 0.7236751 0.7236751 0.7236751 0.7236751 NA
的输出dput(apply(A, 2, function(x) pairwise.t.test(DEP, x, p.adj = 'bonferroni')$p.value))
如下:
Dimnames = list(c("10", "11", "14", "15", "16", "17",
"18", "19", "20", "21", "23", "31", "33", "34", "36", "39", "40",
"42", "43", "45", "46", "47", "49", "50", "51", "52", "54", "59",
"60", "63", "64", "66", "70", "73"), c("6", "10", "11", "14",
"15", "16", "17", "18", "19", "20", "21", "23", "31", "33", "34",
"36", "39", "40", "42", "43", "45", "46", "47", "49", "50", "51",
"52", "54", "59", "60", "63", "64", "66", "70"))), Gender = structure(0.723675133647025, .Dim = c(1L,
1L), .Dimnames = list("1", "0"))), .Names = c("A", "B"
))