如何从Manova输出中提取p值

Question

我正在通过R中的manova公式测试Hotelling ^2测试。我正在测试不同的相同大小，因此有多个Manova表输出。下面是我如何为整个样本生成manova的

attach(iris)
library(Hotelling)
library(corpcor)
s= iris[1:100,1:5]
input= cbind(s$Sepal.Length,s$Sepal.Width, s$Petal.Length, s$Petal.Width )
m= manova(input~ Species, data = s)
summary(m, "Hotelling-Lawley")

我想知道如何从每个表中提取p值。我试着追随，但没有这样的运气：

res$"Pr(>F)"

res$p.value

summary(man)[8]

但是每个返回NULL

Answer 1

在您的示例中，p非常小：

summary(m, "Hotelling-Lawley")$stats

          Df Hotelling-Lawley approx F num Df den Df       Pr(>F)
Species    1         26.33509 625.4583      4     95 2.664857e-67
Residuals 98               NA       NA     NA     NA           NA

对于给定的预测器，例如Species，可以隔离p值，如下所示：

summary(m, "Hotelling-Lawley")$stats["Species", "Pr(>F)"]
[1] 2.664857e-67

Docs here。

我知道这只是iris的一个测试用例，但即使如此:考虑到当p值如此小时，它作为有效的测试统计数据开始失去意义。相反，您可以选择效果大小的度量，甚至是描述性统计，以支持您的结果。