Skip to content

Commit

Permalink
Update run_analysis.sh
Browse files Browse the repository at this point in the history
  • Loading branch information
aadamk committed Sep 20, 2024
1 parent e51b49f commit d1badc7
Showing 1 changed file with 12 additions and 12 deletions.
24 changes: 12 additions & 12 deletions analyses/methylation_analysis/run_analysis.sh
Original file line number Diff line number Diff line change
Expand Up @@ -28,22 +28,22 @@ Rscript --vanilla 01-limma_analysis.R \

# 2) gsameth pathway enrichment on differentially expressed CpG sites from (genebody + promoter) obtained from 01-limma_analysis.R
# using REACTOME
Rscript --vanilla 02-dms_gsameth_analysis.R \
--methyl_mat $methyl_m_file \
--diffexpr_sites "results/limma_output/genebody_promoter_diffexpr_probes_per_cluster.tsv" \
--msigdb "reactome" \
--prefix "genebody_promoter" \
--output_dir "results/dms_gsameth_output/reactome" \
--plots_dir "plots/dms_gsameth_output/reactome"

# using HALLMARK
# Rscript --vanilla 02-dms_gsameth_analysis.R \
# --methyl_mat $methyl_m_file \
# --diffexpr_sites "results/limma_output/genebody_promoter_diffexpr_probes_per_cluster.tsv" \
# --msigdb "hallmark" \
# --msigdb "reactome" \
# --prefix "genebody_promoter" \
# --output_dir "results/dms_gsameth_output/hallmark" \
# --plots_dir "plots/dms_gsameth_output/hallmark"
# --output_dir "results/dms_gsameth_output/reactome" \
# --plots_dir "plots/dms_gsameth_output/reactome"

# using HALLMARK
Rscript --vanilla 02-dms_gsameth_analysis.R \
--methyl_mat $methyl_m_file \
--diffexpr_sites "results/limma_output/genebody_promoter_diffexpr_probes_per_cluster.tsv" \
--msigdb "hallmark" \
--prefix "genebody_promoter" \
--output_dir "results/dms_gsameth_output/hallmark" \
--plots_dir "plots/dms_gsameth_output/hallmark"

# 2) methreg analysis on differentially methylated CpG sites from (genebody + promoter) obtained from 01-limma_analysis.R
# Rscript --vanilla 02-methreg_analysis.R \
Expand Down

0 comments on commit d1badc7

Please sign in to comment.