细菌丰度堆积条形图的制作

Question

我有152个样本(行)的相对细菌丰度的数据框架。我想画一个堆叠的条形图，显示所有样本中每个细菌组的总体丰度(例如，放线菌vs.杆菌vs. Firmicutes等)。我想要它的颜色编码，以及这方面的图例。有没有人能建议一下怎么做？我的问题是，我不确定如何在R中获得用于绘图的列总数。谢谢。

row.names       Actinobacteria  Bacteroidetes   Firmicutes  Fusobacteria    Proteobacteria  Verrucomicrobia Other
1   sample1 0.0084246282    0.41627099  0.55475503  0.000000e+00    7.245180e-04    5.391762e-05    1.977092e-02
2   sample2 0.0168571327    0.13298800  0.80289437  3.560112e-05    4.272135e-03    4.238314e-02    5.696180e-04
3   sample3 0.0020299288    0.53813817  0.42367947  3.311006e-02    7.978327e-04    3.534702e-05    2.209189e-03

Answer 1

我不清楚样本名称是否是数据帧中的行名，所以我简单地重新创建了数据框，将样本名称放在一个变量中，与细菌名称相同：

Sample Actinobacteria Bacteroidetes Firmicutes Fusobacteria Proteobacteria
1 sample1    0.008424628     0.4162710  0.5547550 0.000000e+00   0.0007245180
2 sample2    0.016857133     0.1329880  0.8028944 3.560112e-05   0.0042721350
3 sample3    0.002029929     0.5381382  0.4236795 3.311006e-02   0.0007978327
  Verrucomicrobia       Other
1    5.391762e-05 0.019770920
2    4.238314e-02 0.000569618
3    3.534702e-05 0.002209189

要重现此数据集，可以运行以下命令：

df <- structure(list(Sample = structure(1:3, .Label = c("sample1", 
"sample2", "sample3"), class = "factor"), Actinobacteria = c(0.0084246282, 
0.0168571327, 0.0020299288), Bacteroidetes = c(0.41627099, 0.132988, 
0.53813817), Firmicutes = c(0.55475503, 0.80289437, 0.42367947
), Fusobacteria = c(0, 3.560112e-05, 0.03311006), Proteobacteria = c(0.000724518, 
0.004272135, 0.0007978327), Verrucomicrobia = c(5.391762e-05, 
0.04238314, 3.534702e-05), Other = c(0.01977092, 0.000569618, 
0.002209189)), .Names = c("Sample", "Actinobacteria", "Bacteroidetes", 
"Firmicutes", "Fusobacteria", "Proteobacteria", "Verrucomicrobia", 
"Other"), class = "data.frame", row.names = c("1", "2", "3"))

正如@zx8754所建议的，此数据帧需要重塑，即从宽格式移动到长格式。有关更多信息，请查看此link中的几个示例。

如果上面的数据帧被命名为df，则以下命令将以长格式对其进行整形：

library(reshape2)
df_long <- melt(df, id.vars = "Sample", variable.name = "Phyla")

在这里，我们可以使用ggplot进行绘图：

library(ggplot2)
ggplot(df_long, aes(x = Sample, y = value, fill = Phyla)) + 
    geom_bar(stat = "identity")

这就给出了：

​

​