
有了R語言的基礎,以及ggplot2繪圖基礎,我們的「生信常用分析圖形」的繪製就可以提上日程了!本系列,師兄就開始帶着大家一起學習如何用R語言繪製我們自己的各種分析圖吧!
❝由於本系列的所有分析代碼均為師兄細心整理和詳細注釋而成的,所以本系列中「部分教程需要收費」的哦!付費即可獲取對應小節的繪圖數據和代碼!
此外,如果有小夥伴想要「一次性購買本系列所有代碼和示例數據的」,也可以直接公眾號右下角添加客服微信諮詢 「本系列所有代碼打包僅售:99元」,並且還包括「跟着高分SCI學作圖系列」的「全部示例數據、代碼和付費文章」哦!
最後,已經購買的小夥伴還可以加入師兄的 「R語言生信繪圖」交流群,群內可以討論繪圖上的疑問,我也會在閒暇的時候給小夥伴們答疑!
❞系列內容包括:
高級函數繪製火山圖--ggplot2、ggpurb歷史圖形展示
本期瀑布圖展示
image-20220413115418574圖片意義解讀還是那句話,作圖的前提是理解圖形要表達什麼含義,否則做了也白做!瀑布圖之所以看起來像瀑布,和它的計算規則有關:
「右側柱狀圖:」找出所有樣本中突變個數最多的基因,倒序排列;右側柱狀圖的呈現的是這些基因在所有樣本中的突變次數,次數高的排在上面。
「上方柱狀圖:」將包含這個基因突變的樣本集中排列,且根據樣本中包含的突變個數進行倒序排列;上方柱狀圖呈現的是每個樣本中含有的突變個數,並且不同樣本按照相同突變基因集中排列。
「柱狀圖中不同的顏色」:代表不同突變的鹼基變化類型;(如:C -> T,T -> A等)
「瀑布圖中的不同顏色」:代表不同突變的意義類型;(如:Missense mutation -- 錯義突變,Nonsense mutation -- 無義突變等)
基礎繪圖#安裝並加載所需的R包#BiocManager::install("maftools")library(maftools)#查看示例數據#maf格式的基因突變信息laml.maf<-system.file('extdata','tcga_laml.maf.gz',package='maftools')#臨床表型注釋信息laml.clin<-system.file('extdata','tcga_laml_annot.tsv',package='maftools')#使用read.maf函數讀取數據laml<-read.maf(maf=laml.maf,clinicalData=laml.clin)#查看maf對象laml#獲取maf對象匯總信息getSampleSummary(laml)getFields(laml)#使用plotmafSummary函數可視化maf對象匯總信息plotmafSummary(maf=laml,rmOutlier=TRUE,addStat='median',dashboard=TRUE,titvRaw=FALSE)
image-20220413114318314oncoplot繪製瀑布圖基礎瀑布圖###################使用oncoplot函數繪製基因突變瀑布圖###########3#展示top10變異基因的信息pdf("onco_plot.pdf",height=6,width=7)oncoplot(maf=laml,top=10)dev.off()
image-20220413114651178自定義變異類型的顏色#自定義變異類型的顏色library(RColorBrewer)vc_cols<-brewer.pal(8,"Set1")names(vc_cols)<-levels(laml@data$Variant_Classification)head(vc_cols)##Frame_Shift_DelFrame_Shift_InsIn_Frame_DelIn_Frame_Ins##"#E41A1C""#377EB8""#4DAF4A""#984EA3"##Missense_MutationNonsense_Mutation##"#FF7F00""#FFFF33"pdf("onco_plot1.pdf",height=6,width=7)oncoplot(maf=laml,top=20,colors=vc_cols)dev.off()
image-20220413115249180添加臨床注釋信息#添加臨床注釋信息,按注釋類型進行排序names(laml@clinical.data)##[1]"Tumor_Sample_Barcode""FAB_classification"##[3]"days_to_last_followup""Overall_Survival_Status"oncoplot(maf=laml,top=20,colors=vc_cols,clinicalFeatures="FAB_classification",sortByAnnotation=T)
image-20220413115439983#展示多個臨床注釋信息oncoplot(maf=laml,top=20,colors=vc_cols,clinicalFeatures=c("FAB_classification","Overall_Survival_Status"),sortByAnnotation=T)
image-20220413115418574示例數據和代碼獲取❝本系列「所有代碼和示例數據將會和生信常用圖形系列繪圖」放在一起,公眾號右下角添加師兄微信,「付費99元,即可加入生信繪圖交流群」。群內不僅提供生信常用圖形系列的代碼,還會「提供本系列後續所有Figure的示例數據和代碼」,我會在文章更新後第一時間上傳。
當然了!如果你還想白嫖,師兄的文章中代碼已經寫的很清楚了!但是師兄還是希望你「點個讚再走」唄!
以上就是本期的全部內容啦!「歡迎點讚,點在看」!師兄會儘快更新哦!製作不易,你的打賞將成為師兄繼續更新的十足動力!
❞往期文章跟着Nature Medicine學作圖 -- 箱線+散點圖跟着Nature Communications學作圖 -- 漸變火山圖跟着Nature Communications學作圖 -- 氣泡圖+相關性熱圖跟着Nature Communications學作圖 -- 複雜提琴圖跟着Nature Medicine學作圖 -- 複雜熱圖生信常用分析圖形的繪製01 -- 各種類型的熱圖你學會了嗎?生信常用分析圖形的繪製02 -- 解鎖火山圖真諦!生信常用分析圖形繪製06 -- 富集分析分組柱狀圖+氣泡圖生信常用分析圖形繪製14 -- mantel test相關性圖