2

我有 9 个不同样本中各种门的相对丰度（以 % 为单位）的分子测序数据，我试图将其绘制为颜色编码的条形图（其中每个门对应于不同的颜色）。在 excel 上很简单，但是对于 R 上的一个完整的新手来说，我很挣扎。我的数据采用 excel 格式（格式为制表符），其中第一行是标签（例如样本名称）- 绘制时，条形标签错位，不匹配，R 绘制我的 excel 文件的第一行（名称）作为单独的值（附图片）。到目前为止，我所拥有的是：

attach(data) 
data.1<-as.matrix(data)
par(mfrow=c(1,1))
barplot(data.1, col=c("aquamarine3","azure2","blue2","brown3","cadetblue3","deepskyblue3","firebrick3","gold3","darkorange3","darkorchid3","darkseagreen","darkslateblue","darkviolet","deeppink4"), main=".", xlab="Unit/Treatment", ylab="% Relative abundance")
detach(data)

legend("topright", inset=c(-0.2,0),
       legend = c("Unassigned", "Acidobacteria","Actinobacteria","Bacteroidetes","Chlorobi","Chloroflexi","Firmicutes","Gemmatimonadetes","Planctomycetes","Proteobacteria","Verrucromicrobia","Euryarchaeota","Crenarchaeota","Parvarchaeota"), 
       fill = c("aquamarine3","azure2","blue2","brown3","cadetblue3","deepskyblue3","firebrick3","gold3","darkorange3","darkorchid3","darkseagreen","darkslateblue","darkviolet","deeppink4"))
par(mar=c(5.1, 4.1, 4.1, 8.1), xpd=TRUE)
layout(mat, widths = rep.int(1, ncol(mat)),
       heights = rep.int(1, nrow(mat)), respect = FALSE)

结果，我得到了这个：

条形图尝试，其中 R 将我的样本名称绘制为 x_1 并因此移动其他标签。此外，我的图例涵盖了我的大部分条形图，我似乎无法调整它。

提前非常感谢 - 任何有助于使条形图看起来体面的帮助将不胜感激。

4

0 回答 0