Stata: 空间计量模型溢出效应的动态呈现

连享会-空间计量现场班 (2019年6月27-30日,西安)

点击此处-查看完整推文列表

这篇文章展示了如何创建动画图形,以展示空间自回归(SAR)模型生成的空间溢出效应。阅读本文后,您可以创建如下的动态图形。

连享会-德州犯罪率-空间溢出效应

后文包括三个部分:首先,介绍如何估计 SAR 模型的参数;其次,解释 SAR 模型为什么可以产生空间溢出效应;最后,展示如何创建一个动态图形来说明空间溢出效应。

SAR 模型

我想分析一下德克萨斯州各县的凶杀率与失业率的关系。我猜想一个县的凶杀率会对邻近县的凶杀率造成影响。

我想回答两个问题:

  1. 如何构建这样一个模型,在这个模型中明确的允许一个县的凶杀率依赖于邻近县的凶杀率。
  2. 根据这个模型,如果“达拉斯”(德克萨斯州的一个城市)的失业率上升到 10%,那么与达拉斯县邻近的县的凶杀率将如何变化呢?
建立SAR模型

首先,考虑标准线性模型,即某一个县 i i i 的凶杀率 ( h r a t e i ) ({\bf{hrate}}_{i}) (hratei)是该县失业率 ( u n e m p l o y m e n t i ) ({\bf unemployment}_{i}) (unemploymenti)线性函数

h r a t e i = β 0 + β 1 u n e m p l o y m e n t i + ϵ i {\bf hrate}_i = \beta_0 + \beta_1 {\bf unemployment}_{i} + \epsilon_i hratei=β0+β1unemploymenti+ϵi

SAR 模型构建了某一个县 i i i 的凶杀率 ( h r a t e i ) ({\bf hrate}_i) (hratei) 依赖于其邻近县的凶杀率。此时,需要一些新的符号来定义 SAR 模型。具体的,如果区域 j j j 和区域 i i i 相邻,则设 W i , j W_{i,j} Wi,j 为正数,如果两者不相邻,则 W i , j W_{i,j} Wi,j 为 0 。由于区域 i i i 不可能和自己相邻,所以当 i = j i = j i=j 时, W i , j W_{i,j} Wi,j 也为 0。

有了这个符号,构建某一个县 i i i 的凶杀率依赖于其邻近县凶杀率的 SAR 模型可以写成

h r a t e i = γ 1 ∑ j = 1 N W i , j h r a t e j + β 1 u n e m p l o y m e n t i + β 0 + ϵ i {\bf hrate}_i = \gamma_1\sum_{j=1}^N W_{i,j} {\bf hrate}_{j} + \beta_1 {\bf unemployment}_{i} + \beta_0 + \epsilon_i hratei=γ1j=1NWi

### Stata 中的空间溢出效应分析 #### 使用 `xsmle` 命令进行空间计量经济学建模 为了在 Stata 中处理和分析空间溢出效应,可以利用命令 `xsmle` 来执行空间计量经济模型的估计。此命令允许指定不同的空间权重矩阵以及多种类型的模型结构来捕捉区域间的相互影响。 ```stata * 加载数据集并定义变量 * use "your_dataset.dta", clear * 定义空间权重矩阵 W * spmat use W using "path_to_your_spatial_weights_file.gal" * 运行带有稳健标准误差的空间滞后模型(SAR) * xsmle ln_CO2 ln_ICT ln_URB ln_PGDP ln_IND ln_POP ln_EI, emat(W) model(sar) robust nolog noeffects ``` 上述代码展示了如何加载数据文件、读取预先准备好的空间权重矩阵,并通过 `xsmle` 执行一个稳健的标准差调整后的空间滞后模型(SAR)。这里假设因变量为二氧化碳排放 (`ln_CO2`) 而其他变量作为解释因素参与其中[^3]。 #### 处理动态空间面板数据的方法 当涉及到随时间和地点变化的数据时,则需考虑采用更复杂的动态空间面板模型来进行更为精确地描述现象之间的关系。针对这类情况,通常会应用到如下几种主要的技术: - 准似然函数或偏误修正的最大似然估计; - 广义矩估计(GMM); - 贝叶斯马尔科夫链蒙特卡罗模拟(Bayesian MCMC)[^2]。 这些技术能够更好地控制潜在的时间序列特性及其可能带来的内生性问题,在实际操作过程中可根据具体研究需求选取合适的方式实施计算。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值