[BioC] Probe IDs from Sentrix chips not appearing in Human annotation file
Marc Carlson
mcarlson at fhcrc.org
Wed Apr 29 20:07:49 CEST 2009
Hi Nikhil,
What do you have installed? Can you give us the output of sessionInfo()?
Since you are using bioconductor, you might want to have a look at the
annotation package that is maintained for this platform:
http://www.bioconductor.org/packages/devel/data/annotation/html/illuminaHumanv3.db.html
If you want some help in how to use that package, you might also look at
the AnnotationDbi package:
http://www.bioconductor.org/packages/devel/bioc/html/AnnotationDbi.html
And finally, if you should discover that your HT platform is
significantly different from what has been offered by illumina before,
you could even build a custom annotation package by following the
instructions in the SQLForge vignette (which you will also find along
with AnnotationDbi).
Marc
Nikhil Joshi wrote:
> Hi all,
>
> This is my first post, so please bear with me. I was not able to find
> any information regarding this in the archives, so I am posting here.
> I am doing a differential expression analysis using three Human
> Sentrix chips in R, utilizing the beadarray package. I am then using
> the Human annotation file from Illumina (i.e.
> HumanHT-12_V3_0_R2_11283641_A.txt) to match Probe
> IDs/Array_address_IDs to GI numbers. So here are my questions:
>
> 1. In R, when I load my data using readIllumina
> (backgroundMethod="rma", normalizeMethod="quantile",
> annoPkg="Humanv3"), the data that I get has data for each of the
> arrays, and within an array there is a column called "ProbeID".
> However, this ProbeID seems to correspond to the Array_address_ID
> field in the annotation file, even though there is a field in the
> annotation file called Probe_Id! Can anyone explain why the ProbeID
> in R does NOT correspond to the Probe_Id field in the annotation file,
> but instead corresponds to the Array_address_Id field? Am I looking
> at this incorrectly?
>
> 2. After I do my differential expression analysis, I get data for my
> top results. These results include a ProbeID field. Most of these
> ProbeIDs have entries in the annotation file, however, about 300 do
> not. I looked for the probes in the R data, and I found that they did
> exist and the replication was just as much as all of the other
> probes. Which leads me to believe that they are actual probes, and
> not some error of some kind. So my question is, why don't these
> probes have annotation information associated with them in the
> annotation file?
>
> Any help is highly appreciated. Thank you!
>
> - Nik.
>
More information about the Bioconductor
mailing list