话不多说直接记录library(ggplot2) library(dplyr) library(tidyr) ##包含三列信息基因样本扩增数 df - read.delim(cnv_matrix.txt) plot_data - drop_na(df) p - ggplot(plot_data, aes(x Tumor_Sample_Barcode, y CNV, color Hugo_Symbol, group Hugo_Symbol)) geom_hline(yintercept 3, linetype dashed, linewidth 0.7, colour black) geom_line(linewidth 0.6) geom_point(size 3.5) scale_color_manual(values c(ERBB2 #c41e3a, TOP2A #286090)) ylim(2, 8) labs( x Patient samples, y Copy number (CN), color Symbol ) theme_bw() theme( axis.text.x element_text(angle 45, hjust 1), panel.grid element_blank(), legend.position right ) ggsave(filename cnv_plot.png, plotp, width 5, height 4) ggsave(filename cnv_plot.pdf, plotp, width 5, height 4)出图结果为
话不多说直接记录library(ggplot2) library(dplyr) library(tidyr) ##包含三列信息基因样本扩增数 df - read.delim(cnv_matrix.txt) plot_data - drop_na(df) p - ggplot(plot_data, aes(x Tumor_Sample_Barcode, y CNV, color Hugo_Symbol, group Hugo_Symbol)) geom_hline(yintercept 3, linetype dashed, linewidth 0.7, colour black) geom_line(linewidth 0.6) geom_point(size 3.5) scale_color_manual(values c(ERBB2 #c41e3a, TOP2A #286090)) ylim(2, 8) labs( x Patient samples, y Copy number (CN), color Symbol ) theme_bw() theme( axis.text.x element_text(angle 45, hjust 1), panel.grid element_blank(), legend.position right ) ggsave(filename cnv_plot.png, plotp, width 5, height 4) ggsave(filename cnv_plot.pdf, plotp, width 5, height 4)出图结果为
)
)
