求助ggplot画柱状图

R语言绘图

1. # log10函数把等比数列转化为等差数列
   
2. # log10(60) + test$plus/600 处理y轴截断，ggplot2本身不支持坐标轴截断
   
3. # 上面的函数实际是将0:600 的点拟合到 1.78:2.78上[也就是log10(60):log10(600)]
   
4. test$plus_0  <- ifelse(test$plus <600,log10(60) + test$plus/600, log10(test$plus)) -log10(60)
   
5. test$minus_0 <- ifelse(test$minus<600,log10(60) + test$minus/600,log10(test$minus))-log10(60)
   
6. 
   
7. # 这里假设log10(60)为0
   
8. breaks <- (log10(c(600,6000,60000,600000))-log10(60))
   
9. labels <- c("600","6000","60000","600000")
   
10. 
    
11. g.top  <- ggplot(test) + geom_bar(aes(x=id,y=plus_0),fill="red",stat="identity")+
    
12.   # 调整margin位置
    
13.   theme_bw() +theme(plot.margin = unit(c(38,5,-28,6),units="points")) +
    
14.   # 设置y轴值，标签
    
15.   scale_y_continuous(expand = c(0,0),breaks=breaks,labels=labels) +
    
16.   # 去除panel网格线
    
17.   theme(panel.grid.major.x = element_blank(),panel.grid.minor.x = element_blank(),panel.grid.minor.y = element_blank())+
    
18.   # 删除lable
    
19.   labs(y="")
    
20. g.down <- ggplot(test) + geom_bar(aes(x=id,y=minus_0),fill="blue",stat="identity")+
    
21.   # 调整margin位置，调整y轴lable位置
    
22.   theme_bw() +theme(plot.margin = unit(c(1,5,10,3),units="points"),axis.title.y = element_text(hjust =1.5)) +
    
23.   # 逆转y轴，设置y轴值，标签
    
24.   scale_y_reverse(expand = c(0,0),breaks=breaks,labels=labels) +
    
25.   # 去除panel网格线
    
26.   theme(panel.grid.major.x = element_blank(),panel.grid.minor.x = element_blank(),panel.grid.minor.y = element_blank())+
    
27.   # 设置lable
    
28.   labs(y="Reading Count",x="Chromosome") +
    
29.   # 设置lable大小，黑体 加粗
    
30.   theme(axis.title.y = element_text(face='bold',color='black',size=10)) +
    
31.   theme(axis.title.x = element_text(face='bold',color='black',size=10))
    
32. 
    
33. grid.arrange(g.top,g.down, heights = c(1/2, 1/2))

复制代码

谢谢，太感谢了