我有一个包含 2304 个值的大数据集:384 个不同的基因表达值,一式三份,针对 2 个不同的组(即不同的治疗)。我只想比较两组之间相同的基因:例如,基因 A 处理与基因 A 对照,等等......并获得 p 值。我不是想比较不同的基因!
所以我正在运行这个命令:
data2 %>%
group_by(group) %>%
wilcox_test(expression ~ gene) %>%
adjust_pvalue(method = "BH") %>%
add_significance()
但是,在 R 上运行需要一生的时间,所以我正在为这个特定问题寻找更好的函数?
谢谢
我使用了适合我的数据分布并考虑多个因素的 GLM,而不是使用 Wilcoxon 检验。