How to get all qualified GEO Platform records in R
1
0
Entering edit mode
Yuan Hao ▴ 240
@yuan-hao-3658
Last seen 9.6 years ago
United States
Dear List, I want to get all the Affymetrix datasets (Human) from GEO and then group them according to different types of chips, like hgu133A, hgu133plus2, etc. so that I could know how many datasets available for each type of chips. I looked at the GEOquery package in R/Bioconductor. There are a 'manufacturer' and a 'Title' field in the GEO Platform record which might be useful for my purpose, but I didn't manage to figure out how to get all the Platform records into R workspace first. In the package vignette, it always deals with a single record by giving a specific accession number such as GPLxxxx, but how could I get all the records meet my criteria? Could you kindly provide some help on this? Thank you very much in advance! kind regards, Yuan
hgu133a GEOquery hgu133a GEOquery • 828 views
ADD COMMENT
0
Entering edit mode
@sean-davis-490
Last seen 3 months ago
United States
On Wed, Sep 1, 2010 at 7:39 AM, Yuan Hao <yuan.x.hao@gmail.com> wrote: > Dear List, > > I want to get all the Affymetrix datasets (Human) from GEO and then group > them according to different types of chips, like hgu133A, hgu133plus2, etc. > so that I could know how many datasets available for each type of chips. > > I looked at the GEOquery package in R/Bioconductor. There are a > 'manufacturer' and a 'Title' field in the GEO Platform record which might be > useful for my purpose, but I didn't manage to figure out how to get all the > Platform records into R workspace first. In the package vignette, it always > deals with a single record by giving a specific accession number such as > GPLxxxx, but how could I get all the records meet my criteria? Could you > kindly provide some help on this? Thank you very much in advance! > > Hi, Yuan. Take a look at the GEOmetadb package. There is an example almost matching the one you describe. To get you started, install the GEOmetadb package and then try: library(GEOmetadb) fn = getSQLiteFile() con <- dbConnect(SQLite(), fn) rs <- dbGetQuery(con, paste("select gpl.bioc_package,gsm.gpl,count(gsm)", "from gsm join gpl on gsm.gpl=gpl.gpl", "where gpl.manufacturer='Affymetrix'", "group by gpl.bioc_package,gsm.gpl")) rs Sean [[alternative HTML version deleted]]
ADD COMMENT

Login before adding your answer.

Traffic: 650 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6