从具有固定效应的 lm() 的摘要中排除变量答案

【问题标题】：Excluding variables from summary of lm() with Fixed Effects从具有固定效应的 lm() 的摘要中排除变量
【发布时间】：2015-11-21 17:55:33
【问题描述】：

假设我估计了下面的具有固定效应的OLS模型，总结如下：

m <- lm(x ~ y + factor(z)) 
summary(m)

如果 z 足够大，则汇总输出会很长，因为它将 z 的每个元素的信息显示为系数。有没有办法合并这个输出，这样我就没有 z 的每个元素的系数，而是简单地看到“FE：是”或者只是从表中排除 FE？我正在尝试找到一个通用的解决方案。

【问题讨论】：

请提供一些示例数据。 x、y 和 z 是什么？
“FE：是”是什么意思？您的意思是您只想查看输出中具有统计意义的估计值吗？只有具有一定绝对量级的影响？只有具有一定效果大小的效果？就个人而言，我认为自动化统计决策是非常不明智的......

【解决方案1】：

以下是您的请求的一种可能解释：

summary(m)$coefficients[ !grepl("z", names(coef(m)) ) , ,drop=FALSE] # tested only on example from ?lm

如果不是你想要的，那么请做一个更完整的描述，以及一个具有所需输出的测试用例。

【讨论】：

【解决方案2】：

使用 stargazer 包打印模型摘要，并使用 omit 选项压下输出中的因子名称。要包含表示固定效果 = true 的行，您可以使用选项 Additional.lines 将其包含在表格的底部。 Stargazer 可以输出类型为 html、text 或 latex。

【讨论】：

【解决方案3】：

我认为您确实想简化 z 中的关卡。

FE = z == "FE"
lm(x ~ y + factor(FE))

【讨论】：