linux收入即学习

ID转换

2020-05-12  本文已影响0人  找兔子的小萝卜

circRNA

十行代码搞定

1 下载GPL文件,代码本地下载均可

library(data.table)
a=fread('probeMatrix.txt',data.table = F)
a[1:4,1:4]
b=read.table('ann.txt',sep = '\t',header = T)
tail(head(b,20))
d=merge(a,b,by.x='ID_REF',by.y='ID')
e=read.table('ID.txt',header = T)
head(e)
f=merge(e,d,by='circRNA')
head(f[,1:6])

关键是文件的读取

mRNA

三种方法

方法1 BioconductorR包

gpl 
#http://www.bio-info-trainee.com/1399.html
if(!require(hugene10sttranscriptcluster.db))BiocManager::install("hugene10sttranscriptcluster.db")
library(hugene10sttranscriptcluster.db)
ls("package:hugene10sttranscriptcluster.db")
ids <- toTable(hugene10sttranscriptclusterSYMBOL)
head(ids)

方法2 读取gpl页面的soft文件,按列取子集

方法3 官网下载

方法4 自主注释

上一篇 下一篇

猜你喜欢

热点阅读