从R中的两个线性模型的比较中提取anova()的p值

Question

我很抱歉没有提供足够的信息来正确访问这个问题.无论如何,我想比较不同的线性模型对并将它们的p值记录到表中.我的麻烦在于我无法仅提取p值.

TC <- matrix(data=NA, nrow=1, ncol=6)
    ML5 <- 1:5
for(B in ML5) {
    Coop <- anova(M6,Models5[[B]])$"Pr(>F)"
    TC[1,B] <- Coop
}

 >Coop
   NA 2.041767e-05

那么我如何才将数字放入表中而不是将NA放入表中？

太感谢了!

Answer 1

你只需要引号:

pvalue <- anova(m2,m1)$"Pr(>F)"

您可以使用普通括号子集访问pvalue的第二个元素:

pvalue[2]

所以在你的例子中我相信你会使用

Coop <- anova(M6,Models5[[B]])$"Pr(>F)"[2]

(虽然无法访问M6,但我不能完全确定).

str在这种情况下,该命令非常有用,可以找出您正在处理的对象类型:

str(myanova$"Pr(>F)")