Python implementation for enrichment of terms on already annotated genes
0
0
Entering edit mode
7.4 years ago
Benni ▴ 30

I have a matrix like this:

  1. Column = Genes
  2. Column = Cluster gene belongs to
  3. Column = Annotation1
  4. Column = Annotation2
  5. Column = Annotation3

A new matrix would look like this:

  1. Column = Cluster
  2. Column = Annotation Terms and their pValue.

It checks one column and one term at a time. The background contains all rows (all genes) and the enrichment is checked for all genes in the same cluster.

I need a python modul, that takes those two datasets and calculates the pValue (with benjamini correction). Is there such a modul and which method should I use for the anrichment analysis?

Gene Set Enrichment • 2.2k views
ADD COMMENT
0
Entering edit mode

For over enrichment using KEGG of clusters I use functions in the clusterprofiler package, I am pretty sure it would be easier to adapt this to your needs than write it out in Python. There is also a similar function in the Limma package for simple over enrichment testing.

ADD REPLY

Login before adding your answer.

Traffic: 1833 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6