geom_point 散點分布圖 配色


1.  GATK GenotypeConcordance

p<-ggplot(data=data1,mapping=aes(x = TRUTH_SAMPLE,y=NON_REF_GENOTYPE_CONCORDANCE,colour=VARIANT_TYPE,shape=VARIANT_TYPE ))+geom_point(size=2)+scale_color_brewer(palette = 'Accent')+labs(x="Samples",y="Concordance",title="Zbolt_15X vs Downloded_30X")+theme(  panel.grid=element_blank(),axis.text=(element_text(size=6)),axis.title=(element_text(size=11)),plot.title=element_text(hjust = 0.5))+scale_x_discrete(labels=paste0("",1:50))

p1<-p+geom_hline(yintercept =0.989,color="#44758E",linetype="dashed",size=1)+geom_hline(yintercept =0.804,color="#44758E",linetype="dashed",size=1)+annotate('text',x=25,y=1,label="SNP Mean Concordance:0.989",size=3)+annotate('text',x=25,y=0.815,label="INDEL Mean Concordance:0.804",size=3)

ggsave("./zbolt_test.png",width = 8,height = 6,plot=p1)

 

 

注:

因為x軸是樣本名,不是連續型變量,所以不能用scale_x_continuous(),得用scale_x_discrete(labels=paste0("",1:50)),生成字符串向量。這種樣本名不重要的情況,可以直接用1:length(sample)數值向量替換,便於操作。

geom_point(size=2) 點的大小

scale_color_brewer(palette = 'Accent') 設置點不同變量的配色模版

labs 設置x,y,標題的名稱

axis.text 軸刻度名的大小 就是1 2 3..50的大小,防止太密了重疊。

+scale_x_discrete(labels=paste0("",1:50)) 復寫橫坐標label 標簽。因為之前樣本名字數太多,沒必要

annotate('text',x=25,y=1,label="SNP Mean Concordance:0.989",size=3)+annotate('text',x=25,y=0.815,label="INDEL Mean Concordance:0.804",size=3) 這里x,y 表示文字中心 對應在圖坐標的位置。size是字體大小,可以參照點設置尺寸的大小

 

顏色取色器:在AI打開圖片,用吸管工具點擊要取的區域,在右邊屬性欄點擊填色,#nnmmdd 即與ggplot2通用的十六進制

      選擇顏色:

          1)可以在此AI此界面,想要顏色對應的十六進制,;

          2)根據數據類型,選擇配色建議ggplot配色建議 https://colorbrewer2.org/ 

          3)ggplot2自帶配色模版

 

通過element_text (屬於ggplot2主題系統)來設置圖中各文本元素的格式。如用於:

 

   axis.title.x x軸標簽/標題

  axis.title.y y軸標簽/標題

  axis.text.x x軸刻度標簽

  axis.text.y y軸刻度標簽

  legend.title 圖例標題文本

  legend.text 圖例文本

  plot.title 圖的主標題

  plot.subtitle 圖副標題

  plot.caption 圖的腳注

  plot.tag 繪圖的標簽

 # ggplot 允許多個調控函數共同作用,即圖層疊加。如geom_text ,annotate,element_text 共同作用很正常


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM