我正在尝试运行一个基因集富集分析,如下图所示:http://yulab-smu.top/clusterProfiler-book/chapter7.html
我一直在尝试使用这个特定的代码:
data(df, package="DOSE")
gene <- names(df)[abs(df) > 2]
head(gene)当我运行上面的代码时,我得到了以下错误:
data set ‘df’ not found[1] NA NA NA NA NA NA
[1] FALSE我正在尝试对Log2FoldChange (从-x到+x的值)的数据集应用截止值1。然后,我想将它和相应的名称(ENTREZ ID)放入一个名为gene的数据框中,这样我就可以运行富集分析了。
如果我指定数据的列:
#Define DEG as LFC>1
data(df, package="DOSE")
gene <- names(df$ENTREZ)[abs(df$log2FoldChange) > 2]
head(gene)
gene <- (abs(1) > 1)
head(gene)我得到了:
data set ‘df’ not foundNULL
[1] FALSE上面的'abs()‘函数不起作用(我已经复制并修改了上面的vignette)。
我在这里找到了另一种解释,它可能暗示了一条不同的路线。Ordering multiple columns using cut-off values in R
我是不是错过了一些语法来完成这项工作?任何帮助都将不胜感激。
谢谢
发布于 2021-03-17 21:00:59
出现错误是很正常的,因为data(df, package="DOSE")是无效的命令。您应该将df替换为您的数据帧。
https://stackoverflow.com/questions/66673437
复制相似问题