ported fc analysis to h3ha data

chipseq_workflow/chipseq_utils.sh

@@ -106,14 +106,14 @@ rm algn_counts.txt
 for bamFile in $(ls $*); do
 (
     echo "processing $bamFile"
-    samtools flagstat $bamFile | head -n1 | cut -f1 -d' ' | sed -e 's/$/\t'$(basename $bamFile _mmf.bam)'/' >> algn_counts.txt
+    samtools flagstat $bamFile | head -n1 | cut -f1 -d' ' | sed -e 's/$/\t'$(basename $bamFile .bam)'/' >> algn_counts.txt
 ) &
 done
 wait
 
 echo '
-devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v1.22/R/core_commons.R")
-devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v1.22/R/ggplot_commons.R")
+devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v1.26/R/core_commons.R")
+devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v1.26/R/ggplot_commons.R")
 
 libSizes <- read_tsv("algn_counts.txt", col_names=F) %>% set_names(c("library_size", "sample"))
 libSizes %>% ggplot(aes(sample, library_size)) +