使用ggpubr从双向anova中添加pvalues到boxplot

Question

我有一个包含多个因素的数据集，目前我正在将两个因素与箱形图进行比较。

我已下载ggpubr，正在使用stat_compare_means()添加p值。

但是，我希望所有4次比较都有p值。这只给了我2。

我查看了在线示例，我可以看到如何使用my_comparisons <- list( c("0.5", "1"), c("1", "2"), c("0.5", "2") )进行自己的比较 但这是针对同一因素的不同层次。

因此，例如在牙齿生长示例中，我希望p值将OJ/D1与OJ/D0.5进行比较。

这是我的代码：

OTU<-colnames(df)[j]
ggplot(subset(df,variable==OTU),
   aes(x=Fungi,y=abund,fill=Insect))+geom_boxplot()+
stat_compare_means()+
 ggtitle(paste(OTU,location,sep="_”))

这是我的数据：

> df[1:5,1:5]
      Location Bean Fungi Insect Acaulospora
R-B1       Root Bean    NF     NI           0
R-B2       Root Bean    NF     NI           0
R-B3       Root Bean    NF     NI           0
R-BF-1     Root Bean Fungi     NI           0
R-BF-2     Root Bean Fungi     NI           0

我想要Fungi_Insect vs NF_Insect以及Fungi_Insect与Fungi_NI等的p值。

任何帮助表示感谢。