我目前正在一个项目中,在该项目中,我正在比较多个线性回归。 Stargazer在导出相关模型的摘要方面做得非常出色,但是我一直在想,是否有可能在stargazer命令中实现一个功能,以同时导出多个摘要并以图形方式对其进行排序。我当前的方法如下:
首先,我使用以下方法将相关摘要导出为html格式:
stargazer(development_gr_US_pre_BvG, type = "html", out = "whatever")
stargazer(development_gr_US_post_BvG, type = "html", out = "whatever")
相应的输出是:
<table style="text-align:center"><tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"></td><td><em>Dependent variable:</em></td></tr>
<tr><td></td><td colspan="1" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td>n</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">earliest_filing_year</td><td>1.510<sup>*</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.760)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>-3,007.075<sup>*</sup></td></tr>
<tr><td style="text-align:left"></td><td>(1,523.425)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>14</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.248</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.185</td></tr>
<tr><td style="text-align:left">Residual Std. Error</td><td>11.463 (df = 12)</td></tr>
<tr><td style="text-align:left">F Statistic</td><td>3.947<sup>*</sup> (df = 1; 12)</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
</table>
<table style="text-align:center"><tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"></td><td><em>Dependent variable:</em></td></tr>
<tr><td></td><td colspan="1" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td>n</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">earliest_filing_year</td><td>-4.400<sup>**</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.648)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>8,872.900<sup>**</sup></td></tr>
<tr><td style="text-align:left"></td><td>(1,304.897)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>4</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.958</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.938</td></tr>
<tr><td style="text-align:left">Residual Std. Error</td><td>1.449 (df = 2)</td></tr>
<tr><td style="text-align:left">F Statistic</td><td>46.095<sup>**</sup> (df = 1; 2)</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
</table>
然后我将输出传输到Notepad ++,并通过以下方式对其进行格式化:
<table style="float: left;">
<tr><td></td><td colspan="1" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td>n</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">earliest_filing_year</td><td>1.510<sup>*</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.760)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>-3,007.075<sup>*</sup></td></tr>
<tr><td style="text-align:left"></td><td>(1,523.425)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>14</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.248</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.185</td></tr>
<tr><td style="text-align:left">Residual Std. Error</td><td>11.463 (df = 12)</td></tr>
<tr><td style="text-align:left">F Statistic</td><td>3.947<sup>*</sup> (df = 1; 12)</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
<table style="float: left;">
<tr><td></td><td colspan="1" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td>n</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">earliest_filing_year</td><td>-4.400<sup>**</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.648)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>8,872.900<sup>**</sup></td></tr>
<tr><td style="text-align:left"></td><td>(1,304.897)</td></tr>
<tr><td style="text-align:left"></td><td></td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>4</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.958</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.938</td></tr>
<tr><td style="text-align:left">Residual Std. Error</td><td>1.449 (df = 2)</td></tr>
<tr><td style="text-align:left">F Statistic</td><td>46.095<sup>**</sup> (df = 1; 2)</td></tr>
<tr><td colspan="2" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
</table>
但是那只是我坎bump的解决方法。我一直在想是否有一种很好的方法告诉R立即以这种方式存储多个摘要,或者甚至更好的是立即以这种格式将它们作为PNG提供给我。谢谢您的指导!
stargazer
可以轻松获取更多模型并排显示。在控制台中显示为“文本”,或导出为“ tex”或“ html”文件。因此,您只看到一次协变量标签,对于模型统计数据也是如此。
library("plm")
library("stargazer")
data("Produc", package = "plm")
# Regression
model1 <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
data = Produc,
index = c("state","year"),
method="pooling")
model2 <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp),
data = Produc,
index = c("state","year"),
method="pooling")
stargazer(model1, model2, type = "html", out="models.htm")
输出看起来像这样: