我对一个组变量的多个t检验有疑问。
我有一个基因表达数据框，仅举一个例子：
df <- data.frame(a=runif(100), b=runif(100), c=runif(100)+0.5, d=runif(100)+0.5)
d <- melt(df)
然后我有一个带有分组变量的数据框。像这样：
dg <- data.frame(c(rep("A",50), rep("B",50)))
names(dg)<-"Group"
我发现了几个非常相似的问题，并且实现了对第一个数据帧的每一列和分组变量使用 lapply 进行多次t检验。
但是，我认为有必要调整p值以进行多次比较（我总共有15列）。这是我的问题，我该怎么办？例如与Bonferroni？我可以在本《守则》中实现吗？
此外，我认为也可以使用多个数据列（正态分布/等方差）对假设进行t检验。
这里有人知道如何以一种奇妙的方式将所有这些结合起来吗？完美的脚本应该是这样的脚本，其中统计信息也可以在数据帧或类似的格式中转换。
谢谢您的帮助！