从 R 中的方差分析 table 中提取 p 值
Extracting p-value from anova table in R
**
前面的步骤是这样的:
a <- rnorm(100, mean=2, sd=3)
response <- a * 1.5 + rnorm(100, mean=0, sd=1)
model <- lm(response ~ a)
vartest <- anova(model)
我想将 p 值 提取到与 a 系数相关联的向量 中,该值 < 2.2e-16。
我的密码是:
vartest[1,5]
[1] 1.002182e-63
其中 vartest
产生以下方差 table。
我想知道我是否做错了,或者是否有替代方法可以将值提取到向量中?
我们可以使用 [[
或 $
直接提取列名
out <- vartest[["Pr(>F)"]][1]
is.vector(out)
#[1] TRUE
-检查OP的方法
identical(out, vartest[1,5])
#[1] TRUE
我们可以用str
检查对象的结构
str(vartest)
这将提供有关如何提取组件的想法
** 前面的步骤是这样的:
a <- rnorm(100, mean=2, sd=3)
response <- a * 1.5 + rnorm(100, mean=0, sd=1)
model <- lm(response ~ a)
vartest <- anova(model)
我想将 p 值 提取到与 a 系数相关联的向量 中,该值 < 2.2e-16。
我的密码是:
vartest[1,5]
[1] 1.002182e-63
其中 vartest
产生以下方差 table。
我想知道我是否做错了,或者是否有替代方法可以将值提取到向量中?
我们可以使用 [[
或 $
out <- vartest[["Pr(>F)"]][1]
is.vector(out)
#[1] TRUE
-检查OP的方法
identical(out, vartest[1,5])
#[1] TRUE
我们可以用str
str(vartest)
这将提供有关如何提取组件的想法