cont cs analysis

chipseq_workflow/cs_region_dba.R

@@ -94,6 +94,11 @@ deResults %>% ggplot(aes(pvalue)) +
 #res <- results(dds, contrast=c(c("condition","H1M_Shield", "H3K4_Shield")))
 
 deResults %<>% mutate(is_hit=pvalue<0.01)
+
+write.delim(deResults, file="deResults.txt")
+# deResults <- read.delim("deResults.txt")
+#'  [deResults](deResults.txt)
+
 degs <- deResults %>% filter(is_hit)
 
 ggplot(degs, aes(paste(sample_1, "vs",  sample_2))) + geom_bar() +xlab(NULL) + ylab("# DBGs") +ggtitle("DBG count summary") + coord_flip()