我已经对 92 种蛋白质进行了方差分析，以同时使用
x <- data.frame(p.value=sapply(df[,2:93], function(i) summary(aov(i ~ df$X.Alcohol))[[1]][["Pr(>F)"]][[1]]))
pAlcohol<-x[x$p.value < 0.05/92,,drop=F]
获得组间显着不同的蛋白质。我现在想做一个事后测试来了解这些组的不同之处。这有效：
y <- data.frame(post=sapply(df[,2:93], function(i) TukeyHSD((aov(i ~ df$X.Alcohol)), ordered = TRUE)))
但它适用于所有蛋白质，而不仅仅是那些重要的蛋白质。我希望它只测试 pAlcohol 中的蛋白质。
现在，当它对所有 92 种蛋白质进行事后测试时，它会将它们列在三行上（我有 3 组）。是否可以为每种蛋白质添加 3 行新行？