R 带有交互项的联合显著性检验:包含错误的系数/变量名称
我试图在R中进行联合显著性测试:R 带有交互项的联合显著性检验:包含错误的系数/变量名称,r,linear-regression,significance,R,Linear Regression,Significance,我试图在R中进行联合显著性测试: library(car) data("mtcars") mylm <- lm(mpg ~ qsec + gear + am + am:qsec + am:hp, data=mtcars) linearHypothesis(mylm, c("am + am:qsec + am:hp")) 我想测试的是 am+am:qsec+am:hp=0 我在文档中找到了如何测试所有交互术语: linearHypothesis(mylm, matchCoefs(mylm
library(car)
data("mtcars")
mylm <- lm(mpg ~ qsec + gear + am + am:qsec + am:hp, data=mtcars)
linearHypothesis(mylm, c("am + am:qsec + am:hp"))
我想测试的是
am+am:qsec+am:hp=0
我在文档中找到了如何测试所有交互术语:
linearHypothesis(mylm, matchCoefs(mylm, ":"), verbose=TRUE)
但我想一起测试交互术语和级别术语。这可能吗?请注意
mylm$coefficients
# (Intercept) qsec gear am qsec:am am:hp
# -12.2376256 0.8891289 4.1170265 -19.4050359 1.5298394 -0.0316123
具有qsec:am
而不是am:qsec
。然后
linearHypothesis(mylm, c("am + qsec:am + am:hp"))
确实有效,但这种排序并不明显。比如说,
lm(mpg ~ am:qsec + am:hp, data = mtcars)$coef
# (Intercept) am:qsec am:hp
# 17.1256930 0.7542508 -0.0456892
一旦我看到它,看起来很容易。非常感谢,我在文档中无法理解这一点。是的,正如我想演示的那样,命名似乎并没有遵循一个清晰的模式。
lm(mpg ~ am:qsec + am:hp, data = mtcars)$coef
# (Intercept) am:qsec am:hp
# 17.1256930 0.7542508 -0.0456892