[BioC] missing probe_id in ath1121501.db

Marc Carlson mcarlson at fhcrc.org
Tue Feb 21 22:33:52 CET 2012


Hi Stefan,

I have done some detective work and tracked this down for you.  The 
basic problem here seems to be that for this one probe there is no data 
whatsoever coming from  tair for ANY of the genes that map to it 
(probeset "245651_s_at").  According the sources I have, affy matches 
this probeset matches 5 loci IDs and TAIR does not provide data for any 
of them.  I agree with you though that the probe should still be present 
and will look into a correction for this.

   Marc


On 02/13/2012 12:23 AM, Stefan Bleuler wrote:
> Hi,
>
> I am working with the Bioconductor annotation package for the 
> Arabidopsis thaliana GeneChip ATH1 and have found a strange difference 
> between version 2.7 and versions 2.8 or 2.9.
>
> In 2.7:
> source("http://bioconductor.org/biocLite.R")
> biocLite("ath1121501.db")
> library("ath1121501.db")
> length(as.list(toggleProbes(ath1121501ACCNUM,"all")))
>
> produces the following output:
> [1] 22810
>
> In 2.8 or 2.9 the same code produces the following output:
> [1] 22809
>
> Comparing the lists of probe_ids shows that the probe_id 245651_s_at 
> is present in 2.7 and not present in 2.8 or 2.9. I know that the 
> annotations for a given probe_id can change with updates but why would 
> a probe_id be removed? This probe_id is still present in the latest 
> annotation files provided by Affymetrix.
>
> Thanks for any idea.
>
> Stefan Bleuler
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> 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