調整したp値を比較してggplotに追加

Question

調整したp値（bonferroniなど）をp値の代わりにggplot boxplotに追加すると助かります。私は、審美的に..p.adj..を使用して、ggpubパッケージのstat_compare_meansでそれをやろうとしましたが、比較リストを追加しても機能しません。これについて

ggplot(data= mydf, aes(x=B,y=A)) + 
    geom_boxplot() + 
    stat_compare_means(aes(label=..p.adj..), 
        comparisons = list(c("x","y"),c("x","z"),c("y","z"))) 

Answer 1

方法。

library(ggplot2) 
library(ggpubr) 
mydf <- data.frame(A=1:300, B=rep(c("x","y","z"),100)) 
my_comparisons<- list(c("x","y"),c("x","z"),c("y","z")) 

p<- ggplot(data= mydf, aes(x=B,y=A)) + geom_boxplot() + 
    theme_bw()+ 
    stat_compare_means(aes(label=..p.adj..), comparisons = my_comparisons, 
        label.x = 1.5, label.y = 300) 
# Add p-value 
p + stat_compare_means(label.y = 280, label.x = 1.2) 

# Change method 
p + stat_compare_means(method = "anova", label.y = 260, label.x = 2.2) 

私はまた、あなたがこのpost on R-bloggers

[1]: https://i.stack.imgur.com/6UGw4.png 
[2]: https://i.stack.imgur.com/oriH8.png