egeulgen / pathfindR

pathfindR: Enrichment Analysis Utilizing Active Subnetworks
https://egeulgen.github.io/pathfindR/
Other
178 stars 25 forks source link

Unchanged genes marked as light green in pathway visualization #119

Closed B-McBean closed 2 years ago

B-McBean commented 2 years ago

Hello! This is a really great tool and the documentation is excellent! I am having a small problem when I'm using the wrapper function run_pathfindR (and the visualize_hsa_KEGG function by itself). Genes that do not have changed expression in the pathway are marked as light green, which makes it looks like they have changed expression.

I'm wondering if this is a problem others have experienced and if there's any way to make unchanged genes (or rather, genes not even included on the significantly changed genes list) white instead?

Thanks!

(I've attached pictures to show what I mean in case it isn't clear. In the one with blue and orange, I used the visualize_hsa_KEGG function and changed the low and high to blue and orange (node_cols = c("blue", "gray", "orange")) but there's still green on here for genes in the pathway whose expression does not change. The same is occurring when I use run_pathfindR, which is the other image with the default red and green for high and low.)

hsa04066_pathfindR

hsa04066_pathfindR

egeulgen commented 2 years ago

Hello, This happens with large (>60 genes) KEGG pathways. For these we have to disregard the white color for non-input genes, otherwise KEGG API complains. I'll try to think of a workaround and update you here if I can find one. Best, -E