[BioC] GEOquery, Limma Annotations
    Ovokeraye Achinike-Oduaran 
    ovokeraye at gmail.com
       
    Tue Oct  4 08:18:19 CEST 2011
    
    
  
Hi all,
Can I please get some help with including a code to the script below
to get the actual annotation for the resulting gene list not just the
platform clone_id?
It will be greatly appreciated. Thanks.
Regards,
Avoks
gds157dat = getGEO('GDS157',destdir=".")
gds157eset = GDS2eSet(gds157dat, do.log2=TRUE)
m = pData(gds157eset)$metabolism
design_gds157 = createDesignMatrix(gds157eset)
design_gds157 = model.matrix(~m)
fit = lmFit(gds157eset, design_gds157)
fit2 = eBayes(fit)
results = topTable(fit2, adjust ="BH", number = nrow(gdseset))
excel = write.table(results, file = "C:/Documents and Settings/xxx)
> sessionInfo()
R version 2.13.1 (2011-07-08)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=English_xxx. LC_CTYPE=English_xxx
[3] LC_MONETARY=English_xxx LC_NUMERIC=C
[5] LC_TIME=English_xxx
attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base
other attached packages:
[1] puma_2.4.0      mclust_3.4.10   affy_1.30.0     limma_3.8.3
[5] GEOquery_2.19.4 Biobase_2.12.2
loaded via a namespace (and not attached):
[1] affyio_1.20.0         preprocessCore_1.14.0 RCurl_1.6-10.1
[4] tools_2.13.1          XML_3.4-2.2
    
    
More information about the Bioconductor
mailing list