R语言GEO数据处理(四)
# 3. id转换 -----------------------------------------------------------------
##方法一:使用R包转换
index = gset[[1]]@annotation
if(!require("hgu133a.db"))
BiocManager::install("hgu133a.db")
library(hgu133a.db)
ls("package:hgu133a.db")
ids <- toTable(hgu133aSYMBOL)
head(ids)
length(unique(ids$symbol))
table(sort(table(ids$symbol)))
###id转换
library(tidyverse)
exp <- as.data.frame(exp)
exp <- exp %>%
mutate(probe_id=rownames(exp)) %>%
inner_join(ids,by="probe_id") %>%
select(probe_id, symbol, everything())
exp2 <- exp[!duplicated(exp$symbol),]
rownames(exp2) <- exp2$symbol
exp2 <- exp[,-(1:2)]
exp2[1:3,1:3]
##方法二:使用soft文件注释
##方法三:手工注释
赞 (0)