在R的stargazer中如何省略固定效应模型的因子变量系数?
解决Stargazer输出固定效应模型时隐藏城市虚拟变量的问题
当你用lm()拟合包含城市固定效应的模型后,默认用stargazer()输出会把所有城市虚拟变量都列出来,结果显得杂乱又冗长。要只保留截距项和pop_growth的系数,其实只需要一个简单的参数设置就能搞定——就是用omit="city"来过滤掉所有带"city"关键词的虚拟变量。
完整操作步骤
- 先拟合你的固定效应模型:
m1 <- lm(GDP ~ pop_growth + factor(city), data = your_data)
- 调用
stargazer()时加上omit参数:
stargazer(m1, omit = "city", type = "text") # 按需把type改成"html"或"latex"
原理说明
stargazer的omit参数支持部分匹配变量名,而我们用factor(city)生成的虚拟变量,变量名都是city1、city2这种带"city"前缀的格式,所以指定omit="city"就能精准过滤掉所有城市固定效应的虚拟变量,只留下你关心的截距和pop_growth系数。
内容的提问来源于stack exchange,提问作者Mackenzie




