爾云間?? 一個專門做科研的團隊?? 關注我們

原創(chuàng) ? 小果 ? 生信果

各位小伙伴大家好，小果又和大家見面了，小果最近發(fā)現(xiàn)通過繪制百分比柱狀圖我們能夠直觀的看到一個類別在另一類別不同分組的分布情況

如下圖

這個代碼也很簡單，小果給大家演示一下

rm(list=ls());gc();

setwd("D:/rstudy/")#建立工作目錄

library(RColorBrewer)

library(scales) # 將y軸轉化成百分比格式

library(tidyverse)

library(ggplot2)

mydata <- read.csv("sur_m.csv",header = T)

table(mydata$Molecular.subtype, mydata$riskscore)

library(scales)

ggplot(mydata, aes(x = Molecular.subtype, fill = riskscore)) +

geom_bar(width = 0.3, position = "fill") + # 百分比柱狀圖

scale_fill_brewer(palette = "Blues") + # 調色板{RColorBrewer}

scale_y_continuous(labels = percent) + # 修改y軸刻度為百分比

guides(fill=guide_legend(title = "")) +

labs(title = "",

x = "Molecular.subtype",

y = "riskscore") +

theme_minimal() +

coord_flip() # 倒轉x與y軸

這樣就結束了，是不是很簡單

好了，以上就是小果的分享，是不是干貨滿滿呢，小伙伴快去試試吧。

推薦閱讀

• 多種花樣都能懂，富集圖的看法

• 相關性圖譜之相關性熱圖學習

• 小果教你快速優(yōu)雅使用Genecards數(shù)據(jù)庫

• 從UCSC下載TCGA數(shù)據(jù)

• 簡單高效利用Batch Entrez批量獲取基因別名

關注小果，小果將會持續(xù)為你帶來更多生信干貨哦。

“生信果”，生信入門、R語言、生信圖解讀與繪制、軟件操作、代碼復現(xiàn)、生信硬核知識技能、服務器、生物信息學的教程，以及基于R的分析和可視化等原創(chuàng)內容，一起見證小白和大佬的成長。