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

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