Author: | Haibao Tang (tanghaibao), Brent Pedersen (brentp), Aurelien Naldi (aurelien-naldi) |
---|---|
Email: | [email protected] |
License: | BSD |
Contents
This package contains a Python library to
- process over- and under-representation of certain GO terms, based on Fisher's exact test. Also implemented several multiple correction routines (including Bonferroni, Sidak, and false discovery rate).
- process the obo-formatted file from Gene Ontology website. The data structure is a directed acyclic graph (DAG) that allows easy traversal from leaf to root.
Python version >= 2.6, try install this package first. Within this folder:
easy_install .
.obo
file for the most current gene ontology:wget http://geneontology.org/ontology/obo_format_1_2/gene_ontology.1_2.obo
fisher module for calculating Fisher's exact test:
easy_install fisher
And put it in the current folder where you want to do your analysis.
If you need to plot the ontology lineage, you need the following to tools to be installed.
Graphviz, for graph visualization.
pygraphviz, Python binding for communicating with Graphviz:
easy_install pygraphviz
run.sh
contains example cases, which calls the utility scripts in the scripts
folder.
see find_enrichment.py
for usage. It takes as arguments files containing:
- gene names in a study
- gene names in population (or other study if --compare is specified)
- an association file that maps a gene name to a GO category.
please look at data/
folder to see examples on how to make these files. when ready, the command looks like:
python scripts/find_enrichment.py --alpha=0.05 --indent data/study data/population data/association
and can filter on the significance of (e)nrichment or (p)urification. it can report various multiple testing corrected p-values as well as the false discovery rate.
The "e" in the "Enrichment" column means "enriched" - the concentration of GO term in the study group is significantly higher than those in the population. The "p" stands for "purified" - significantly lower concentration of the GO term in the study group than in the population.
see plot_go_term.py
for usage.
plot_go_term.py
can plot the lineage of a certain GO term, by:
python scripts/plot_go_term.py --term=GO:0008135
this will plot the following image.