tctcab
DEG.brca.edgeR.filt.TSS<-DEG.brca.edgeR$Mycontrast[which(abs(DEG.brca.edgeR$Mycontrast$logFC) >= 6), ]
str(DEG.brca.edgeR.filt.TSS)
#57 genes
TCGAVisualize_volcano(DEG.brca.edgeR$Mycontrast$logFC,
DEG.brca.edgeR$Mycontrast$FDR,
filename = "LuminalABvsNormal_FC6.TSS.edgeR.pdf",
xlab = "logFC",
names = rownames(DEG.brca.edgeR$Mycontrast),
show.names = "highlighted",
x.cut = 1,
y.cut = 0.01,
highlight = rownames(DEG.brca.edgeR$Mycontrast)[which(abs(DEG.brca.edgeR$Mycontrast$logFC) >= 6)],
highlight.color = "orange")