从具有固定效果的 lm() 摘要中排除变量
Excluding variables from summary of lm() with Fixed Effects
假设我用固定效应估计下面的 OLS 模型并总结如下:
m <- lm(x ~ y + factor(z))
summary(m)
如果 z 足够大,则汇总输出会非常冗长,因为它将 z 的每个元素的信息显示为系数。有没有办法合并这个输出,这样我就没有 z 的每个元素的系数,而是简单地看到 "FE: Yes" 或简单地从 table 中排除 FE?我正在尝试找到一个通用的解决方案。
以下是对您的请求的一种可能解释:
summary(m)$coefficients[ !grepl("z", names(coef(m)) ) , ,drop=FALSE]
# tested only on example from ?lm
如果不是您想要的,请提供更完整的描述以及具有所需输出的测试用例。
我认为你真的想简化 z 中的关卡。
FE = z == "FE"
lm(x ~ y + factor(FE))
使用 stargazer 包打印模型摘要,并使用 omit 选项压低输出中的因子名称。要包含一行 fixed effects = true,您可以使用选项 additional.lines 将其包含在 table 的底部。 Stargazer 可以将类型输出为 html、文本或乳胶。
假设我用固定效应估计下面的 OLS 模型并总结如下:
m <- lm(x ~ y + factor(z))
summary(m)
如果 z 足够大,则汇总输出会非常冗长,因为它将 z 的每个元素的信息显示为系数。有没有办法合并这个输出,这样我就没有 z 的每个元素的系数,而是简单地看到 "FE: Yes" 或简单地从 table 中排除 FE?我正在尝试找到一个通用的解决方案。
以下是对您的请求的一种可能解释:
summary(m)$coefficients[ !grepl("z", names(coef(m)) ) , ,drop=FALSE]
# tested only on example from ?lm
如果不是您想要的,请提供更完整的描述以及具有所需输出的测试用例。
我认为你真的想简化 z 中的关卡。
FE = z == "FE"
lm(x ~ y + factor(FE))
使用 stargazer 包打印模型摘要,并使用 omit 选项压低输出中的因子名称。要包含一行 fixed effects = true,您可以使用选项 additional.lines 将其包含在 table 的底部。 Stargazer 可以将类型输出为 html、文本或乳胶。