[Bioc-devel] add function to "annotate" package
Robert M. Flight
rflight79 at gmail.com
Wed Jul 13 20:47:08 CEST 2011
Would it make sense to add one extra level of convenience to the
"getAnnMap" function in "annotate", that represents what I am sure is
one of the largest number of uses of the bimaps as an environment,
which is getting information about a particular set of genes or other
identifier?
getAnnotation <- function(id,annPackage,mapID,doUnlist=TRUE){
annEnv <- getAnnMap(mapID, annPackage, load=TRUE)
tmpDat <- mget(id, annEnv, ifnotfound=NA)
if (doUnlist){
tmpDat <- unlist(tmpDat, use.names=FALSE)
}
return(tmpDat)
}
# simple example, but I do this for a lot of different things frequently
library(annotate)
testEntrez <- c("1","10","100","1000")
getAnnotation(testEntrez,"org.Hs.eg.db","SYMBOL")
Thanks,
-Robert
Robert M. Flight, Ph.D.
University of Louisville Bioinformatics Laboratory
University of Louisville
Louisville, KY
PH 502-852-1809 (HSC)
PH 502-852-0467 (Belknap)
EM robert.flight at louisville.edu
EM rflight79 at gmail.com
Williams and Holland's Law:
If enough data is collected, anything may be proven by
statistical methods.
More information about the Bioc-devel
mailing list