[BioC] subsetting an eset object with a defined list

Celine Carret ckc at sanger.ac.uk
Thu Oct 5 12:27:51 CEST 2006


Hi,
I'm stuck remembering the way of subsetting an eset object (19 affy 
arrays) to accomodate only a defined list of genes that I'm interested in.
I read.table(list) into R, transformed it as.matrix and use the functions:
exprs(eset)[list] or
exprs(eset)[as.character(list)]

those functions recognize that I have 904 genes in that list, but gives 
me 904 NA
What I would like to obtain is my 904 genes with the values for the 19 
arrays, in order to assess the DE only on this subgroup as the other 
genes are not supposed to be expressed at this stage of the cultures I'm 
looking at.
Could someone please point me in the right direction?
Thank you very much
Celine



sessionInfo()
Version 2.3.0 (2006-04-24)
ia64-unknown-linux-gnu

attached base packages:
[1] "splines"   "tools"     "methods"   "stats"     "graphics"  "grDevices"
[7] "utils"     "datasets"  "base"    

other attached packages:
affycoretools       GOstats      Category      hgu95av2          KEGG
      "1.4.0"       "1.6.0"       "1.4.1"      "1.12.0"      "1.12.0"
     multtest    genefilter      survival        xtable          RBGL
     "1.10.2"      "1.10.1"        "2.24"       "1.3-2"       "1.8.1"
     annotate            GO         graph         Ruuid         limma
     "1.10.0"      "1.12.0"      "1.10.6"       "1.9.2"      "2.7.10"
         affy        affyio       Biobase
     "1.10.0"      "0.99.3"      "1.10.1"


-- 
----------------------------------------------------
Celine Carret PhD
Pathogen Microarrays group
The Wellcome Trust Sanger Institute
Hinxton, Cambridge CB10 1SA, UK.
tel. +44 (0)1223 494 940
fax. +44 (0)1223 494 919
email: ckc at sanger.ac.uk
http://www.sanger.ac.uk/PostGenomics/PathogenArrays/



More information about the Bioconductor mailing list