Hi,

I'm working with a fungus that don't have its genome sequenced. I have
a list of accession numbers from a blast of peptides obtained from a
proteomics experiment, and each number come from a different organism.

For what I read until now, I found how to recover GO terms for
organism specific databases as in:

org.Mm.egGO[mappedLkeys(egids)]

is there a way to do this for a diverse group of accession numbers
from different organisms?

Assuming that I have my list of annotated GO terms the most simple way
to calculate the probability of "enrichment" of, lets say GOterm1
would be by:

phyper(9-1, 26, 825-26, 18, lower.tail=F)

where, 9 is number of times I saw the term on my sample,
26 is the total number of genes inside this classe (associated to
GOterm1)??,
825 is the total number of genes inside all GO categories found?,
and 18 is my number of samples.

If any one knows a reference explaining each step I would appreciate
the reference, so that I could understand in details before using a
package.

thanks in advance for your help

	[[alternative HTML version deleted]]

