2017-09-26 254 views
1

我试图用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

你能不能帮我调试我的代码? 谢谢

回答

1

将对象添加到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,这将帮助你在正确的可读的方式编写代码。祝你好运!