用首先下载GEO数据 有的下载会有几个对象 有的就只有一个 但是不影响 因为GSE就是数据的集合,这几个不同的对象只是测序平台不同而已 这里选择第二个对象
gset = getGEO('GSE12417', destdir=".",getGPL = F)
e2=gset[[2]] #取下载geo数据list中的第一个对象 可以看见它测序用哪一个GPL平台
#也可以取其他的对象,反正记住对应的测序GPL平台
expr=e2@assayData[["exprs"]] #获取表达矩阵
expr=exprs(e2) #用function直接获得了表达矩阵 二者一样
附三个对象:可以看到第二个对象测序平台为GPL96

去geo数据库搜索GSE12417,往下翻
点击GPL96进入该界面
download到当前工作路径下 ,然后读入
anno=fread("GPL96-57554.txt",header = T,data.table = F)#一定要data.table=F
gene.1=anno[,c("ID","Gene Symbol")] #获得探针对应symbol的data.frame
#由于一个探针可能对应多个sy
2704




被折叠的 条评论
为什么被折叠?



