r - 如何从 Manova 输出中提取 p 值

Question

我正在通过 R 中的 manova 公式测试 Hotelling T^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

可以为给定的预测变量分离 p 值，例如Species，像这样：

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

文档在这里。

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