Update 01-limma_analysis.R

d3b-center · Oct 9, 2024 · abfde6e · abfde6e
1 parent a3bbb46
commit abfde6e
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/analyses/methylation_analysis/01-limma_analysis.R b/analyses/methylation_analysis/01-limma_analysis.R
@@ -23,6 +23,7 @@ output_dir <- opt$output_dir
 dir.create(output_dir, showWarnings = F, recursive = T)
 
 # read m-values
+cat('Reading m-values')
 methyl_m_values_full <- readRDS(opt$methyl_mat) %>% 
   dplyr::slice_head(n = 500000) %>%
   na.omit() %>%
@@ -31,6 +32,7 @@ methyl_m_values_full <- methyl_m_values_full %>%
   tibble::column_to_rownames('Probe_ID')
 
 # read annotation
+cat('Reading methylation feature annotations')
 methyl_annot_full <- data.table::fread(opt$methyl_annot) %>%
   dplyr::filter(!duplicated(Probe_ID))