繁荣背后何尝没有隐患

Posted on by

R语言的繁荣是毋庸置疑的，至少在科研界的地位蒸蒸日上，极大的占领了原来属于各种商业绘图软件的市场。不仅仅是在于其免费的属性，更重要的是随心所欲地自由定制。

但是参与的玩家多了之后，也会出现一些冲突。最近在运行一些三五年前的代码报错了，引发了我的思考。

### 正常的ID转换并不会报错

有时候我会使用 clusterProfiler 包的函数 bitr进行ID转换，代码如下：

> library(org.Mm.eg.db)
> library(clusterProfiler)
> gene <- bitr(rownames(need_DEG), fromType = "SYMBOL",
+ toType = "ENTREZID",
+ OrgDb = org.Mm.eg.db)
'select()' returned 1:1 mapping between keys and columns
Warning message:
In bitr(rownames(need_DEG), fromType = "SYMBOL", toType = "ENTREZID", :
 19.2% of input gene IDs are fail to map...

但是五年前我是不用clusterProfiler 包的，之前的代码是select函数

现在select函数就报错

如下所示：

image-20201112203109436

非常的诡异，首先它居然在没有赋值的情况下就把我的输入变量给修改了，不可思议！

实在是太不安全了！！！

略微思考了一下，猜测应该是这个select函数名字太大众了，所以在很多包里面都有，出现了冲突！

搜索看看select函数来自于哪里

可以看到优先级最高多的是dplyr包：

Help on topic 'select' was found in the following packages:

Subset columns using their names and types
(in package dplyr in library )
Ridge Regression
(in package MASS in library )
Objects exported from other packages
(in package tidygraph in library )
Objects exported from other packages
(in package clusterProfiler in library )
AnnotationDb objects and their progeny, methods etc.
(in package AnnotationDbi in library )

dplyr包设计的有问题啊，都不复制就修改了变量。鄙视它

我加上了select函数真正的来源

代码虽然复杂了一点：

library(org.Mm.eg.db)
gene_up=as.character(na.omit(AnnotationDbi::select(org.Mm.eg.db,keys = gene_up,columns = 'ENTREZID',keytype = 'SYMBOL')[,2]))
gene_down=as.character(na.omit(AnnotationDbi::select(org.Mm.eg.db,keys = gene_down,columns = 'ENTREZID',keytype = 'SYMBOL')[,2]))
head(gene_up)

但是这次就不报错了。

教程的隐患？

现在每次分析数据，都引入一大波R包，优先级关系很烦人，如果每个函数都需要显式调用，我以前的代码都会很麻烦！

​

Comments are closed.