如何在R中使用ggbetweenstats处理不同时间点的分组数据？

我在R 4.3.2版本中使用ggstatsplot包的ggbetweenstats函数，想要绘制一张展示2个不同处理组在3个不同时间点的数据分布情况的图表。

目前使用grouped_ggbetweenstats函数会生成基于处理组的facet_wrap分面效果，得到两个独立的图表。我可以通过ggplot2的geom_boxplot结合ggsignif包添加统计信息，但更喜欢ggstatsplot生成的图表样式。

示例数据 #

Treatment <- sample(c("A", "B"), 100, replace=T) 
Timepoint <- sample(c("1", "2", "3"), 100, replace=T) 
Value <- sample(20:80, size=100, replace=T) 
df <- data.frame(Treatment, Timepoint, Value) 
library(ggstatsplot)

当前尝试的代码（分面效果） #

这段代码会生成两个独立的分面图表：

grouped_ggbetweenstats(data=df, x=Timepoint, y=Value,
                       grouping.var=Treatment, results.subtitle=F) 

期望的图表样式 #

我想要的是同一面板内按处理组分色叠加的样式，类似以下ggplot2代码生成的效果：

library(ggplot2)

ggplot() + 
  geom_boxplot(data=df, aes(x=Timepoint, y=Value, color=Treatment)) 

请问是否可以用ggstatsplot包直接生成这种样式的图表，还是需要自行组合geom_boxplot、geom_point、geom_violin、geom_signif等函数？

内容的提问来源于stack exchange，提问作者Spencer Epp