### R code from vignette source 'vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw'

###################################################
### code chunk number 1: Userguide.Rnw:38-42
###################################################
library(org.Bt.eg.db)
genes <- c('280705', '280706', '100327208')
symbols <- org.Bt.egSYMBOL[genes]
toTable(symbols)


###################################################
### code chunk number 2: Userguide.Rnw:55-61
###################################################
symbols <- c('BCSE','TF','TG')
entrez <- org.Bt.egSYMBOL2EG[symbols]
entrez <- toTable(entrez)[,'gene_id']
entrez
refseq <- org.Bt.egREFSEQ[entrez]
toTable(refseq)


###################################################
### code chunk number 3: Userguide.Rnw:71-74
###################################################
library(AnnotationFuncs)
translate(c(280705, 280706, 100327208, 123), org.Bt.egSYMBOL)
translate(c('BCSE','TF','TG'), from=org.Bt.egSYMBOL2EG, to=org.Bt.egREFSEQ)


###################################################
### code chunk number 4: Userguide.Rnw:101-108
###################################################
library(AnnotationFuncs)
library(org.Bt.eg.db)
genes <- c(280705, 280706, 100327208)
translate(genes, org.Bt.egGENENAME)
translate(genes, org.Bt.egCHR)
translate(genes, org.Bt.egENSEMBL)
translate(genes, org.Bt.egREFSEQ, remove.missing=F)


###################################################
### code chunk number 5: Userguide.Rnw:115-118
###################################################
symbols <- c("SERPINA1","KERA","CD5")
translate(symbols, org.Bt.egSYMBOL2EG)
translate(symbols, revmap(org.Bt.egSYMBOL))


###################################################
### code chunk number 6: Userguide.Rnw:124-129
###################################################
translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egGENENAME)
# As a curiosity, if you specify another organism, you are warned:
library(org.Hs.eg.db)
translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Hs.egGENENAME)
warnings()


###################################################
### code chunk number 7: Userguide.Rnw:136-140
###################################################
symbols <- c("SERPINA1","KERA","CD5")
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ)
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, reduce='first')
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, reduce='last')


###################################################
### code chunk number 8: Userguide.Rnw:145-146
###################################################
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, return.list=F)


###################################################
### code chunk number 9: Userguide.Rnw:159-165
###################################################
symbols <- c("SERPINA1","KERA","CD5")
refseq <- translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egREFSEQ)
mRNA <- pickRefSeq(refseq, priorities=c('NM','XM'))	# Equals pickRefSeq.mRNA
mRNA
proteins <- pickRefSeq(refseq, priorities=c('NP','XP')) # Equals pickRefSeq.Protein
proteins


###################################################
### code chunk number 10: Userguide.Rnw:177-184 (eval = FALSE)
###################################################
## symbols <- c("SERPINA1","KERA","CD5")
## GOs <- translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egGO)
## # Pick biological process:
## pickGO(GOs, category='BP')
## # Pick only those biological processes for Entrez Gene ID 280730
## # which have been inferred from sequence similarity or electronic annotation:
## pickGO(translate(280730, org.Bt.egGO), category='BP', evidence=c('ISS','IEA'))


###################################################
### code chunk number 11: Userguide.Rnw:197-198
###################################################
toLatex(sessionInfo())