ggplot2的水平條形圖

我試圖用ggplot2從這個數據框中做條形圖來比較Uniprot和macki列的數據。ggplot2的水平條形圖

Famille Uniprot Macki 
SN_02  26   7 
SN_03  21   22 
SN_04  16   7 
SN_05  4   0 
SN_09  10   0 
SN_10  5   0 
SN_17  6   0 
SN_19  13   4 
SN_20  3   1 
SN_31  2   0 
SN_32  5   3 
SN_33  3   0 
SN_34  1   0 
SN_37  3   0

這是我的R代碼裏面：

bar <- c(26, 21, 16,4,10,5,6,13,3,2,5,3,1,3,7,22,7,0,0,0,0,4,1,0,3,0,0,0) 
lab <- c("SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17","SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37", "SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37") 
pct <- round(bar/sum(bar)*100) 
lab <- paste(lab,"%", sep = "") 
lab <- paste(lab, pct) 
lab0 <- c("SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37", "SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37") 
nlab <- length(lab) 
type <- rep("Uniprot:MACKI", each=nlab/2) 
library(ggplot2) 
theme_set(theme_bw()) 
dd <- data.frame(lab0, type, bar, pct) 
dd$lab0 <- reorder(dd$lab0,-dd$bar) 
ggplot(dd,aes(x=lab0,y=bar,fill=lab0))+geom_bar(aes(alpha=factor(bar)), stat = "identity", 
position=position_dodge(width=1))+scale_alpha_discrete(range=c(0.5,1.0))+ 
geom_text(aes(label=paste0(pct,"%"),group=interaction(lab0,type)),hjust=-0.5, position=position_dodge(width=1)) 
+ coord_flip()+expand_limits(y=20)+labs(x="",y="total")

這是情節我已經得到了this horizontal bar plot

你能不能幫我調試我的代碼？謝謝

來源

2017-09-26 Anicet Ebou

將對象添加到ggplot2是通過+，它必須在同一行。在最後一行看一下 - 你必須把+放在那裏。修改代碼工作

bar <- c(26, 21, 16,4,10,5,6,13,3,2,5,3,1,3,7,22,7,0,0,0,0,4,1,0,3,0,0,0) 
lab <- c("SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17","SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37", "SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37") 
pct <- round(bar/sum(bar)*100) 
lab <- paste(lab,"%", sep = "") 
lab <- paste(lab, pct) 
lab0 <- c("SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37", "SN_02", "SN_03", "SN_04", "SN_05", "SN_09", "SN_10", "SN_17", "SN_19", "SN_20", "SN_31", "SN_32", "SN_33", "SN_34", "SN_37") 
nlab <- length(lab) 
type <- rep("Uniprot:MACKI", each=nlab/2) 
library(ggplot2) 
theme_set(theme_bw()) 
dd <- data.frame(lab0, type, bar, pct) 
dd$lab0 <- reorder(dd$lab0,-dd$bar) 
ggplot(dd,aes(x=lab0,y=bar,fill=lab0))+ 
    geom_bar(aes(alpha=factor(bar)) , stat = "identity", position = position_dodge(width=1)) + 
    scale_alpha_discrete(range=c(0.5,1.0))+ 
    geom_text( 
    aes(label=paste0(pct,"%") , group=interaction(lab0,type)), 
    hjust = -0.5, 
    position = position_dodge(width=1)) + 
    coord_flip()+ 
    expand_limits(y = 20) + 
    labs(x = "",y = "total")

後，我勸你也跟着r-code-formatting-guide，這將幫助你在正確的可讀的方式編寫代碼。祝你好運！

來源

2017-09-26 13:00:10 Gonzo

ggplot2的水平條形圖

回答

相關問題