空间计量经济学分析

空间计量经济学分析

空间依赖、空间异质性

? ?

?

? ? ? ?

传统的统计理论是一种建立在独立观测值假定基础上的理论。然而，在现实世界中，特别是遇到空间数据问题时，独立观测值在现实生活中并不是普遍存在的（Getis, 1997）。

对于具有地理空间属性的数据，一般认为离的近的变量之间比在空间上离的远的变量之间具有更加密切的关系（Anselin & Getis，1992）。正如著名的Tobler地理学第一定律所说：“任何事物之间均相关，而离的较近事物总比离的较远的事物相关性要高。”（Tobler，1979）

地区之间的经济地理行为之间一般都存在一定程度的Spatial Interaction，Spatial Effects）：Spatial Dependence and Spatial Autocorrelation）。

一般而言，分析中涉及的空间单元越小，离的近的单元越有可能在空间上密切关联（Anselin & Getis, 1992）。

然而，在现实的经济地理研究中，许多涉及地理空间的数据，由于普遍忽视空间依赖性，其统计与计量分析的结果值得进一步深入探究（Anselin & Griffin, 1988）。

可喜的是，对于这种地理与经济现象中常常表现出的空间效应（特征）问题的识别估计，空间计量经济学提供了一系列有效的理论和实证分析方法。

一般而言，在经济研究中出现不恰当的模型识别和设定所忽略的空间效应主要有两个来源（Anselin，1988）：空间依赖性（Spatial Dependence）和空间异质性（Spatial Heterogeneity）。

空间依赖性

? 空间依赖性（也叫空间自相关性）是空间效应识别的第一个来源，它产生于空间组织观测单元之间缺乏? ? ? ? ? ? ? ? ?

依赖性的考察（Cliff & Ord, 1973）。

Anselin & Rey（1991） 区别了真实（Substantial）空间依赖性和干扰（Nuisance）空间依赖性的不同。 真实空间依赖性反映现实中存在的空间交互作用（Spatial Interaction Effects）， 比如区域经济要素的流动、创新的扩散、技术溢出等，

它们是区域间经济或创新差异演变过程中的真实成分，是确确实实存在的空间交互影响，

如劳动力、资本流动等耦合形成的经济行为在空间上相互影响、相互作用，研发的投入产出行为及政策在地理空间上的示范作用和激励效应。

干扰空间依赖性可能来源于测量问题，比如区域经济发展过程研究中的空间模式与观测单元之间边界的不匹配，造成了相邻地理空间单元出现了测量误差所导致。

测量误差是由于在调查过程中，数据的采集与空间中的单位有关，如数据一般是按照省市县等行政区划统计的，这种假设的空间单位与研究问题的实际边界可能不一致，这样就很容易产生测量误差。

空间依赖不仅意味着空间上的观测值缺乏独立性，而且意味着潜在于这种空间相关中的数据结构，也就是说空间相关的强度及模式由绝对位置（格局）和相对位置（距离）共同决定。

空间相关性表现出的空间效应可以用以下两种模型来表征和刻画：当模型的误差项在空间上相关时，即为空间误差模型；当变量间的空间依赖性对模型显得非常关键而导致了空间相关时，即为空间滞后模型（Anselin，1988）。

空间异质性

? 空间异质性（空间差异性），是空间计量学模型识别的第二个来源。

? 空间异质性或空间差异性，指地理空间上的区域缺乏均质性，存在发达地区和落后地区、中心（核心）? ? ? ? ?

和外围（边缘）地区等经济地理结构，从而导致经济社会发展和创新行为存在较大的空间上的差异性。 空间异质性反映了经济实践中的空间观测单元之间经济行为（如增长或创新）关系的一种普遍存在的不稳定性。

区域创新的企业、大学、研究机构等主体在研发行为上存在不可忽视的个体差异，譬如研发投入的差异导致产出的技术知识的差异，

这种创新主体的异质性与技术知识异质性的耦合将导致创新行为在地理空间上具有显著的异质性差异，进而可能存在创新在地理空间上的相互依赖现象或者创新的局域俱乐部集团。

对于空间异质性，只要将空间单元的特性考虑进去，大多可以用经典的计量经济学方法进行估计。 但是当空间异质性与空间相关性同时存在时，经典的计量经济学估计方法不再有效，而且在这种情况下，

第 1 页 共 7 页

问题变得异常复杂，区分空间异质性与空间相关性比较困难。

? 空间变系数的地理加权回归模型（Geographical Weighted Regression，简记为GWR）是处理空间异质性的一种良好的估计方法。

空间依赖性及集群的识别与检验

? 将空间依赖性考虑进来以后，在建立模型进行创新分析研究之前，一般必须先进行空间相关性的预检验 ? 如果空间效应在发挥作用，则需要将空间效应纳入模型分析框架之中，并采用适合于空间计量经济模型? ? ? ? ? ? ?

估计的方法进行估计

如果没有表现出空间效应，则可直接采用一般估计方法（如OLS）估计模型参数。

在引入空间变量或者经过空间过滤的空间计量模型建立之后，其效果的好坏还需要通过空间相关检验进行判断，一般可通过对真实值和模型估计值之间的残差进行空间相关性检验实现。

如果参数经过检验在空间上没有表现出相关性，则表明在引入空间变量或者经过考虑了空间效应后的模型已经成功地处理了空间相关性。

计算和检验一个地区的创新行为在地理空间上有没有表现出空间自相关（依赖）性（Cliff & Ord，1981；Cressie，1993）、是否存在集群现象（吴玉鸣，2006a，2006b）。 空间统计和空间计量经济学的方法有许多种，最著名也最常用的有：Moran’s I（Moran，1950）、Geary’s C、Getis指数（Ord & Getis，1995）。

以上方法各有其功用，各有其优缺点，同时亦有其适用范围及限制。

一般来说，这些方法在功用上可大致分为两大类：全域空间自相关（Global Spatial Autocorrelation）和局域空间自相关（Local Spatial Autocorrelation）。

全域空间相关性检验与分析

? 全域空间自相关（Global Spatial Autocorrelation）是从区域空间的整体上刻画区域创新活动空间分布的

集群情况。

? 在实际的空间相关分析应用研究中，由于Moran’s I指数和Geary’s C比率的作用基本相同，其不同之处是Moran’s I主要针对于全域空间相关性分析，而Geary’s C指数则适用于局域空间关联分析。

? 在许多实证研究中，Moran’s I 和Geary’s C是常用方法，已在大量文献中出现，尤其是前者。因此，以下介绍常用的Moran’s I指数的计算及检验过程。 ? Moran’s I定义如下： nn

??Wij(Yi?Y)(Yj?Y) i?1j?1Moran'sI? nn2 S??Wij i?1j?1

S??

21nn其中， ， ，i 表示第i个地区的观测值（如专利数），为地区总数（如省域），为二进制的邻近空间权值矩阵，表示其中的任一元素，采用邻近标准或距离标准，其目的是定义空间对象的相互邻近关系。

? 如果Moran’s I的正态统计量的Z值均大于正态分布函数在0.05（0.01）水平下的临界值1.65（1.96），表明区域创新在空间分布上具有明显的正向相关关系，

? 正的空间相关代表相邻地区的类似特征值出现集群（Clustering）趋势。

i?1i?1?(Yi?Y)Y?1nn?YiY目前计量研究方法及其局限性

? 目前有关研究的计量方法主要是传统的回归分析方法（如多元统计分析、回归分析、数据包络分析DEA

等方法），其实质上都是线性的变量之间相互关系的一种测量方法，适合于企业或产业部门时间序列层面的经验研究，未考虑区域（或截面单元）之间的空间关联，局限性比较明显。 ? 区域之间的经济行为会相互影响，这使其存在显著的外部效应，导致地区之间的经济行为存在溢出效应。 ? 经济产出不仅受到本地投资的强度、而且还会受到周边其他地区的投资活动产生的溢出效应及政策的影响。

? 经济集群行为可以通过检验一个代表地区间考虑交互作用的生产函数，即该地区的经济活动及其自身的

第 2 页 共 7 页

特征与其他地区的经济活动的关系，来考察区域经济行为的集群行为。

? 可通过纳入空间效应（空间相关和空间差异）的空间计量经济模型——空间回归模型，包括常系数的空间滞后模型（Spatial Lag Model，SLM）与空间误差模型两种（Spatial Error Model，SEM）（Anselin，1988；Anselin，Florax，1995；吴玉鸣，2005）和变系数的地理加权回归模型（Geographical Weighted Regression，GWR）来实现。

空间计量经济模型及估计技术

?

?

?

?

? ? ? ? ?

? ?

?

空间计量经济学是计量经济学的一个分支，研究的是在横截面数据（Cross-sectional Data）和面板数据（Panel Data）的回归模型中如何处理空间交互作用（空间自相关）和空间结构（空间非均匀性）（Anselin，1988）。

最近发展起来的空间统计学和空间计量经济学不仅解决了标准统计方法在处理空间数据时的失误问题,更重要的是为测量这种空间联系及其性质、并在建模时明确地引入空间联系变量以估算与检验其贡献提供了全新的手段（应龙根，宁越敏，2005）。

实际上，早在1970年代欧洲就展开了空间计量经济学研究，并将它作为一个确定的领域。Paelinck & Klaassen（1979）定义了这个领域，包括：空间相互依赖在空间模型中的任务；空间关系不对称性；位于其他空间的解释因素的重要性；过去的和将来的相互作用之间的区别；明确的空间模拟。

Anselin（1988）将空间计量经济学定义为：处理由区域科学模型统计分析中的空间所引起的特殊性的技术总称。换句话说，空间计量经济学研究的是明确考虑空间影响（空间自相关和空间不均匀性）的方法。

目前，空间计量经济学研究包括以下四个感兴趣的领域： 计量经济模型中空间效应的确定； 合并了空间影响的模型的估计； 空间效应存在的说明、检验和诊断； 空间预测。

空间计量经济学模型有多种类型（Anselin，et al. 2004）。

首先介绍纳入了空间效应（空间相关和空间差异）、适用于截面数据的空间常系数回归模型，包括空间滞后模型（Spatial Lag Model，SLM）与空间误差模型（Spatial Error Model，SEM）两种，以及空间变系数回归模型——地理加权回归模型（Geographical Weighted Regression，GWR）。 适用于时间序列和截面数据合成的空间面板数据计量经济学模型将在以后予以介绍。

空间滞后模型（Spatial Lag Model，SLM）主要是探讨各变量在一地区是否有扩散现象（溢出效应）。其模型表达式为：y??Wy?X???

空间滞后模型SLM

? ?

参数?反映了自变量对因变量的影响，空间滞后因变量Wy是一内生变量，反映了空间距离对区域行为的作用。

? 区域行为受到文化环境及与空间距离有关的迁移成本的影响，具有很强的地域性（Anselin et al.，1996）。 ? 由于SLM模型与时间序列中自回归模型相类似，因此SLM也被称作空间自回归模型（Spatial Autoregressive Model，SAR）。

空间误差模型SEM

? 空间误差模型（Spatial Error Model，SEM）的数学表达式为：

y?X???

???W???

? 式中，?为随机误差项向量，?为n*1阶的截面因变量向量的空间误差系数，?为正态分布的随机误差

向量。

? SEM中参数?反映了自变量X对因变量y的影响。参数?衡量了样本观察值中的空间依赖作用，即相邻地区的观察值y对本地区观察值y的影响方向和程度。

? 存在于扰动误差项之中的空间依赖作用，度量了邻近地区关于因变量的误差冲击对本地区观察值的影响程度。由于SEM模型与时间序列中的序列相关问题类似，也被称为空间自相关模型（Spatial Autocorrelation Model，SAC）。

第 3 页 共 7 页

估计技术

? 鉴于空间回归模型由于自变量的内生性，对于上述两种模型的估计如果仍采用OLS，系数估计值会有偏

或者无效，需要通过IV、ML或GLS、GMM等其他方法来进行估计。

? Anselin（1988）建议采用极大似然法估计空间滞后模型（SLM）和空间误差模型（SEM）的参数。

空间自相关检验与SLM、SEM的选择

?

?

?

?

?

判断地区间创新产出行为的空间相关性是否存在，以及SLM和SEM那个模型更恰当，一般可通过包括Moran’s I检验、两个拉格朗日乘数（Lagrange Multiplier）形式LMERR、LMLAG及其稳健（Robust）的R-LMERR、R-LMLAG）等形式来实现。

由于事先无法根据先验经验推断在SLM和SEM模型中是否存在空间依赖性，有必要构建一种判别准则，以决定哪种空间模型更加符合客观实际。

Anselin和Florax（1995）提出了如下判别准则：如果在空间依赖性的检验中发现LMLAG较之LMERR在统计上更加显著，且R-LMLAG显著而R-LMERR不显著，则可以断定适合的模型是空间滞后模型；相反，如果LMERR比LMLAG在统计上更加显著，且R-LMERR显著而R-LMLAG不显著，则可以断定空间误差模型是恰当的模型。

除了拟合优度R2检验以外，常用的检验准则还有：自然对数似然函数值（Log likelihood，LogL）、似然比率（Likelihood Ratio，LR）、赤池信息准则（Akaike information criterion，AIC）、施瓦茨准则（Schwartz criterion，SC）。

对数似然值越大，AIC和SC值越小，模型拟合效果越好。这几个指标也用来比较OLS估计的经典线性回归模型和SLM、SEM，似然值的自然对数最大的模型最好。

空间变系数回归模型及估计

? 就目前国内外的研究来看，大多直接假定横截面单元是同质的，即地区或企业之间没有差异。

? 传统的OLS只是对参数进行“平均”或“全域”估计，不能反映参数在不同空间的空间非稳定性（吴玉鸣，李

建霞，2006；苏方林，2007）。

? 当用横截面数据建立计量经济学模型时，由于这种数据在空间上表现出的复杂性、自相关性和变异性，使得解释变量对被解释变量的影响在不同区域之间可能是不同的，假定区域之间的经济行为在空间上具有异质性的差异可能更加符合现实。

? 空间变系数回归模型（Spatial Varying-Coefficient Regression Model）中的地理加权回归模型（Geographical Weighted Regression，GWR）是一种解决这种问题的有效方法。

地理加权回归模型GWR

? 地理加权回归模型是一种相对简单的回归估计技术，它扩展了普通线性回归模型。

? 在扩展的GWR模型中，特定区位的回归系数不再是利用全部信息获得的假定常数，而是利用邻近观测

值的子样本数据信息进行局域（Local）回归估计而得的、随着空间上局域地理位置变化而变化的变数，GWR模型可以表示为：

kyi??0?ui,vi?????u,v?xjiij?1ij??i

? 式中，?j系数的下标j表示与m?1观测值联系的阶待估计参数向量，是关于地理位置?ui,vi? 的

k+1元函数。GWR可以对每个观测值估计出k个参数向量的估计值，?是第i个区域的随机误差，满足零均值、同方差、相互独立等球形扰动假定。

? 实际上，以上模型可以表示为在每个区域都有一个对应的估计函数，其对数似然函数可以表示为：

logL?L???0?u,v?,?,?k?u,v?|M????12?2n2 ??k?ui,vi?xi?????i?1??yi??0?ui,vi???k?j?1?? 式中，?为常数，M???yi,xij,?ui,vi?,i?1,2,?,n,j?1,2,?,k?。由于极大似然法（ML）的解不是

唯一的，Hastie和Tibshirani（1993）认为用该方法求解是不恰当的。Tibshirani和Hastie（1987）提出

了局域求解法，原理与方法如下：

第 4 页 共 7 页

? 对于第s个空间位置???us,vs?,s?1,2,?,n??，任取一空间位置?u0,v0?与其位置邻近，构造一个简单的

回归模型：

kyi??0???j?1jxij??i

? 式中，每个?j为常数且为GWR模型中?j?us,vs?的近似值，通过考虑与点?u0,v0?相邻近的点来校正经

典回归模型中的解。一个基本的方法就是采用加权最小二乘法（WLS），寻找合适的使得下式最小：

?Wdyi??0????0i?i?1?n??x?jij? j?1?k2? 式中，d0i为位置?u0,v0?和?ui,vi?之间的空间距离，

??u,v?的估计值，可得GWR模型在空间位置?u,v?上的估计值? W?d0i?为空间权值。令 ??j为?jssss???0?us,vs?,??1?us,vs?,?,??k?us,vs??。对上式求?j的一阶偏导数，并令其等于0，可得：

2??j??X?W0X??X?W?120Y?

?? 式中，W0为??W?d01,?,W?d02,?,?,W?d0n???的对角线矩阵。可以看出，?j(j?1,2,?,k)

的GWR估计值是随着空间权值矩阵Wij的变化而变化的，因此Wij的选择至关重要，一般由观测值的空间（经纬度）坐标决定。

? 实际研究中常用的空间距离权值计算公式有三种（LeSage，2004）。

? （1）高斯距离权值（Gaussian Distance）

Wij???dij/???

? （2）指数距离权值（Exponential Distance）

Wij?exp??dij/q?

3? （3）三次方距离权值（Tricube Distance）

3? Wij?1???/dij?? ????? 其中dij为第i个区域与第j个区域间的地理距离,?为标准正态分布密度函数，q为观测值i到第q个最近邻

居之间的距离，?为距离向量的标准差，?为衰减参数（窗宽）。

? 在空间权值矩阵中，d和θ非常关键。如果d较大，则局域模型的解越趋向于全域模型的解；如果d等于所研究空间任意两点间的最大距离，则全域和局域两个模型将相等，反之则相反。

? 若θ趋于无穷大，任意两点的权重将趋于1，则被估计的参数变成一致时，GWR就等于以OLS估计的经典线性回归；反之，当带宽变得很小时，参数估计将更加依赖于邻近的观测值。计算适当的窗宽或衰减函数的原理方法很多，最小二乘法仍然是一般常用的方法，其原理是：

n2iD????yi?1?i?????y??0

?i???是用窗宽θ计算所得的yi的拟合值。 ? 式中，y空间权值矩阵的确定

? 在区域经济管理研究中，将空间效应因素引入经济管理过程的研究，建立空间计量经济模型进行空间统

计分析时，一般要用空间权值矩阵来表达空间相互作用。

? 空间权值矩阵是一种与被解释变量的空间自回归过程相联系的矩阵。在实际的区域分析中，该矩阵的选择设定是外生的，原因是n×n维的W包含了关于区域i和区域j之间相关的空间连接的外生信息，不需要通过模型来估计得到它，只需通过权值计算出来就行了。

第 5 页 共 7 页

本文来源：https://www.bwwdw.com/article/lo5h.html