5B

BIN_DATA = read.table("/Ets1/ProcessedData/Ets1_p300_binDATA.txt",sep = "\t", header=T,fill = TRUE)

Qlevel = c("Q10","Q9","Q8","Q7","Q6","Q5","Q4","Q3","Q2","Q1")
Qlevel = as.data.frame(Qlevel)
Qlevel$ratio_DIFF = 0
p300_Ets1_lineDF = Qlevel
label = c("J_Q10","I_Q9","H_Q8","G_Q7","F_Q6","E_Q5","D_Q4","C_Q3","B_Q2","A_Q1")

label = as.data.frame(label)
p300_Ets1_lineDF = cbind(p300_Ets1_lineDF,label)
p300_Ets1_lineDF$group = "ALL"

for(i in 1:nrow(p300_Ets1_lineDF)){

p300_Q10 = BIN_DATA[which(BIN_DATA$p300_Q == p300_Ets1_lineDF[i,1]),]
p300_Ets1_Q10 = p300_Q10[which(p300_Q10$Ets1_Q == p300_Ets1_lineDF[i,1]),]

A = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q1"),])/nrow(p300_Ets1_Q10)*100,2)
B = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q2"),])/nrow(p300_Ets1_Q10)*100,2)
C = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q3"),])/nrow(p300_Ets1_Q10)*100,2)
D = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q4"),])/nrow(p300_Ets1_Q10)*100,2)
E = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q5"),])/nrow(p300_Ets1_Q10)*100,2)

F = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q6"),])/nrow(p300_Ets1_Q10)*100,2)

G = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q7"),])/nrow(p300_Ets1_Q10)*100,2)
H = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q8"),])/nrow(p300_Ets1_Q10)*100,2)
I = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q9"),])/nrow(p300_Ets1_Q10)*100,2)
J = round(nrow(p300_Ets1_Q10[which(p300_Ets1_Q10$DIFF_Q == "Q10"),])/nrow(p300_Ets1_Q10)*100,2)

p300_Ets1_lineDF[i,2] = (G+H+I+J)/(A+B+C+D+E)
    }

p300_Ets1_lineDF[,1] = c("EP10","EP9","EP8","EP7","EP6","EP5","EP4","EP3","EP2","EP1")

options(repr.plot.width = 5, repr.plot.height = 3, repr.plot.res = 1000, repr.plot.pointsize = 40)
ggplot(p300_Ets1_lineDF,aes(x=factor(label), y=ratio_DIFF), color="black", size=0.2) +
geom_hline(yintercept=1,size=0.2)+
geom_line(group="ALL") +
geom_point(shape=21,color="black", fill="white", size=2) +
scale_x_discrete(labels = c("EP1","EP2","EP3","EP4","EP5","EP6","EP7","EP8","EP9","EP10"))+ 
theme_classic(base_size = 10) +
ggtitle(NULL)+
xlab("Ets1 & p300 binding level") +
ylab("Gap of toward\nAB compartment signal\n(toward A/toward B)") +
ylim(0,2)+
theme(plot.title=element_text(face="bold",hjust=0.5,size=25,color = "black"),
      axis.text.x=element_text(face="bold",size=10),
      axis.text.y=element_text(face="bold",size=10),
      axis.title.x = element_text(face="bold",size = 13,color = "black"),
      axis.title.y = element_text(face="bold",size = 13,color = "black"),
      legend.title=element_text(face="bold",size=15), 
      legend.text=element_text(face="bold",size=15))

Image