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