[BioC] AnnBuilder doesn't work on gene information ("GENENAME", "SYMBOL", etc)

Luo Weijun luo_weijun at yahoo.com
Tue Jul 11 06:42:16 CEST 2006


Thank you so much, guys, 

For John¡¯s question, My base map type is "ll", since
my probes are mapped to Entrez Gene IDs directly.
(please all see below for more details)

For Nianhua¡¯s question, below is all my codes for
this work.  

library(AnnBuilder)
load('/Users/luow/project/microarraydata/annotation/hs95av2Entrezg7.Rdata')
myBase=cbind(hs95av2Entrezg7,hs95av2Entrezg7)
myBase[,1]=paste(hs95av2Entrezg7,'_at',sep='')
write.table(myBase,file='/Users/luow/project/microarraydata/annotation/hs95av2Entrezg7Base.txt',sep='\t',row.names=F,col.names=F)
myBase='/Users/luow/project/microarraydata/annotation/hs95av2Entrezg7Base.txt'
myBaseType <- "ll"

myDir <-
'/Users/luow/project/microarraydata/annotation/'
ABPkgBuilder(baseName = myBase, baseMapType =
myBaseType,
pkgName = "hs95av2Entrezg7", pkgPath = myDir,
organism = "Homo sapiens", version = "1.1.0", author =
list(authors = "Weijun",
maintainer = "Weijun <luo_weijun at yahoo.com>"), fromWeb
=T)


Here is part of my base file:
>
a=read.delim('/Users/luow/project/microarraydata/annotation/hs95av2Entrezg7Base.txt',sep='\t',
head=F)
> head(a)                                             
                                                 
        V1    V2
1 10001_at 10001
2 10002_at 10002
3 10003_at 10003
4 10004_at 10004
5 10005_at 10005
6 10006_at 10006

Here is what I got when I created and installed the
package. And Nianhua is right, I got problem with
hs95av2Entrezg7CHRLOC environment, no entry at all. I
will need to update my AnnBuilder. But exactly how can
I do that, using update.packages function or, is there
anything for updating bioc packages specially?

> library(hs95av2Entrezg7)
> hs95av2Entrezg7()


Quality control information for  hs95av2Entrezg7 
Date built: Created: Sat Jul  8 20:23:41 2006  
 
Number of probes: 8359 
Probe number missmatch: None 
Probe missmatch: None 
Mappings found for probe based rda files: 
         hs95av2Entrezg7CHRLOC found 0 of 8359
         hs95av2Entrezg7ENZYME found 1313 of 8359
         hs95av2Entrezg7LOCUSID found 8292 of 8359
         hs95av2Entrezg7PATH found 2784 of 8359 
Mappings found for non-probe based rda files:
         hs95av2Entrezg7CHRLENGTHS found 25
         hs95av2Entrezg7ENZYME2PROBE found 651
         hs95av2Entrezg7ORGANISM found 1
         hs95av2Entrezg7PATH2PROBE found 180
         hs95av2Entrezg7PFAM found 6997
         hs95av2Entrezg7PROSITE found 5451 

Please let me know if you have any ideas, suggestions,
or other questions. I really appreciate all the kind
help you guys offer!
Weijun



--- Nianhua Li <nli at fhcrc.org> wrote:

> Hi, Weijun,
> 
> The humanLLMappings package was generated by
> function map2LL rather than 
> ABPkgBuilder, therefore it is a different case. I
> used ABPkgBuilder 
> function generate a coupld of pakcages recently, and
> didn't see 
> xxGENENAME or xxSYMBOL environments have any
> problem. So, I suspect the 
> problem is specific to your input data. Could you
> please provide more 
> details? e.g. myBaseType and maybe also a small part
> of myBase
> Also, you may want to upgrade AnnBuilder to 1.10.1.
> The UCSC Genome 
> database updated their schema (one table) recently
> for human and mouse 
> data, and AnnBuilder has been updated
> correspondingly. If you see 
> problems with xxCHRLOC environment, please update
> your AnnBuilder.
> 
> thanks
> 
> nianhua
> 
> Nianhua Li
> Computational Biology, PHS, FHCRC
> 
> _______________________________________________
> 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