R语言之可视化（33）绘制差异基因分析统计图

GSE1 <- read.csv('diffSig.csv', header = T, row.names = 1)

GSE1 $Gene <- rownames(GSE1 )

GSE1 $GSE <- 'GSE1 '

colnames(GSE1 )

GSE1  <- subset(GSE1  select = c("Gene",    "GSE", "logFC") )

head(GSE1 )

# --------------------------------------------------------

# --------------------------------------------------------


data <- rbind(dt1, dt2, dt3....)

table(data$GSE)


data1 <- subset(data, data$logFC > 0)

data2 <- subset(data, data$logFC < 0)


data1$logFC <-NULL
data2$logFC <-NULL

dt1 <- as.data.frame(table(data1$GSE))

dt1$sig <- 'Up-regulation'

dt2 <- as.data.frame(table(data2$GSE))

dt2$Freq <- -(dt2$Freq)

dt2$sig <- 'Down-regulation'

p <- ggplot() + 
  geom_bar(data = dt1, aes(x=Var1, y=Freq, fill=sig),stat = "identity",position = 'dodge') +
  geom_text(data = dt1, aes(x=Var1,  y=Freq, label=Freq, vjust=-0.25))+
 geom_bar(data = dt2, aes(x=Var1, y=Freq, fill=sig),stat = "identity",position = 'dodge') +
  geom_text(data = dt2, aes(x=Var1,  y=Freq, label=Freq, vjust=-0.25))+
  scale_fill_manual(values=c("#0072B5","#BC3C28"))+
  theme(axis.text.x = element_text(angle = 45, hjust = 0.5, vjust=0.5))+xlab('')+ylab('')

pdf(file = "deg_show.pdf",width=6,height=6)
p
dev.off()