Make pan-Epi and pan-St DEG lists

r.scripts/sc_DEG.R

@@ -180,6 +180,9 @@ DotPlot(sc10x.Group.NoStress,DEG,x.lab.rot=TRUE,plot.legend=TRUE,dot.scale=4)
 dev.off()
 rm(DEG)
 
+DEG_Epi <- setdiff(rownames(DEG.Epi.vs.St),c(DEG_BE,DEG_LE,DEG_OE1,DEG_OE2,DEG_NE))
+DEG_St <- setdiff(rownames(DEG.St.vs.Epi),c(DEG_Fib,DEG_SM,DEG_Endo))
+
 for (i in ls(pattern="^DEG")){
   write.table(get(i),file=paste0("./analysis/",opt$g,"/global/DEG/",i,".csv"),row.names=TRUE,col.names=NA,append=FALSE,sep=",")
   write.table(get(i) %in% surface,file=paste0("./analysis/",opt$g,"/global/DEG/",i,".surface.csv"),row.names=TRUE,col.names=NA,append=FALSE,sep=",")