close

有了R語言的基礎,以及ggplot2繪圖基礎,我們的「生信常用分析圖形」的繪製就可以提上日程了!本系列,師兄就開始帶着大家一起學習如何用R語言繪製我們自己的各種分析圖吧!


由於本系列的所有分析代碼均為師兄細心整理和詳細注釋而成的,所以本系列中「部分教程需要收費」的哦!付費即可獲取對應小節的繪圖數據和代碼!

此外,如果有小夥伴想要「一次性購買本系列所有代碼和示例數據的」,也可以直接公眾號右下角添加客服微信諮詢 「本系列所有代碼打包僅售:99元」,並且還包括「跟着高分SCI學作圖系列」的「全部示例數據、代碼和付費文章」哦!

最後,已經購買的小夥伴還可以加入師兄的 「R語言生信繪圖」交流群,群內可以討論繪圖上的疑問,我也會在閒暇的時候給小夥伴們答疑!

系列內容包括:

各種類型的熱圖你學會了嗎?
普通熱圖
環形熱圖
解鎖火山圖真諦!
plot函數就能畫火山圖?
高級函數繪製火山圖--ggplot2、ggpurb
經典富集分析及氣泡圖、柱狀圖繪製
氣泡圖繪製
柱狀圖繪製
富集分析圈圖
富集分析弦圖
繪製一張可以打動審稿人的桑基圖
生存分析 -- KM曲線圖
基礎PCA圖
雲雨圖
韋恩圖
環形互作網絡圖
相互作用網絡圖
聚類樹美化
富集分析氣泡圖進階版
mantel test相關性圖
詞雲圖
瀑布圖
等等,想到再繼續補充!!!
歷史圖形展示

本期瀑布圖展示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學作圖 -- 複雜熱圖
跟着高分SCI學作圖 -- 氣泡火山圖
跟着高分SCI學作圖--複雜熱圖02
生信常用分析圖形的繪製01 -- 各種類型的熱圖你學會了嗎?
生信常用分析圖形的繪製02 -- 解鎖火山圖真諦!
生信常用分析圖形的繪製03 -- 富集分析圈圖
生信常用分析圖形的繪製04 -- 桑基圖
生信常用分析圖形的繪製05 -- 弦圖
生信常用分析圖形繪製06 -- 富集分析分組柱狀圖+氣泡圖
生信常用分析圖形繪製07 -- PCA圖
生信常用分析圖形繪製08--雲雨圖
生信常用分析圖形繪製09 -- Venn圖
生信常用分析圖形繪製10 -- 環形互作網絡圖
生信常用分析圖形繪製11 -- 相互作用網絡圖
生信常用分析圖形繪製12 -- 聚類樹美化
生信常用分析圖形繪圖13 -- 富集進階氣泡圖
生信常用分析圖形繪製14 -- mantel test相關性圖
生信常用分析圖形繪製15 -- 詞雲圖
AI繪圖、R語言基礎及ggplot2繪圖教程合集
arrow
arrow
    全站熱搜
    創作者介紹
    創作者 鑽石舞台 的頭像
    鑽石舞台

    鑽石舞台

    鑽石舞台 發表在 痞客邦 留言(0) 人氣()