基因探针的ID转换

通过r语言将基因探针的id转换转换成对应的基因

rm(list=ls())
options(stringsAsFactors = F)
a=read.table('tz.txt')
colnames(a)="probe_id"
ids=toTable(hgu133aSYMBOL)##每个平台所对应的包都不一样,需要从网上下载所对应的包
head(ids)
b=merge(a,ids,by='probe_id')
c=ids[match(a$probe_id,ids$probe_id),]
b==c
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容