按天分组分割条形图

Question

我使用以下代码制作了以下条形图：

MD1<-read.csv("MD_qual_OTU_sorted.csv")

MD1<-data.frame(Samples=c("A","B","C","D","E","F","G","H","I","J","K","L","M", "N","O","P","Q", "R"), Number.of.OTUs=c(13,10,9,9,15,11,7,7,9,9,5,10,10,7,15,17,8,9))
par(las=1)

barplot(MD1[,2],names.arg=MD1[,1], ylab='OTU Count', yaxt='n', xlab='MD samples', main='Total OTU count/Sample',density=c(90,90, 90, 90, 90, 90, 10, 10, 10, 10, 10, 10, 40, 40, 40, 40, 40, 40), col=c("yellow","yellow","pink", "pink","green","green","red","red", "purple", "purple", "blue", "blue", "orange", "orange","cyan", "cyan","chartreuse4", "chartreuse4" ))

usr <- par("usr")
par(usr=c(usr[1:2], 0, 20))
axis(2, at=seq(0,20,5))

我想将样本A-F分成另一组（第3天），G-L（第5天）和M-R（第15天）

发布了类似的问题，但我不知道如何整理我输入数据的方式以便能够使用这些解决方案。

Answer 1

你可以考虑使用ggplot2，使用facet_wrap和facet_grid可以很容易地分开绘图。

library(ggplot2)


#create a grouping variable
MD1$Day <- rep(c("Day 03","Day 05","Day 15"),
               each=6)

p1 <- ggplot(MD1, aes(x=Samples,y=Number.of.OTUs)) +
  geom_bar(stat="identity") + facet_wrap(~Day,
                                         scales="free_x")
p1

或者，如果您想使用base-R并接近原始图像：     #add colors / densities

MD1$col <-  c("yellow","yellow","pink", "pink","green","green","red","red", 
              "purple", "purple", "blue", "blue", "orange", "orange","cyan", "cyan","chartreuse4", "chartreuse4" )
MD1$density <- c(90,90, 90, 90, 90, 90, 10, 10, 10, 10, 10, 10, 40, 40, 40, 40, 40, 40)

#set 1 row three cols for plotting
par(mfrow=c(1,3))
#split and plot
lapply(split(MD1, MD1$Day),function(x){
  barplot(x[,2],
          names.arg=x[,1],
          ylab='OTU Count',
          ylim=c(0,20),
          main=unique(x$Day),
          col=x$col,
          density=x$density)
})