clusterProfiler enrichGO on only T-cell related GO terms
1
0
Entering edit mode
Asma • 0
@asma-24876
Last seen 5 hours ago
Canada

I have a set of differentially expressed (DE) genes for each cluster calculated using Seurat. Now, I want to identify the enriched GO terms per cluster related to T-cell function. To do this, I extracted three lists (one per category) from the AMIGO database, each containing genes and associated GO terms specific to T-cell function. How can I use these three gene lists as a database in compareCluster with enrichGO, instead of using org.DB?

clusterProfiler • 64 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 hours ago
United States

You might be able to do that, but in general I would doubt that the author of the package included functionality for a very specific edge case like that. The easy way would be to just use the correct OrgDb package and have enrichGO handle all the complexities for you, after which you could just filter to the three terms you care about.

Alternatively, if you want to 'roll yer own', do note that all any of these GO packages do is compute a hypergeometric test using phyper. If you have the set of genes in each term, it should be trivial to compute the required values for phyper - see the help page - and get the p-values for the three terms.

Login before adding your answer.

Traffic: 354 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