[BioC] GOstats hyperGTest question
ivan.borozan at utoronto.ca
ivan.borozan at utoronto.ca
Fri Jan 26 00:37:56 CET 2007
Hi all,
I got following results using hyperGTest(params) with a given list of genes
> summary(hgOver)
GOBPID Pvalue OddsRatio ExpCount Count Size
1 GO:0030185 0.000000e+00 -73.314685 0.02692165 2 1
2 GO:0006067 0.000000e+00 -110.746479 0.05384330 3 2
3 GO:0006069 0.000000e+00 -110.746479 0.05384330 3 2
ect ...
If for example I look at genes that are associated with the first GO
term (i.e GO:0030185) I get:
> probeSetSummary(hgOver)[[1]]
EntrezID ProbeSetID selected
1 3043 144221 0
2 3043 148425 0
3 3043 3108408 0
4 3043 5708746 0
My question is how are Counts (in this case Count = 2) in the above
summary(hgOver) table obtained ?
Looking at probeSetSummary(hgOver)[[1]] I can see one EntrezID
(EntrezID = 3043) and 4 ProbeSetID associated with this particular
node (i.e GO:0030185).
R version 2.4.0 (2006-10-03)
i686-redhat-linux-gnu
locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] "splines" "tools" "methods" "stats" "graphics" "grDevices"
[7] "utils" "datasets" "base"
other attached packages:
xtable Hu19K8Build17102006 GOstats
Category
"1.4-3" "1.1.0" "2.0.4"
"2.0.3"
genefilter KEGG RBGL
GO
"1.12.0" "1.14.1" "1.10.0"
"1.14.0"
graph multtest survival
sma
"1.12.0" "1.12.0" "2.29"
"0.5.15"
annotate Biobase human19K601042005
"1.12.0" "1.12.1" "1.1.0"
all the best,
Ivan
More information about the Bioconductor
mailing list