R语言绘制表达矩阵的热图与火山图

火山图

仅需要pval 和log2FoldChange

v_circ$color <- ifelse(v_circ$pval <0.05 & abs(v_circ$log2FoldChange)>= 1,ifelse(v_circ$log2FoldChange > 1,'red','blue'),'gray')
color <- c(red = "red",gray = "gray",blue = "blue")

circ_vo <- ggplot(v_circ, aes(log2FoldChange, -log10(pval), col = color)) +
  geom_point() +
  theme_bw() +
  scale_color_manual(values = color) +
  labs(x="log2 (fold change)",y="-log10 (p-value)") +
  geom_hline(yintercept = -log10(0.05), lty=4,col="grey",lwd=0.6) +
  geom_vline(xintercept = c(-1, 1), lty=4,col="grey",lwd=0.6) +
  theme(legend.position = "none",
        panel.grid=element_blank(),
        axis.title = element_text(size = 16),
        axis.text = element_text(size = 14))

circ_vo + ggtitle('circRNA')
circ_vo

热图

circ_heat <- pheatmap(log2(h_circ+1), legend = FALSE, cluster_cols = FALSE, scale = 'row',clustering_method = "average", 
                      show_rownames = FALSE, color = colorRampPalette(c("navy", "white", "firebrick3"))(50), main = 'circRNA')
circ_heat

需要用log（x+1) 标准化表达量