[BioC] problem with lumiRatAll.db

Flora Fan fanx0038 at umn.edu
Wed Jan 20 16:40:20 CET 2010


Hi,

I am trying to use lumi package to analyze Illumina RatRef-12_V1. 
However, about 30% of the probe IDs don't have annotations in 
lumiRatAll.db_1.6.1. Is it because the mapping didn't pass the threshold 
or am I doing something wrong?

For example, for probe ID ILMN_1355285:
nuID = Niuq09_43JW4K1SOe4
probe seq = GAGGTGGGGTCATTCTTGTGATCTAGCCCGTGAAGGTCCCAGATGCTGTG
manual blast result as well as Illumina annotation file 
RatRef-12_V1_0_R3_11222119_A.txt shows that it mapped to A4galt.

But when I do getSYMBOL, no symbol was returned.
 > getSYMBOL("Niuq09_43JW4K1SOe4","lumiRatAll.db")
Niuq09_43JW4K1SOe4
                NA
or
 >  x<- lumiRatAllACCNUM
 >  mapped_probes <- mappedkeys(x)
 >  ind<-c(1:16134)[mapped_probes == "Niuq09_43JW4K1SOe4"]
It returned integer(0)

How can I annotate these IDs? Can I keep the original gene symbols that 
were generated from GenomeStudio version 1.0.6 as identifiers in lumi 
package?

Here is the sessionInfo()
R version 2.9.1 (2009-06-26)
i386-apple-darwin8.11.1

locale:
en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base    

other attached packages:
 [1] lumiRatAll.db_1.6.1    lumiRatV1_1.2.0        beadarray_1.12.1     
 [4] hwriter_1.1            lumiRatIDMapping_1.2.1 limma_2.18.0         
 [7] lumi_1.10.2            RSQLite_0.7-1          DBI_0.2-4            
[10] preprocessCore_1.6.0   mgcv_1.5-5             affy_1.22.0          
[13] annotate_1.22.0        AnnotationDbi_1.6.0    Biobase_2.4.1        

loaded via a namespace (and not attached):
[1] affyio_1.12.0   grid_2.9.1      lattice_0.17-25 nlme_3.1-92   
[5] tools_2.9.1     xtable_1.5-5  

Your help will be greatly appreciated.
-Flora

-- 
Flora 
Biostatistics and Informatics
Masonic Cancer Center, University of Minnesota



More information about the Bioconductor mailing list