4

GEOquery is a great R package to retrieve and analyze the Gene Expression data stored in NCBI Gene Expression Omnibus (GEO) database. I have used the following code provided from GEO2R service of GEO database (that generates the initial R script to analyze your desired data automatically) to extract some GEO series of experiments:

gset <- getGEO("GSE10246", GSEMatrix =TRUE)
if (length(gset) > 1) idx <- grep("GPL1261", attr(gset, "names")) else idx <- 1
    gset <- gset[[idx]]
gset # displays a summary of the data stored in this variable

The problem is that I can not retrieve the sample titles from it. I have found some function Columns() that works on GDS datasets and returns the sample names, but not on GSE.

Please note I am not interested in sample accession IDs (i.e. GSM258609 GSM258610, etc), what I want is the sample human readable titles.

Is there any idea? Thanks

4

1 回答 1

3

gset <- getGEO("GSE10246", GSEMatrix =TRUE)

gset是一个简单的列表,它的第一个元素是 an ExpressionSet,示例信息在phenoDataorpData中,所以也许您正在寻找

pData(gset[[1]])

查看?ExpressionSet更多。

于 2013-01-10T23:05:34.377 回答