[BioC] unexpected change in number of exons from/on mouse exon arrays
Michal Okoniewski
michal.okoniewski at fgcz.ethz.ch
Wed Jun 25 15:31:43 CEST 2008
Hi Wolfgang,
I do not remember changing ever the CDFs for exon chips.
1411189 is the number of probesets for human chips.
1201338 is probably the right number for mouse. For rat it is some 1.1M
So most likely you must have been running human CDF before...
It is a bit like the joke (sorry, I could not help it ;) ):
Teacher: ... and bear in mind that the boiling point for water is 90
degrees...
Student: But professor, we were already told that the water boils in 100
degrees!
T (checking his notes): Right, 90 degrees is the right angle...
all the best,
Michal
Wolfgang RAFFELSBERGER wrote:
> Dear Bioconductors,while re-running the basic rma protocol on some Affy mouse exon arrays (Mouse Exon 1.0 alias MoEx-1.0) I noticed that all at sudden I get 1201338 exons(/elements) instead of 1411189 in Oct 2007. Since the R and the BioC release(s) changed, however, the annotation package (on www.bioconductor.org) is still the same (i.e. mouseexonpmcdf_1.1).a) Is this correct that mouse Affy exon 1.0 arrays now give just 1201338 exons ?b) How can one explain this change if the annotation underneith hasn't changed ?Here the code I use :library(exonmap)library(affy)dat.raw <- read.exon() dat.raw at cdfName <- "mouseexonpmcdf" eset.rma <- rma(dat.raw) nrow(exprs(eset.rma)) # get (june 2008) : 1201338 ; BUT was (oct 2007) 1411189 eset.rma # gives :ExpressionSet (storageMode: lockedEnvironment)assayData: 1201338 features, 24 samples element names: exprs phenoData sampleNames: A1xx.CEL, A2xx.CEL, ..., Z2xx.CEL (24 total) varLabels and varMetadata description: sample: NA treatment: NAfeatureData featureNames: 4304920, 4304921, ..., 5617580 (1201338 total) fvarLabels and fvarMetadata description: noneexperimentData: use 'experimentData(object)'Annotation: moex10stv1 sessionInfo() # gives :R version 2.7.0 (2008-04-22) x86_64-unknown-linux-gnu locale:Cattached base packages:[1] splines tools stats graphics grDevices utils datasets [8] methods base other attached packages: [1] mouseexonpmcdf_1.1 exonmap_1.6.02 plier_1.10.0 [4] RMySQL_0.6-0 DBI_0.2-4 RColorBrewer_1.0-2 [7] simpleaffy_2.14.05 gcrma_2.12.1 matchprobes_1.12.0 [10] genefilter_1.20.0 survival_2.34-1 affy_1.16.0 [13] preprocessCore_1.0.0 affyio_1.6.1 Biobase_2.0.1 loaded via a namespace (and not attached):[1] AnnotationDbi_1.0.6 RSQLite_0.6-8 annotate_1.18.0 Thanks for any helpful comments,Wolfgang . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .Wolfgang Raffelsberger, PhDLaboratoire de BioInformatique et Génomique IntégrativesIGBMC1 rue Laurent Fries, 67404 Illkirch Strasbourg, FranceTel (+33) 388 65 3300 , Fax (+33) 388 65 3276http://www-bio3d-igbmc.u-strasbg.fr/~wraffwolfgang.raffelsberger@igbmc.fr
>
> [[alternative HTML version deleted]]
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list