如何从lm摘要输出中省略NA系数？

Question

我有许多定性相互作用的回归，有些组合没有任何观察结果。如何从summary.lm()输出中没有观察到排除交互系数？

Answer 1

如果不破解summary.lm()方法就很难（不可能？），但如果你愿意只看系数表，你可以使用na.omit(coef(summary(...))))。

示例：

set.seed(101)
d <- data.frame(y=rnorm(100),x1=rnorm(100))
## add redundant variables
d$x2 <- d$x3 <- d$x1
m1 <- lm(y~.,data=d)  ## fit y to all other vars

summary(m1)
na.omit(coef(summary(m1)))
##               Estimate Std. Error    t value Pr(>|t|)
## (Intercept) -0.0329789 0.09341769 -0.3530263 0.724827
## x1           0.1002849 0.09341569  1.0735337 0.285668