条形图以及分组条形图

生信喵

一、条形图

1. # 实践：绘制人染色体长度分布图
   
2. x <- read.csv(file = "homo_length.csv",header = T)
   
3. head(x)
   
4. x <- x[1:24,]
   
5. 
   
6. barplot(height = x$length)  
   
7. barplot(height = x$length,names.arg = x$chr)  
   
8. barplot(height = x$length,names.arg = x$chr,las =3)  
   
9. barplot(height = x$length,names.arg = x$chr,las =3,col = rainbow(4))  
   
10. barplot(height = x$length,names.arg = x$chr,las =3,col = 'red')  
    
11. barplot(height = x$length,names.arg = x$chr,las =3,col = c('red','green'))
    
12. library(RColorBrewer)  
    
13. display.brewer.all()
    
14. cols <- brewer.pal(n = 6,name = "RdPu")  
    
15. barplot(height = x$length,names.arg = x$chr,col = cols)
    
16. cols <- brewer.pal(n = 6,name = "Dark2")  
    
17. barplot(height = x$length,names.arg = x$chr,col = cols)
    
18. cols <- brewer.pal(n = 6,name = "Blues")  
    
19. barplot(height = x$length,names.arg = x$chr,col = cols)
    
20. cols <- brewer.pal(n = 6,name = "Set1")
    
21. barplot(height = x$length,names.arg = x$chr,col = cols,las=3)
    
22. barplot(height = x$length,names.arg = x$chr,col = cols,horiz = T,las=2)
    
23. par('mar')
    
24. barplot(height = x$length,names.arg = x$chr,col = cols,horiz = T,las=2,
    
25.         border = F,width = c(1,2),space = 0.1)
    
26. barplot(height = x$length,names.arg = x$chr,las=2,
    
27.         border = F,width = c(1,2),space = 1,density = 12,angle = c(45,135))
    

复制代码

人染色体长度分布图


二、分组条形图

1. #绘制分组条形图
   
2. x <- read.csv("sv_distrubution.csv",header = T,row.names = 1)  
   
3. head(x)
   
4. #barplot(x)  
   
5. barplot(as.matrix(x))  
   
6. barplot(t(as.matrix(x)))  
   
7. # barplot(t(as.matrix(x)),col = rainbow(4))  
   
8. barplot(t(as.matrix(x)),col = brewer.pal(4,name = 'Set2'))
   
9. barplot(t(as.matrix(x)),col = brewer.pal(4,name = 'Set2'),beside = T)  
   
10. barplot(t(as.matrix(x)),col = brewer.pal(4,name = 'Set2'),
    
11.         horiz = F,density = 30,angle = c(45,135))
    
12. barplot(t(as.matrix(x)),col = brewer.pal(4,name = 'Set2'))
    
13. barplot(t(as.matrix(x)),col = brewer.pal(4,name = 'Set2'),legend.text = colnames(x),ylim = c(0,800),  
    
14.         main = "SV Distribution",xlab="Chromosome Number",ylab="SV Numbers")  
    

复制代码

堆叠条形图