Question: ClusterProfiler dotplot method does not include all input groups.
0
gravatar for markm014
6 weeks ago by
markm0140
markm0140 wrote:

I am using the clusterProfiler methods enrichGO and dotplot to do some downstream differential expression analysis. I invoke enrichGO with the following call:

formula_res = compareCluster(`EntrezGene ID~Direction+Treatment, data=filtered_data, fun="enrichGO", OrgDb="org.Hs.eg.db", ont="BP", qvalueCutoff= 0.2, readable=TRUE)

where filtered_data is a dataframe with the follwoing columns:

Gene, EntrezGene ID, Gene Stable ID, Name, Treatment, Direction

Then, plot via the call to dotplot:

p = dotplot(formula_res, showCategory=10, title="test", font.size=9, includeAll = TRUE)

I am surprised that if I have in my data frame 8 unique combinations of treatment and direction (Ex: four different treatments with both up and down regulated genes for each) my dotplot often outputs less groups. For instance, I often get a plot with six groups shown despite having eight "input groups" (as defined by unique combination of direction + treatment).

I have read the documentation and tried setting showCategory to NULL instead of 10, which the documentation says should force a return of all groups, but still not all treatments are included. Am I missing something obvious on this one? I want to verify my results are not erroneous.

ADD COMMENTlink written 6 weeks ago by markm0140

If it doesn't find any enriched terms for a group that group probably doesn't get added to the result. Try running enrichGO on some of the groups that didn't appear in your results.

ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by rpolicastro1.7k

Thank you for taking the time to answer. The thing is, when I run the enrichGo function on each of the 8 input groups separately with the same parameters, they all 8 have significant GO terms. Sorry, I should have included this above!

ADD REPLYlink written 5 weeks ago by markm0140
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 740 users visited in the last hour