3

我正在制作一个带有 ggpairs{GGally) 的散点图矩阵，如下所示，但我想在上面的面板中显示我的 aov 结果中每个术语的 p 值，而不仅仅是出现的整体和按物种的相关值与包裹。

我怎样才能从这个 aov 结果中得到正确的列到我的上图中？我可以编写一个自定义函数来做到这一点，以及如何？甚至可以使用ggpairs吗？谢谢。

library(GGally);library(ggplot2)
pm <- ggpairs(data = iris, 
          mapping = aes(color = Species),
          columns = c("Sepal.Length", "Sepal.Width", "Petal.Length", "Petal.Width"))

pm

虹膜数据 ggpairs 图

result <- aov(Sepal.Length ~  Sepal.Width*Petal.Length, data = iris)
print(summary(result))

                      Df Sum Sq Mean Sq F value   Pr(>F)    
Sepal.Width                1   1.41    1.41    12.9 0.000447 ***
Petal.Length               1  84.43   84.43   771.4  < 2e-16 ***
Sepal.Width:Petal.Length   1   0.35    0.35     3.2 0.075712 .  
Residuals                146  15.98    0.11                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
4

0 回答 0