以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]]
![](https://cdn.static.attains.cn/app/developer-bbs/upload/1723128412456105877.jpeg)
3 用GEOquery里的pdata函数获取样本信息
看一下pdata的结构,很明显是数据框
代码语言:javascript
复制
pdata<-pData(Gset)
class(pdata)
View(pdata)
![](https://cdn.static.attains.cn/app/developer-bbs/upload/1723128412545529835.jpeg)
4 dim查看行列 colnames查看列名
代码语言:javascript
复制
dim(pdata)
colnames(pdata)
方法二
1.百度搜索GEO
![](https://cdn.static.attains.cn/app/developer-bbs/upload/1723128413144943677.jpeg)
2.进入sample
![](https://cdn.static.attains.cn/app/developer-bbs/upload/1723128413538495211.jpeg)
3.输入GEO号search一下,导出即可
![](https://cdn.static.attains.cn/app/developer-bbs/upload/1723128413848566864.jpeg)