image.png
这道题很简单，但是应该更重要的是学会用R的方法，我后面会用R把这个系列的题目重做一遍，这里先用python做。
ID转换主要就是构建对应关系的字典
先看一下文件的内容
对应关系文件
head geneid2symbol.txt1 A1BG2 A2M3 A2MP19 NAT110 NAT211 NATP12 SERPINA313 AADAC14 AAMP15 AANAT
geneID文件
head my_geneID.txt1175 2263 1581 4091 1783 2051 2248 22941 1410 881
脚本
import collections
geneDict=collections.OrderedDict()with open ("geneid2symbol.txt") as fh: for line in fh:
lineL=line.strip().split("\t")
gene_id=lineL[0]
symbol=lineL[1]
geneDict[gene_id]=symbolfor line in open("my_geneID.txt"):
line=line.strip()
print(geneDict[line])运行结果
python3 exchange.py AP2S1 FGFR2 CYP7A1 SMAD6 DYNC1LI2 EPHB6 FGF3 SHANK2 CRYAB CCIN CHKB TSPAN2 YAF2 ZFC3H1 ATP6V0E2 LOC101929524 PCCA HLA-DRB6 TMC6 JTB POLD1 LHCGR WBP4 FGL2 DEPDC5 EAPP RRAS B3GALT5 DESI2SCN7A
其他的，geneID,symbol,ensembleID,entrez ID,probe ID 等等之间的转换，方法是一样的
作者：天秤座的机器狗
链接：https://www.jianshu.com/p/527fcdbb2d11
點(diǎn)擊查看更多內(nèi)容
為 TA 點(diǎn)贊
評(píng)論
評(píng)論
共同學(xué)習(xí)，寫下你的評(píng)論
評(píng)論加載中...
作者其他優(yōu)質(zhì)文章
正在加載中
感謝您的支持，我會(huì)繼續(xù)努力的～
掃碼打賞，你說多少就多少
贊賞金額會(huì)直接到老師賬戶
支付方式
打開微信掃一掃，即可進(jìn)行掃碼打賞哦