saezlab / transcriptutorial

This is a tutorial to guide the analysis of RNAseq dataset using footprint based tools such as DOROTHEA, PROGENY and CARNIVAL
https://saezlab.github.io/transcriptutorial/
GNU General Public License v3.0
55 stars 30 forks source link

Enrichment Analysis significance cutoff #21

Closed Semple-bio closed 1 year ago

Semple-bio commented 3 years ago

Good day. There are 2 adj. p-value cutoff used in the 'Enrichment Analysis' part, where the secondary is redundant as the first one is more stringent.

First cutoff, at 0.001 level:

#data for plotting
PathwaysSelect <- sig_pathways_df %>%
    dplyr::select(pathway, `p-value`, `Adjusted p-value`) %>%
    dplyr::filter(`Adjusted p-value` <= 0.001) %>%
    dplyr::rename(pvalue = `p-value`, AdjPvalu = `Adjusted p-value`) %>% 
    dplyr::mutate(pathway = as.factor(pathway))
Second cutoff, with no effect:
ggdata = PathwaysSelect %>% 
    dplyr::filter(AdjPvalu <= 0.05) %>% 
    dplyr::group_by(source) %>% 
    dplyr::arrange(AdjPvalu) %>%
    dplyr::slice(1:5)