按组更改 siberDensityPlot 颜色

Question

我正在使用 SIBER 包来开发包含三组的绘图，我想根据图例为每个组着色。 SIBER pdf 说要使用 clr 命令，其中指出

“用于对每个框区域进行着色的颜色矩阵。默认为灰度 grDevices::gray((9:1)/10)，复制与 dat 中一样多的列。当用户指定时，行包含 probs 中指定的每个置信区域的颜色，columns 代表 dat 中数据列的 ecah。这样，每个组都可以有蓝色、红色和黄色的阴影。”

但是，如果您在帮助网站（https://rdrr.io/cran/SIBER/man/siberDensityPlot.html）上运行代码以获得“彩色”图而不是灰度图，其中矩阵为指定颜色：

# A basic default greyscale density plot
Y <- matrix(stats::rnorm(1000), 250, 4)
siberDensityPlot(Y)

# A more colourful example
my_clrs <- matrix(c("lightblue", "blue", "darkblue",
"red1", "red3", "red4",
"yellow1", "yellow3", "yellow4",
"turquoise", "turquoise3", "turquoise4"), nrow = 3, ncol = 4)
siberDensityPlot(Y, clr = my_clrs)

结果不是一个用三种蓝色阴影代表每个 CI 的箱线图，然后是另一个用红色阴影、黄色阴影等绘制的箱线图……您为每个 CI 创建一个用 3 种蓝色阴影表示的图，并在所有不同的“组”中重复或面板中的箱线图。

我的问题是：有谁知道如何制作图形，以便在箱线图中获得代表不同 CI 的颜色阴影，然后为另一组的 CI 获得不同颜色的阴影？附图有望进一步解释我想要实现的目标。

Answer 1

看起来像是某种错误。我看到您发布的链接生成了全蓝色图，但是当我在电脑上运行代码时，我得到了这个图。我正在运行 Windows 10 R 版本 4.0.2 以及来自 CRAN 的最新 SIBER 软件包。该链接似乎正在运行 linux 和 R 3.4.4。

Answer 2

很抱歉进入这个主题，但也是相关的.. 函数 siberDensityPlot 确实创建了我正在寻找的图，但它本身计算置信区间。这是我的问题。我可以重新创建这个图，但使用我自己的 50%、75% 和 95% 间隔值以及平均值吗？也许还有其他功能..??