如何在R中使用Stargazer进行多次回归?

时间:2016-02-09 21:55:35

标签: r stargazer dummy-variable

我正在尝试使用R中的Stargazer包创建一个回归表。 我有几个回归,只有虚拟变量不同。 我希望它报告自变量的系数,常数等,并且如果在回归中包括某些固定效应(即虚拟变量),则说“是”或“否”。 这些是我的回归:

m1 <- lm(data=merge1,log(total_units)~log(priceIndex))
m2 <- lm(data=merge1,log(total_units)~log(priceIndex)+factor(fips_state_code))
m3 <- lm(data=merge1,log(total_units)~log(priceIndex)+factor(fips_state_code)+factor(month))
m4 <- lm(data=merge1,log(total_units)~log(priceIndex)+factor(fips_state_code)+factor(month)+factor(year))
m5 <- lm(data=merge1,log(total_units)~log(priceIndex)+factor(fips_state_code)+factor(month)+time*factor(fips_state_code))
m6 <- lm(data=merge1,log(total_units)~log(priceIndex)+factor(fips_state_code)+factor(month)+factor(year)+time*factor(fips_state_code))

我尝试过stargazer命令的几种变体来让表格看起来正确,但它从来没有。 当我运行此命令时,对于五个或六个回归中的每一个,它按预期工作:

stargazer(m2,type="text",
          omit = c("fips_state_code","month","year","time"),
          omit.labels = c("State FE?","Month of year FE?","Year FE?","State time trend?"))

即,对于m2,它在“状态FE”旁边显示“是”?所有其他问题旁边都有“否”。对于m3,它在“状态FE”旁边显示“是”?和“一年中的FE?”和其他问题旁边的“否”。

但是当我运行此命令时,表格会针对所有回归的所有问题报告“否”:

stargazer(m1,m2,m3,m4,m5,m6,type="text",
          omit = c("fips_state_code","month","year","time"),
          omit.labels = c("State FE?","Month of year FE?","Year FE?","Time FE?"))

有人知道发生了什么吗?无论我是单独还是一起做各回归,它应该是一样的吗?不是吗?

我也得到了其他奇怪的结果......当我运行以下内容时:

stargazer(m3,m4,m5,m6,type="html",
          omit = c("fips_state_code","month","year","time"),
          omit.labels = c("State FE?","Month of year FE?","Year FE?","State time trend?"))

我明白了:

<table style="text-align:center"><tr><td colspan="5" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"></td><td colspan="4"><em>Dependent variable:</em></td></tr>
<tr><td></td><td colspan="4" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td colspan="4">log(total_units)</td></tr>
<tr><td style="text-align:left"></td><td>(1)</td><td>(2)</td><td>(3)</td><td>(4)</td></tr>
<tr><td colspan="5" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">log(priceIndex)</td><td>2.962<sup>***</sup></td><td>-0.746<sup>***</sup></td><td>0.142</td><td>-1.947<sup>***</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.206)</td><td>(0.249)</td><td>(0.224)</td><td>(0.276)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>-0.094</td><td>11.248<sup>***</sup></td><td>8.570<sup>***</sup></td><td>15.030<sup>***</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.652)</td><td>(0.779)</td><td>(0.709)</td><td>(0.868)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td></tr>
<tr><td colspan="5" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">State FE?</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td></tr>
<tr><td style="text-align:left">Month of year FE?</td><td>Yes</td><td>Yes</td><td>Yes</td><td>Yes</td></tr>
<tr><td style="text-align:left">Year FE?</td><td>No</td><td>No</td><td>No</td><td>No</td></tr>
<tr><td style="text-align:left">State time trend?</td><td>No</td><td>No</td><td>No</td><td>No</td></tr>
<tr><td colspan="5" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>2,853</td><td>2,853</td><td>2,853</td><td>2,853</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.968</td><td>0.974</td><td>0.984</td><td>0.985</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.967</td><td>0.974</td><td>0.983</td><td>0.984</td></tr>
<tr><td style="text-align:left">Residual Std. Error</td><td>0.340 (df = 2806)</td><td>0.305 (df = 2800)</td><td>0.244 (df = 2771)</td><td>0.235 (df = 2766)</td></tr>
<tr><td style="text-align:left">F Statistic</td><td>1,830.438<sup>***</sup> (df = 46; 2806)</td><td>2,019.427<sup>***</sup> (df = 52; 2800)</td><td>2,041.885<sup>***</sup> (df = 81; 2771)</td><td>2,084.570<sup>***</sup> (df = 86; 2766)</td></tr>
<tr><td colspan="5" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td colspan="4" style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
</table>

0 个答案:

没有答案