0

我在 R 中使用来自 gpubr 的 stat_compare_means() 绘图函数，并且试图为 wilcox.test 指定左尾（或替代 =“less”）规范。不幸的是，我无法获得让我能够这样做的功能。

任何人都可以协助指定替代方案吗？

我在 gpubr github 帐户上列出了一个问题，但尚未收到对此的反馈。

任何想法都会很棒。

MWE如下：

library(dplyr)
library(ggpubr)
library(ggplot2)

Mockdata <- bind_cols(data.frame(date = rep(seq(as.Date("2003/01/01"), as.Date("2003/12/01"), "months"), 6)) %>% arrange(date),
                  data.frame(Universe = rep(c("SPEURO", "SP500", "World"),24), Measure = rep(c("EBITDA", "GM"), 36), Value = rnorm(72, 5, 2))) %>% mutate_at(.vars= vars(Measure, Universe), .funs = funs(as.character))        

Median <- Mockdata %>% group_by(Universe, Measure) %>% summarise(Median.Measure = median(Value)) %>% filter(Universe == "JALSHAll")

# Boxplot:
ggboxplot(Mockdata, x = "Universe", y = "Value", color = "Universe", add = 
"jitter", legend = "none", facet.by = "Measure", )  +
 stat_compare_means(aes(label = paste0("KW: p = ", ..p.format..)))
4

1 回答 1

0

使用 ggsignif 包中的 geom_signif 查看以下可能的解决方案：

library(ggsignif) 

ggboxplot(Mockdata, x = "Universe", y = "Value", color = "Universe", add = "jitter", legend = "none", facet.by = "Measure", ) +                        
stat_compare_means(aes(label = paste0("KW: p = ", ..p.format..)), label.y.npc = "bottom") +  
geom_signif(test = "wilcox.test", test.args = list(alternative = "less", paired = TRUE), 
          comparisons = list(c("SPEURO", "SP500"), c("SPEURO", "World")),
          vjust = 0, step_increase = 0.12, map_signif_level = TRUE)

可以在参数中指定其他测试test.args参数。

于 2017-10-08T18:06:08.340 回答