Difference between revisions of "Using Bioconductor To Analyse Microarray Data"

From Bridges Lab Protocols
Jump to: navigation, search
m
m (Obtaining GEO Datasets)
Line 24: Line 24:
 
**measurements - '''GSM'''
 
**measurements - '''GSM'''
 
**platforms - '''GPL'''
 
**platforms - '''GPL'''
**series - '''GSE''')
+
**series - '''GSE'''
 
<pre>
 
<pre>
 
gds <- getGEO("GDS162")  #load GDS162 dataset
 
gds <- getGEO("GDS162")  #load GDS162 dataset
 
eset <- GDS2eset(gds)  #convert to expression set
 
eset <- GDS2eset(gds)  #convert to expression set
 
</pre>
 
</pre>
*see [http://www2.warwick.ac.uk/fac/sci/moac/currentstudents/peter_cock/r/geo/] or [http://www.bioconductor.org/packages/1.8/bioc/html/GEOquery.html] for more information.
+
*see [[http://www2.warwick.ac.uk/fac/sci/moac/currentstudents/peter_cock/r/geo/ Peter Cock's Page]] or [[http://www.bioconductor.org/packages/1.8/bioc/html/GEOquery.html GEOquery Documentation]] for more information.

Revision as of 23:59, 26 July 2009


Software Requirements

  • R, get from [CRAN]
  • Bioconductor, get from [Bioconductor]
  • Bioconductor packages. Install as needed:
    • Biobase
    • GEOquery - [1]
source("http://www.bioconductor.org/biocLite.R")
biocLite("PACKAGE")

Obtaining GEO Datasets

  • Open a R terminal
  • Load Biobase and GEOquery packages
libary(Biobase)
library(GEOquery)
  • Can load:
    • datasets - GDS
    • measurements - GSM
    • platforms - GPL
    • series - GSE
gds <- getGEO("GDS162")  #load GDS162 dataset
eset <- GDS2eset(gds)  #convert to expression set