以GSE111229为例
方法一
1 下载并保存GEO数据
下载有error可以:设置镜像、访问外国网站、rm(list=ls())一下重试…
代码语言:javascript
复制
library(GEOquery)
GSE_name = 'GSE111229'
options( 'download.file.method.GEOquery' = 'libcurl' ) #windows系统
gset <- getGEO( GSE_name, getGPL = F )
save( gset, file = 'gset.Rdata' )
2 加载GEO数据
gset包含下载的所有信息
由于gset是列表,故将其转为可操作的数据结构Gset
代码语言:javascript
复制
load("gset.Rdata")
Gset <- gset[[1]]
3 用GEOquery里的pdata函数获取样本信息
看一下pdata的结构,很明显是数据框
代码语言:javascript
复制
pdata<-pData(Gset)
class(pdata)
View(pdata)
4 dim查看行列 colnames查看列名
代码语言:javascript
复制
dim(pdata)
colnames(pdata)
方法二
1.百度搜索GEO
2.进入sample
3.输入GEO号search一下,导出即可