抑制 R 中的固定效应系数答案

【问题标题】：Suppress fixed effects coefficients in R抑制 R 中的固定效应系数
【发布时间】：2014-12-21 08:26:31
【问题描述】：

在使用 summary() 函数时，有没有办法抑制线性模型中固定效应的系数（例如，stata 中的吸收函数的等效项）。例如，我希望汇总函数只输出截距和 x，而不是因子的系数和标准误：

frame <- data.frame(x = rnorm(100), y = rnorm(100), z = rep(c("A", "B", "C", "D"),25))
summary(lm(y~x + as.factor(z), data = frame)) 

Call:
lm(formula = y ~ x + as.factor(z), data = frame)

Residuals:
Min      1Q  Median      3Q     Max 
-2.2417 -0.6407  0.1783  0.5789  2.4347 

Coefficients:
          Estimate Std. Error t value Pr(>|t|)  
(Intercept)   -0.25829    0.19196  -1.346   0.1816  
 x              0.09983    0.09788   1.020   0.3104

...

谢谢。

【问题讨论】：

标签： r summary lm

【解决方案1】：

我强烈推荐stargazer 包：

frame <- data.frame(x = rnorm(100), y = rnorm(100), z = rep(c("A", "B", "C", "D"),25))
summary(lm(y~x + as.factor(z), data = frame)) 
model <- lm(y~x + as.factor(z), data = frame)
library(stargazer)
stargazer(model, type="text")
stargazer(model, omit = "z", type="text")

【讨论】：

【解决方案2】：

您可以执行以下任何操作：

mod <- lm(y~x + as.factor(z), data = frame)
coef(mod)[c("(Intercept)", "x")]
# (Intercept)           x 
#  0.12357491 -0.06430765 
coef(mod)[grepl("Intercept|x", names(coef(mod)))]
# (Intercept)           x 
#  0.12357491 -0.06430765 
coef(mod)[1:2]
# (Intercept)           x 
#  0.12357491 -0.06430765 
mod$coefficients[1:2]
# (Intercept)           x 
#  0.12357491 -0.06430765

【讨论】：

有什么可以显示整个输出的吗？比仅仅获取系数更通用的东西？
@robertevansanders - 我回答了手头的问题。 OP想要拦截和x。如果您有任何问题，您可以向社区提出。感谢您的反对。
我会改变它，但他的问题不连贯。首先他说他只想抑制固定效应但有总结，然后他说他只想有截距和系数。我认为他的问题意味着第一件事，因为这就是我要找的。我一直在寻找如何做第一件事，但这些问题的许多答案都没有做到，直到我找到了“lfe”包。
似乎最好通过在问题下发表评论然后等待回复来消除这种混淆。 OP检查了这个答案，所以它解决了他的直接问题。不太确定你想让我怎么做。
我的错，我会改回来的。（做一个小的编辑，这样我就可以了）我只是对没有做第一部分的答案投了反对票。出乎意料地难以找到......！