[BioC] Many missing GO terms from GOHyperG call using mgu74av2
John Zhang
jzhang at jimmy.harvard.edu
Wed Jun 29 17:53:56 CEST 2005
>In the course of using GOHyperG, I checked my results against results from
NETAFFX and found large discrepancies in the number of GO terms returned for
some probeids.
Unless we build the annotation packages frequently (every week?), there will be
descripancies. You may try to build your own annotation packages to minimize the
descripancies.
>
>Here is an example code snippet that illustrates my point:
>
>library(affy)
>library(mgu74av2)
>library(GOstats)
>library(annaffy)
>smallList <- c("160102_at")
>myLL <- unlist(mget(smallList, mgu74av2LOCUSID))
>sum(duplicated(myLL))
>length(myLLunique<- !duplicated(mget(smallList, mgu74av2LOCUSID)))
>bphyper <- GOHyperG(myLL[myLLunique], lib="mgu74av2", what="BP")
>
>sort(names(bphyper$go2Affy))
>[1] "GO:0006457" "GO:0007582" "GO:0008150" "GO:0008152" "GO:0009987"
"GO:0019538"
>[7] "GO:0043170" "GO:0050875"
>
>#From NETAFFX for 160102_at
>#"GO:0006457" "GO:0007582" "GO:0008150" "GO:0008152" "GO:0009987" "GO:0019538"
>#"GO:0043170" "GO:0050875"
>#"GO:0044238" "GO:0044237" "GO:0044260" "GO:0044260" "GO:0044267
>
>Notice that NETAFFX returns 5 additional GO terms that are missing from the
results from the GOHyperG call.
>
>Here is my sessionInfo:
>
>> sessionInfo()
>R version 2.1.0, 2005-04-18, i386-pc-mingw32
>
>attached base packages:
>[1] "splines" "tools" "methods" "stats" "graphics"
>[6] "grDevices" "utils" "datasets" "base"
>
>other attached packages:
> annaffy KEGG hgu95av2 GOstats multtest genefilter
> "1.0.18" "1.6.8" "1.8.4" "1.1.2" "1.6.0" "1.6.1"
> survival xtable RBGL annotate GO XML
> "2.17" "1.2-5" "1.3.7" "1.5.15" "1.6.8" "0.95-6"
> graph Ruuid cluster mgu74av2 affy reposTools
> "1.5.0" "1.5.0" "1.9.8" "1.8.4" "1.6.7" "1.5.19"
> Biobase
> "1.5.12"
>
>Is this a problem with the mgu74av2 metadata, or with GOHyperG?
>
>Thanks much,
>Dick
>*******************************************************************************
>Richard P. Beyer, Ph.D. University of Washington
>Tel.:(206) 616 7378 Env. & Occ. Health Sci. , Box 354695
>Fax: (206) 685 4696 4225 Roosevelt Way NE, # 100
> Seattle, WA 98105-6099
>http://depts.washington.edu/ceeh/ServiceCores/FC5/FC5.html
>http://staff.washington.edu/~dbeyer
>
>_______________________________________________
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch
>https://stat.ethz.ch/mailman/listinfo/bioconductor
Jianhua Zhang
Department of Medical Oncology
Dana-Farber Cancer Institute
44 Binney Street
Boston, MA 02115-6084
More information about the Bioconductor
mailing list