Entering edit mode
7.1 years ago
Björn
▴
110
I run following command in edgeR
tfit<-treat(vfit, lfc = 2)
dt<-decideTests(tfit)
summary(dt) # using log fold change of 2
DM RRMB RRMC RRMP TR
-1 254 0 0 0 25
0 523 792 792 792 531
1 15 0 0 0 236
de.common<-which(dt[,1]!=0 & dt[,2]!=0)
length(de.common)
head(tfit$genes[de.common], n=20)
vennDiagram(dt[,1:3], circle.col=c("blue", "salmon", "green")) # between DM and RRMB
Now my question how can I create venn diagram ONLY between DM (dt:1) and TR (dt:5). Thanks