update calibration image

Author: stemangiola

dev/simulate_read_counts.R

@@ -4,6 +4,21 @@ library(ppcSeq)
 library(furrr)
 plan(multicore)
 
+my_theme =
+	theme_bw() +
+	theme(
+		panel.border = element_blank(),
+		axis.line = element_line(),
+		panel.grid.major = element_line(size = 0.2),
+		panel.grid.minor = element_line(size = 0.1),
+		text = element_text(size=12),
+		legend.position="bottom",
+		aspect.ratio=1,
+		strip.background = element_blank(),
+		axis.title.x  = element_text(margin = margin(t = 10, r = 10, b = 10, l = 10)),
+		axis.title.y  = element_text(margin = margin(t = 10, r = 10, b = 10, l = 10))
+	)
+
 FDR_threshold = 0.2
 
 res_1 =
@@ -66,7 +81,17 @@ es =
 				 	)
 		)
 
-(es %>% mutate(fp = fp %>% divide_by(10)) %>% unnest(`false positive predicted`) %>% ggplot(aes(x=fp, y= `false positive predicted`)) + geom_jitter() + xlab("false positive") + my_theme) %>% ggsave(plot=., "dev/false_positive_study.pdf", device="pdf")
+(
+	es %>%
+		mutate(fp = fp %>% divide_by(10)) %>%
+		unnest(`false positive predicted`) %>%
+		ggplot(aes(x=fp, y= `false positive predicted`)) +
+		geom_jitter(width = 0.01) +
+		geom_smooth(method = "lm") +
+		xlab("false positive") +
+		my_theme
+) %>%
+	ggsave(plot=., "dev/false_positive_study.pdf", device = cairo_pdf, width=89, height = 90, units = "mm" )
 
 save(list=c("es", "input_2", "res_1"), file="dev/false_positive_study_3_runs.RData")