跨多列非配对t检验的R代码

跨多列非配对t检验的R代码,r,statistics,t-test,R,Statistics,T Test,我希望在以下代码中调整p值(fdr方法): >lstOut <- combn(unique(df1$ATTRIBUTE_Group), 2, FUN = function(x) {+ dat1 <- subset(df1, ATTRIBUTE_Group == x[1])+ dat2 <- subset(df1, ATTRIBUTE_Group == x[2])+ Map(function(x, y) tryCatch(t.test(x, y,var.equal=TRUE

我希望在以下代码中调整p值(fdr方法):

>lstOut <- combn(unique(df1$ATTRIBUTE_Group), 2, FUN = function(x) {+ dat1 <- subset(df1, ATTRIBUTE_Group == x[1])+ dat2 <- subset(df1, ATTRIBUTE_Group == x[2])+ Map(function(x, y) tryCatch(t.test(x, y,var.equal=TRUE), error = function(e) NA), dat1[-1], dat2[-1])+ }, simplify = FALSE)
>names(lstOut) <- combn(as.character(unique(df1$ATTRIBUTE_Group)), 2, FUN = paste, collapse="_")
>out <- map_depth(lstOut, .depth = 2, tidy)
>out2 <- map_dfr(out, ~bind_rows(.x, .id = 'colname'), .id = 'classCompare')

>lstOut我们可以使用
p.adjust

library(dplyr)
out3 <- out2 %>%
              group_by(colname) %>%
              mutate(padj = p.adjust(p.value, method = "BH"))
库(dplyr)
出局3%
分组依据(colname)%>%
突变(padj=p.adjust(p.value,method=“BH”))

可能是
out2%>%groupby(colname)%>%mutate(padj=p.adjust(p.value,method=“BH”))
我不确定这是否有效。我的p.values changedIt没有创建新列“padj”,即
out3%groupby(colname)%>%mutate(padj=p.adjust(p.value,method=“BH”);out3$padj