最新地理加权回归(-GWR)

空间计量经济学打破大多数经典统计和计量分析中相互独立的基本假设，主要解决如何在横截面数据和面板数据的回归模型中处理空间相互作用（空间自相关）和空间结构（空间不均匀性）分析的问题。

空间计量经济理论认为一个地区空间单元上的某种经济地理现象或某一属性值与邻近地区空间单元上同一现象或属性值是相关的。

也就是说，各区域之间的数据存在与时间序列相关相对应的空间相关。

空间计量模型所研究的空间效应包括空间自相关和空间差异性。

空间相关性在空间回归模型中体现在误差项和因变量的滞后项，因此，空间计量的两个模型分别是空间自回归模型(Spatial Auto Regressive Model , SAR) 与空间误差模型(Spatial Error Model , SEM)，空间自回归模型研究各变量在一个地区是否有扩散效应，空间误差模型考察邻接地区关于因变量的误差冲击对本地区观察值的影响。

其表达式分别为：其中，Y 为因变量；W 为n n ⨯阶的空间权重矩阵，权数系数可以根据实际情况决定，一般用邻接矩阵；Wy 为空间滞后因变量，反映了空间距离对区域行为的作用；ρ为空间自回归系数，反映相邻区域的观测值Wy 对本地区观察值y 的影响方向和程度；X 为k n ⨯的外生解释变量向量(包括常数项)，β为变量系数，反映了自变量X 对因变量Y 的影响；ε为误差成分；λ为1⨯n 的因变量向量的空间误差系数，衡量了相邻地区的观察值Y 对本地区观察值Y 的影响方向和程度；γ为正态分布的随机误差向量。

上述两种模型的估计如果仍采用OLS ，往往导致各种结果和推论不够完整、科学。

本文采用极大似然法估计参数。

常用检验准则有拟合优度R 2 和对数似然值LogL 。

拟合优度和对数似然值越大，模型拟合效果越好, 对数似然值最大的模型最好。

( 一) 空间权重矩阵的选取空间权重矩阵 w 表征了空间单位之间的相互信赖性与关联程度。

实证研究中，通常采用相邻规则与距离规则来定义空间加权矩阵。

为了研究需要，本文从地理位置特征与社会经济特征两个不同角度分别建立包括相邻规则与距离规则的空间加权矩阵，以便更准确地把握房价的区域相关关系。

1. 地理位置特征加权矩阵。

本文采用两种常用的地理位置特征矩阵体现房价的空间相关关系:第一种是空间相邻加权矩阵 W1，其中的元素 wi ，j= 1表示两个地区拥有共同的边界，wi ，j= 0 表示两个地区没有共同的边界，然后对矩阵进行标准化处理。

为了避免“单个岛屿效应”，设定海南省与广东省、广西壮族自治区有共同边界。

第二种是空间距离加权矩阵 W2，其元素 wi ，j=1 / d2iji ≠ j0{i = j ，即两地区之间距离越远，相互之间的影响程度越小，两地区之间的距离 di ，j 为两地区省会城市之间的距离于是进一步用更能反映经济变量之间的空间依赖性的地理加权回归(GWR)方法，以全国30个省市为例,建立模型0(,)(,)i i i k i i ik i ky u v u v X ββε=++∑，1,2,i =…,n ，其中(,)i i u v 是第i 个样本点的空间坐标；利用加权最小二乘法来估计，估计出30个省市的模型参数，并就此分析了各省市之间的差异。

最后比较了普通回归与地理加权回归的优劣，得出了教育支出促进经济增长, 不同地区间促进的效果不同的结论。

长期以来, 在主流的经济学理论中, 空间事物无关联及均质性假定的局限,以及普遍使用忽视空间效应的普通最小二乘法(OL S) 进行模型估计, 使得在实际应用中往往存在模型的设定偏差问题, 进而导致经济学研究得出的各种结果和推论不够完整、科学, 缺乏应有的解释力。

经典计量经济学中的线性回归模型的经典假定, 以及回归模型的系数β是一个常数假定, 面对异常复杂的经济系统和因素变量之间的交互影响, 尤其是碰到横截面数据之间存在空间自相关性和空间异质性时, 经典计量的线性回归模型就显得有些力不从心, 需要发展新的方法来弥补这种不足。

【1】空间计量经济学(Anselin ,1988) 理论认为,一个地区空间单元上的某种经济地理现象或某一属性值与邻近地区空间单元上同一现象或属性值是相关的。

几乎所有的空间数据都具有空间依赖性或空间自相关性的特征, 空间依赖的存在打破了大多数经典统计和计量分析中相互独立的基本假设。

也就是说, 各区域之间的数据存在与时间序列相关、相对应的空间相关。

空间统计和空间计量经济方法是在继承和发展完善经典统计和计量方法的基础上, 将经典统计和计量方法应用于与地理位置及空间交互作用相关的地理空间数据, 通过地理位置与空间联系建立的统计与计量关系, 以统计和计量方法识别和度量空间变动的规律与空间模式的决定因素。

(一)空间经济计量学介绍空间统计和空间计量经济学理论与方法继承和发展了经典统计和计量理论方法，将经典统计和计量方法应用于与地理位置及空间交互作用相关的地理空间数据，通过地理位置与空间联系建立统计与计量关系，以统计和计量方法识别和度量空间变动规律及空间模式的决定因素。

空间经济计量学主要研究存在空间效应的问题。

空间效应主要包括空间相关和空间差异性。

在研究中涉及空间相邻、空间相邻矩阵等概念。

1.空间相关空间相关指在样本观测中，位于位置i的观测与其它j≠i的观测有关，即存在空间相关的原因有两方面：相邻空间单元存在测量误差，空间交互影响的存在。

测量误差是由于调查过程中，数据的采集与空间中的单位有关，如数据是按盛市、县等统计的，但设定的空间单位与研究问题不一致，存在测量误差。

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

2. 空间差异性空间差异性指空间上的区域缺乏均一性，如存在中心区和郊区、先进和后进地区等。

例如，我国沿海地区和中西部地区经济存在较大差别。

对于空间差异性，只要将空间单元的特性考虑进去，大多可以用经典经济计量学方法解决。

但当空间差异性与空间相关共同存在时，经典经济计量学方法不再适用，而且这时问题可能变得非常复杂，因为这时要区分空间差异性与空间相关可能非常困难。

3. 时空数据空间模型在模型中考虑时间维增加了描述的复杂性，但综合时间空间的模型在实际工作中非常有用。

在经典的经济计量学模型中，这是综合截面和时间序列数据的情形。

如果数据不存在空间相关，则可以采用PanelData 模型。

Anselin(1988)将似不相关(SUR)模型扩展到空间的情形，提出空间SUR 模型。

【2】(五) 空间计量经济学中的空间自相关分析【5】根据空间统计和空间计量经济学原理方法，首先应采用空间统计分析Moran 指数法检验因变量（被解释变量）是否存在空间自相关性或集聚现象，如果存在，则需要在空间计量经济学理论方法支持下，建立空间计量经济模型，进行区域经济增长集聚的空间计量估计和检验。

Moran’s I 定义如下： 11211()()n n ij i j i j n niji j W Y Y Y Y I S W ====--=∑∑∑∑ 其中，211()n i i S Y Y n ==-∑，11ni i Y Y n ==∑， 表示第i 地区的观测值，n 为地区总数，ij W 为二进制的邻接空间权值矩阵，表示其中的任一元素，采用邻接标准或距离标准，其目的是定义空间对象的相互邻接关系。

一般邻接标准的为1,ij W ⎧=⎨⎩ 当区域i 和区域j 相邻；0， 当区域i 和区域j 不相邻；式中，1,2,i =…，n ；j=1,2,?…,n ；≠m=n 或m n 。

Moran ’s I 可看作各地区观测值的乘积和，其取值范围为11I -≤≤。

若各地区间经济行为为空间正相关，I 的数值应当较大；负相关则较小。

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

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

本文即主要采用的这种模型对全国30个省市教育与经济增长之间的关系进行了探究与分析。

1. GWR 基本模型考虑如下的全局回归模型:0i k ik i ky X ββε=++∑ 1,2,i =…,n地理加权回归( GWR) 扩展了传统的回归框架, 容许局部而不是全局的参数估计, 扩展后模型的参数是位置i 的函数，扩展后的模型如下:0(,)(,)i i i k i i ik i k y u v u v X ββε=++∑ 1,2,i =…,n其中(,)i i u v 是第i 个样本点的空间坐标, (,)k i i u v β 是连续函数(,)k u v β 在i点的值。

如果(,)k u v β 在空间保持不变, 则GWR 模型就变为全局模型 。

因此GWR 方程认可空间变化关系可能是存在的, 并且提供了一种可度量的方法。

由上面可知GWR 模型中的参数在每个回归点是不同的,就不能用最小二乘方法(OL S) 估计参数。

Fotheringham , Brunsdon , Charlton (1996) 依据“接近位置i 的观察数据比那些离位置远一些的数据对(,)k u v β的估计有更多的影响”的思想, 利用加权最小二乘法来估计参数。

因此,其结果是区域性的并非全域性的参数估计,从而就能够探测到空间数据的空间非平稳性。

我们知道，普通最小二乘法可以得到全局的参数估计向量：1(')'X X X Y β∧-=值得提到的是，使用最小二乘估计的前提条件是：()E Y X β=，2()(n n Var Y I I δ=是单位矩阵）成立。

这里第二个条件不满足，可以改为21()i Var Y W δ-=，因为1i W ->0,存在n 阶非奇异对称阵B ，使得12i W B -=。

令1*Y B Y -=，1*X B X -=，则11(*)()*E Y B E Y B X X ββ--===112(*)()n Var Y B Var Y B I δ--==于是，我们得到回归点i 的参数估计向量可以表示如下：1(,)('(,))'(,)i i i i i i u v X W u v X X W u v Y β∧-=其中(,)i i W u v 是n n ⨯的加权矩阵,对角线上的每个元素都是关于观测值所在位置j 与回归点i 的位置之间距离的函数,其作用是权衡不同空间位置j ( j = 1 ,2 , ⋯, n) 的观测值对于回归点i 参数估计的影响程度,而非对角元素为0. 矩阵(,)i i W u v 可以表示为如下形式：12W W (,) W i i i i in W u v ⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭=O 记做12W (W ,W ,W )i i i in diag =…,。