转载:GPL平台注释,如何从gene_assinment中挖出gene symbol

GPL平台注释,如何从gene_assinment中挖出gene symbol

原文出处

https://www.bioinfo.info/?/question/6112

在做芯片分析时,得到了差异基因的list,但是读入GPL6244的注释文件发现gene symbol不是单独的一列,而是隐藏在gene_assignment这一列中间。怎么才能注释得到gene symbol结果?

回复1:这个你可以把 gene_assignment这一列取出来,然后使用 // 分割, 取分割以后的第1列即可。R里面limma包有个函数好像叫  strsplit2() ,你可以试试。

回复2:

library(tidyverse)

symbol<-str_split_fixed(annotation_file$gene_assignment,pattern = "//",3)[,2]

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

友情链接更多精彩内容