R语言DESeq获取显著差异基因名报错，求替代方法

嘿，这个错误我之前也踩过坑！咱们先搞清楚问题出在哪：

你收到的Error in results[which(resA0vsG0$padj < 0.01), ] : object of type 'closure' is not subsettable错误，是因为你直接用了results这个DESeq2的函数名去做子集化操作，但results是函数，不是你存储差异分析结果的变量——你实际的结果存在resA0vsG0里对吧？

先确认你已经把差异分析结果正确赋值给变量了（要是没这一步，先补上）：

# 把差异分析结果存入resA0vsG0变量
resA0vsG0 <- results(dds1and2, contrast = c("Time","A0","G0"))

接下来给你几种靠谱的获取显著差异基因名称的方法：

方法一：直接子集化提取行名 #

这是最直接的方式，对resA0vsG0做筛选后提取行名（行名就是基因名）：

# 筛选padj < 0.01的基因，提取名称
sig_gene_names <- rownames(resA0vsG0[which(resA0vsG0$padj < 0.01), ])

方法二：用subset函数更直观 #

如果你觉得which的写法有点绕，用subset函数可读性更强：

# 先筛选出显著差异的结果子集
sig_res <- subset(resA0vsG0, padj < 0.01)
# 再提取基因名称
sig_gene_names <- rownames(sig_res)

方法三：排除NA值干扰 #

有时候部分基因的padj会是NA（比如表达量太低没统计意义），可以加上排除条件避免这些基因混入：

# 同时排除padj为NA的情况
sig_gene_names <- rownames(resA0vsG0[which(!is.na(resA0vsG0$padj) & resA0vsG0$padj < 0.01), ])

这样就能顺利拿到你想要的显著差异基因名称啦！

内容的提问来源于stack exchange，提问作者cmk