1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
|
library(tidyverse)
com <- read.csv("common.csv")
com$DEG <- ifelse(com$RNA_log2FC > 1 & com$H3K9ac_log2FC > 1, "UP-UP", ifelse(com$RNA_log2FC < -1 & com$H3K9ac_log2FC < -1, "DOWN-DOWN", ifelse(com$RNA_log2FC > 1 & com$H3K9ac_log2FC < -1, "UP-DOWN", ifelse(com$RNA_log2FC < -1 & com$H3K9ac_log2FC > 1, "DOWN-UP", "NONE"))))
ggplot(com, aes(x = RNA_log2FC, y=H3K9ac_log2FC, colour=DEG)) + geom_point(alpha=0.85, size=1) + geom_vline(xintercept=c(-1,1),lty=4,col="black",lwd=0.5) + geom_hline(yintercept=c(-1,1),lty=4,col="black",lwd=0.5) + scale_color_manual(values=c('steelblue','#03BF7D','gray', '#AEAF31','brown')) + theme_bw()+ guides(color = "none") + theme(plot.margin = unit(c(1,1,1,1), "cm"), panel.grid = element_blank()) + labs(x="RNA-Seq_log2FC")
ggsave("H3K9ac_com.pdf", width = 15, height = 15, units = "cm" )
|