如何從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