I have a list of 291 genes which I had enriched for pathway analysis using the "enrichPathway" function of ReactomePA. I then visualised the enrichment as pathways using the dotplot, emapplot and cnetplot functions.
In the emapplot visualisation, I get broadly two distinct connected sets of pathways (mRNA degradation and cytokine signalling). For aesthetic and publication purposes, I require only one of these connected pathway sets to be shown (cytokine signalling) by literally cutting out the other pathways that I am not interested in (mRNA degradation). I tried to do this by:
1. Converting the enrichResult object to a data frame and then editing this to contain only the pathways I want to visualise. But, I am unable to convert back the edited data frame into enrichResult object for making an emapplot.
2. Picking out genes from the pathway and repeating enrichment. But, every gene is now split into further pathways which is not what I want.
Is there a way this can be accomplished? Maybe by being able to pick the pathway IDs and making an emapplot?