2
$\begingroup$

So I'm doing a pathway analysis using GeneTrail database. One thing I am not able to get is this similarity measure being calculated and plotted:

figure

I want to know how this is being made, what is the output being plotted because each gene is may be involved in multiple pathways.

I want to be able to make it in R again. Any suggestion or help would be appreciated. How to plot the enrichment result in this form ?

$\endgroup$
3
$\begingroup$

This is a a symmetric matrix with the pathways in both axes. So for each pair a similarity index is calculated (usually it takes into account the number of genes in each pathway and the total number of genes in both pathways). There are several similarity indexes available. A typical one is Jaccard or Sørensen-Dice index, usually 0 means no similarity and 1 equal pathways.

You can do this kind of analysis in R using BioCor (Disclosure: I am the author of this package) which uses the Dice similarity. You will need the name of the genes of each pathway, after calculating the similarity of pathways using mpathSim, you can then plot it with a heatmap (using the heatmap function).

| improve this answer | |
$\endgroup$
  • $\begingroup$ yes i remember from my last post you have given the link im going to use it and let you know ..i think this looks bit more comprehensive to put all the pathway together..if i have any issue with your library i will ask you.. $\endgroup$ – krushnach Chandra Jan 29 '19 at 8:59

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.