I use GO::TermFinder and want to provide a list of genes for the background population instead of using the default one, does anybody known how to do that?
I use GO::TermFinder and want to provide a list of genes for the background population instead of using the default one, does anybody known how to do that?
You can use this (information from http://search.cpan.org/~sherlock/GO-TermFinder-0.86/lib/GO/TermFinder.pm):
Usage (use a subset of genes as the background population):
my $termFinder = GO::TermFinder->new(annotationProvider=> $annotationProvider,
ontologyProvider => $ontologyProvider,
population => \@genes,
aspect => <P|C|F>);
What is important for you is the "population"
"This argument allows a client to indicate the population that should used to calculate a background distribution of GO terms. In the absence of population argument, then the background distribution will be drawn from all genes in the annotationProvider. This should be provided as an array reference, and no ambiguous names should be provided (see AnnotationProvider for details of name ambiguity). This option is particularly pertinent in a case where for example you assayed only 2000 genes in a two hybrid experiment, and found 20 interesting ones. To find significant terms, you need to do it in the context of the genes that you assayed, not in the context of all genes with annotation."
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Thank you very much, and it really works well!
Great! It is puzzling me now.