如何在Stata中制作这样的系数估计图？

Question

graph from Cantoni 2014 如何在Stata中绘制这样的系数估计？

这是一个横截面的面板数据：

city year gdp church A 1996 2000 23 A 1997 2200 25 A 1998 2500 30 B 1996 2800 20 ...

我想看看gdp对教堂数量的影响，我使用xtreg church gdp, fe robust的固定效果模型 如何制作系数估算图？

谢谢！

Answer 1

1. 对＆＃34;横截面面板数据的任何引用＆＃34;是错的。当在每个时间点对或面板数据采样不同的单位时，数据集可以是 横截面数据。从上下文我推断这是一个面板数据集。

2. 教堂数量为正离散变量，而OLS则为连续响应变量。这可能会导致问题，如越界预测。所以我建议使用泊松回归。

3. 国内生产总值经常非常倾斜，因此我会记录日志。

查看以下代码。

我生成的数据集在church和log_gdp之间有一些正相关，然后是预测计数的图表：

clear
set obs 10000
egen year = seq(), from(1990) to(2015)
bysort year: egen city = seq(), from(1) to(400)
gen gdp = rnormal(3000,1000)
gen log_gdp=log(gdp)
gen disturb=1+int((10)*runiform())
gen church = int(log_gdp+disturb)

xtset city year
xtpoisson church log_gdp, fe robust
margins, at(log_gdp=(0.35(0.1)8.8)) predict(iru0)
marginsplot

产生这个：