How do I creat an expressionset if there is no GDS listed on the GEO database?
1
0
Entering edit mode
pgurkirat • 0
@pgurkirat-20566
Last seen 3.1 years ago

Hello I have been using the GEOquery package to create expression sets from data taken from NCBI gene expression omnibus (GEO) database. However, I have ran into a problem. Most of the studies that are useful to me do not have a GDS listed. They have GSE, GSM, and GPL but, I use the command gds<-getGEO("GDS4061") to start with, and then do eset <- GDS2eSet(gds,do.log2=TRUE) to convert GDS into an expression set, so I assume the GDS class is necessary. Can somebody suggest how to overcome this problem please?

limma microarray GEOquery • 573 views
0
Entering edit mode

Cross-posted on Biostars: https://www.biostars.org/p/375683/

1
Entering edit mode
@james-w-macdonald-5106
Last seen 4 minutes ago
United States

Have you read the vignette? You can just as easily get data from the GSE entry, or if you want the raw data you can use getGEOSuppFiles. I have personally never used the pipeline you outline, opting instead for getting the GSE directly.