Skip to content

Commit

Permalink
Update 01-dge_analysis_deseq.R
Browse files Browse the repository at this point in the history
  • Loading branch information
aadamk authored Oct 9, 2024
1 parent 69824aa commit f0c536f
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions analyses/dge_pathway_analysis/01-dge_analysis_deseq.R
Original file line number Diff line number Diff line change
Expand Up @@ -92,7 +92,7 @@ for (i in 1:length(clusters)) {
output_df <- rbind(output_df, deseq_output)

# pathway enrichment using REACTOME
cat('Running Reactome pathway enrichment')
cat('Running Reactome pathway enrichment \n')
reactome_pathways <- msigdbr::msigdbr(category = "C2", subcategory = "CP:REACTOME")
reactome_pathways <- reactome_pathways %>%
dplyr::select(gs_name, gene_symbol) %>%
Expand All @@ -108,7 +108,7 @@ for (i in 1:length(clusters)) {
)

# pathway enrichment using HALLMARK
cat('Running Hallmark pathway enrichment')
cat('Running Hallmark pathway enrichment \n')
hallmark_pathways <- msigdbr::msigdbr(category = "H", subcategory = NULL)
hallmark_pathways <- hallmark_pathways %>%
dplyr::select(gs_name, gene_symbol) %>%
Expand All @@ -125,6 +125,6 @@ for (i in 1:length(clusters)) {
}

# write output to tsv
cat('Writing outputs')
cat('Writing outputs \n')
write_tsv(x = output_df,
file = file.path(results_dir, "diffexpr_output_per_cluster.tsv"))

0 comments on commit f0c536f

Please sign in to comment.