垂直堆叠核密度图

Question

我有一个来自贝叶斯模型的输出矩阵 (C)，其中包含 3000 行，其中包含给定鸟类繁殖行为（列；唱歌、孵化、雏鸟等）的周数 (1-13)发生。我已经使用以下代码可视化了最有可能发生行为的那一周的内核密度估计：

G <- mcmc_dens(C, pars = c("Singing", "Building", "Incubating", "Nestlings", "Empty Nest", "Fledglings Observed", "Fledgling/Adult Interactions", "Fledgling Foraging"))
G <- G + theme(axis.title = element_text(face="plain",size=12)) + labs(x ="Week") + scale_x_continuous(breaks = 1:13)

...产生这些数字：

我想将这些数字堆叠在一起，以便我有一个具有相同 x 轴的数字，您可以在其中轻松查看哪些行为同时达到峰值，但我不知道如何使用 { {1}}（即我希望唱歌的图表高于建筑物，唱歌和建筑物都高于孵化，依此类推，以便我有八个垂直对齐的图表）。

来自矩阵 C 的数据样本（不包括所有列）：

mcmc_dens

Answer 1

想通了！ mcmc_dens 有参数 facet_args 将每个图形变成它自己的面（我花了这么长时间，因为我不熟悉面）。修改我的原始代码的第一行给了我我要找的图：

pars <- c("Singing", "Building", "Incubating", "Nestlings", 
          "Empty", "Fledglings", "Interactions", "Foraging")
G <- mcmc_dens(C, pars=pars, facet_args=list(ncol=1, strip.position="left"))

这是现在的图像：