setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
o1=read.table("RNA-Poly_0htotal.txt", header = T)
o2=read.table("RNA-Poly_2htotal.txt", header = T)
o3=read.table("RNA-Poly_6htotal.txt", header = T)
o4=read.table("RNA-Poly_IFNtotal.txt", header = T)
o5=read.table("RNA-Poly_IL4total.txt", header = T)

setwd("~/giagkas/RNAseq_POLYseq/")
m1=read.table("RNA-Poly_0htotal.txt", header = T)
m2=read.table("RNA-Poly_2htotal.txt", header = T)
m3=read.table("RNA-Poly_6htotal.txt", header = T)
m4=read.table("RNA-Poly_IFNtotal.txt", header = T)
m5=read.table("RNA-Poly_IL4total.txt", header = T)

#0h abund
o=o1
m=m1

o=subset(o, (o$wtRPKM.RNA>=1)&(o$PValue.RNA<0.05)&(abs(o$logFC.RNA)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$PValue.RNA<0.05)&(abs(m$logFC.RNA)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "0h_Abund_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()






#2h abund
o=o2
m=m2

o=subset(o, (o$wtRPKM.RNA>=1)&(o$PValue.RNA<0.05)&(abs(o$logFC.RNA)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$PValue.RNA<0.05)&(abs(m$logFC.RNA)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "2h_Abund_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()




#6h abund
o=o3
m=m3

o=subset(o, (o$wtRPKM.RNA>=1)&(o$PValue.RNA<0.05)&(abs(o$logFC.RNA)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$PValue.RNA<0.05)&(abs(m$logFC.RNA)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "6h_Abund_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()

#IFN abund
o=o4
m=m4

o=subset(o, (o$wtRPKM.RNA>=1)&(o$PValue.RNA<0.05)&(abs(o$logFC.RNA)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$PValue.RNA<0.05)&(abs(m$logFC.RNA)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "IFN_Abund_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()


#IL4 abund
o=o5
m=m5

o=subset(o, (o$wtRPKM.RNA>=1)&(o$PValue.RNA<0.05)&(abs(o$logFC.RNA)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$PValue.RNA<0.05)&(abs(m$logFC.RNA)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "IL4_Abund_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()






#0h poly
o=o1
m=m1

o=subset(o, (o$wtRPKM.RNA>=1)&(o$wtRPKM.POLY>=1)&(o$PValue.POLY<0.05)&(abs(o$logFC.POLY)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$wtRPKM.POLY>=1)&(m$PValue.POLY<0.05)&(abs(m$logFC.POLY)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "0h_Poly_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()



#2h poly
o=o2
m=m2

o=subset(o, (o$wtRPKM.RNA>=1)&(o$wtRPKM.POLY>=1)&(o$PValue.POLY<0.05)&(abs(o$logFC.POLY)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$wtRPKM.POLY>=1)&(m$PValue.POLY<0.05)&(abs(m$logFC.POLY)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "2h_Poly_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()




#6h poly
o=o3
m=m3

o=subset(o, (o$wtRPKM.RNA>=1)&(o$wtRPKM.POLY>=1)&(o$PValue.POLY<0.05)&(abs(o$logFC.POLY)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$wtRPKM.POLY>=1)&(m$PValue.POLY<0.05)&(abs(m$logFC.POLY)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "6h_Poly_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()




#IFNpoly
o=o4
m=m4

o=subset(o, (o$wtRPKM.RNA>=1)&(o$wtRPKM.POLY>=1)&(o$PValue.POLY<0.05)&(abs(o$logFC.POLY)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$wtRPKM.POLY>=1)&(m$PValue.POLY<0.05)&(abs(m$logFC.POLY)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "IFN_Poly_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()



#IL4poly
o=o5
m=m5

o=subset(o, (o$wtRPKM.RNA>=1)&(o$wtRPKM.POLY>=1)&(o$PValue.POLY<0.05)&(abs(o$logFC.POLY)>log2(1.5)))
m=subset(m, (m$wtRPKM.RNA>=1)&(m$wtRPKM.POLY>=1)&(m$PValue.POLY<0.05)&(abs(m$logFC.POLY)>log2(1.5)))
om=subset(o, o$GeneName%in%m$GeneName)
mo=subset(m, m$GeneName%in%o$GeneName)
nrow(om)==nrow(mo)
library(VennDiagram)
setwd("~/giagkas/RNAseq_POLYseq/RNAseq/DEG/old-Orsalia/")
png(filename = "IL4_Poly_OrsVsMArtin.png", width = 800, height = 800)
grid.newpage()
draw.pairwise.venn(area1 = nrow(o), area2 = nrow(m), 
                   cross.area  = nrow(om), category = c("Orsalia", "Martin"), lty = "blank", 
                   fill = c("skyblue", "pink1"))

dev.off()