轉(zhuǎn)載:GPL平臺注釋,如何從gene_assinment中挖出gene symbol

GPL平臺注釋,如何從gene_assinment中挖出gene symbol

原文出處

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

在做芯片分析時,得到了差異基因的list,但是讀入GPL6244的注釋文件發(fā)現(xiàn)gene symbol不是單獨的一列,而是隱藏在gene_assignment這一列中間。怎么才能注釋得到gene symbol結(jié)果?

回復1:這個你可以把 gene_assignment這一列取出來,然后使用 // 分割, 取分割以后的第1列即可。R里面limma包有個函數(shù)好像叫? strsplit2() ,你可以試試。

回復2:

library(tidyverse)

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

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時請結(jié)合常識與多方信息審慎甄別。
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務。

友情鏈接更多精彩內(nèi)容