1、第1页共7页空间计量经济学分析空间依赖、空间异质性*传统的统计理论是一种建立在独立观测值假定基础上的理论。然而，在现实世界中，特别是遇到空间数据问题时，独立观测值在现实生活中并不是普遍存在的（Getis,197）。*对于具有地理空间属性的数据，一般认为离的近的变量之间比在空间上离的远的变量之间具有更加密切的关系（Anseli&Getis，192）。正如著名的Toblr地理学第一定律所说：“任何事物之间均相关，而离的较近事物总比离的较远的事物相关性要高。”（le，197）*地区之间的经济地理行为之间一般都存在一定程度的SpatilInteraction，SpatilEfects）：SpatilD
2、epndceandSpatilAutocrelation）。*一般而言，分析中涉及的空间单元越小，离的近的单元越有可能在空间上密切关联（Anli&Getis,192）。*然而，在现实的经济地理研究中，许多涉及地理空间的数据，由于普遍忽视空间依赖性，其统计与计量分析的结果值得进一步深入探究（Anseli&Grifn,198）。*可喜的是，对于这种地理与经济现象中常常表现出的空间效应（特征）问题的识别估计，空间计量经济学提供了一系列有效的理论和实证分析方法。*一般而言，在经济研究中出现不恰当的模型识别和设定所忽略的空间效应主要有两个来源（Anseli，198）：空间依赖性（SpatilDepndc
3、e）和空间异质性（SpatilHetrogneity）。空间依赖性*空间依赖性（也叫空间自相关性）是空间效应识别的第一个来源，它产生于空间组织观测单元之间缺乏依赖性的考察（Clif&Ord,1973）。*Anseli&Rey（19）区别了真实（Substanil）空间依赖性和干扰（Nuisance）空间依赖性的不同。真实空间依赖性反映现实中存在的空间交互作用（SpatilInteractionEfet），*比如区域经济要素的流动、创新的扩散、技术溢出等，它们是区域间经济或创新差异演变过程中的真实成分，是确确实实存在的空间交互影响，*如劳动力、资本流动等耦合形成的经济行为在空间上相互影响、相互作
4、用，研发的投入产出行为及政策在地理空间上的示范作用和激励效应。*干扰空间依赖性可能来源于测量问题，比如区域经济发展过程研究中的空间模式与观测单元之间边界的不匹配，造成了相邻地理空间单元出现了测量误差所导致。*测量误差是由于在调查过程中，数据的采集与空间中的单位有关，如数据一般是按照省市县等行政区划统计的，这种假设的空间单位与研究问题的实际边界可能不一致，这样就很容易产生测量误差。*空间依赖不仅意味着空间上的观测值缺乏独立性，而且意味着潜在于这种空间相关中的数据结构，也就是说空间相关的强度及模式由绝对位置（格局）和相对位置（距离）共同决定。*空间相关性表现出的空间效应可以用以下两种模型来表征和刻
5、画：当模型的误差项在空间上相关时，即为空间误差模型；当变量间的空间依赖性对模型显得非常关键而导致了空间相关时，即为空间滞后模型（Anseli，198）。空间异质性*空间异质性（空间差异性），是空间计量学模型识别的第二个来源。空间异质性或空间差异性，指地理空间上的区域缺乏均质性，存在发达地区和落后地区、中心（核心）和外围（边缘）地区等经济地理结构，从而导致经济社会发展和创新行为存在较大的空间上的差异性。*空间异质性反映了经济实践中的空间观测单元之间经济行为（如增长或创新）关系的一种普遍存在的不稳定性。*区域创新的企业、大学、研究机构等主体在研发行为上存在不可忽视的个体差异，譬如研发投入的差异导致
6、产出的技术知识的差异，*这种创新主体的异质性与技术知识异质性的耦合将导致创新行为在地理空间上具有显著的异质性差异，进而可能存在创新在地理空间上的相互依赖现象或者创新的局域俱乐部集团。*对于空间异质性，只要将空间单元的特性考虑进去，大多可以用经典的计量经济学方法进行估计。但是当空间异质性与空间相关性同时存在时，经典的计量经济学估计方法不再有效，而且在这种情况下，第2页共7页问题变得异常复杂，区分空间异质性与空间相关性比较困难。*空间变系数的地理加权回归模型（GeographiclWeightdRegrsion，简记为GWR）是处理空间异质性的一种良好的估计方法。空间依赖性及集群的识别与检验*将空
7、间依赖性考虑进来以后，在建立模型进行创新分析研究之前，一般必须先进行空间相关性的预检验如果空间效应在发挥作用，则需要将空间效应纳入模型分析框架之中，并采用适合于空间计量经济模型估计的方法进行估计*如果没有表现出空间效应，则可直接采用一般估计方法（如OLS）估计模型参数。在引入空间变量或者经过空间过滤的空间计量模型建立之后，其效果的好坏还需要通过空间相关检验进行判断，一般可通过对真实值和模型估计值之间的残差进行空间相关性检验实现。*如果参数经过检验在空间上没有表现出相关性，则表明在引入空间变量或者经过考虑了空间效应后的模型已经成功地处理了空间相关性。*计算和检验一个地区的创新行为在地理空间上有没
8、有表现出空间自相关（依赖）性（Clif&Ord，198；Cresi，193）、是否存在集群现象（吴玉鸣，206a，206b）。*空间统计和空间计量经济学的方法有许多种，最著名也最常用的有：MoransI（Moran，1950）、Gearys、Gti指数（Ord&Getis，195）。*以上方法各有其功用，各有其优缺点，同时亦有其适用范围及限制。一般来说，这些方法在功用上可大致分为两大类：全域空间自相关（GlobalSptialAutocrelation）和局域空间自相关（LocalSptialAutocrelation）。全域空间相关性检验与分析*全域空间自相关（Glbltiltcrlti）是
9、从区域空间的整体上刻画区域创新活动空间分布的集群情况。*在实际的空间相关分析应用研究中，由于MoransI指数和GearysC比率的作用基本相同，其不同之处是MoransI主要针对于全域空间相关性分析，而erys指数则适用于局域空间关联分析。*在许多实证研究中，MoransI和GearysC是常用方法，已在大量文献中出现，尤其是前者。因此，以下介绍常用的rsI指数的计算及检验过程。*oransI定义如下：*其中，表示第个地区的观测值（如专利数），为地区总数（如21()niiSY=1niiYY=i i省域），为二进制的邻近空间权值矩阵，表示其中的任一元素，采用邻近标准或距离标准，其目的是定义空间
10、对象的相互邻近关系。*如果MoransI的正态统计量的Z值均大于正态分布函数在0.5（.01）水平下的临界值1.65（.96），表明区域创新在空间分布上具有明显的正向相关关系，*正的空间相关代表相邻地区的类似特征值出现集群（Clustering）趋势。目前计量研究方法及其局限性*目前有关研究的计量方法主要是传统的回归分析方法（如多元统计分析、回归分析、数据包络分析DEA等方法），其实质上都是线性的变量之间相互关系的一种测量方法，适合于企业或产业部门时间序列层面的经验研究，未考虑区域（或截面单元）之间的空间关联，局限性比较明显。*区域之间的经济行为会相互影响，这使其存在显著的外部效应，导致地区之
11、间的经济行为存在溢出效应。经济产出不仅受到本地投资的强度、而且还会受到周边其他地区的投资活动产生的溢出效应及政策的影响。*经济集群行为可以通过检验一个代表地区间考虑交互作用的生产函数，即该地区的经济活动及其自身的121()() niji jijnijijWYYMoransIS=第3页共7页特征与其他地区的经济活动的关系，来考察区域经济行为的集群行为。*可通过纳入空间效应（空间相关和空间差异）的空间计量经济模型空间回归模型，包括常系数的空间滞后模型（SpatilLgModel，SL）与空间误差模型两种（SpatilEroMdel，SE）（Anseli，198；Anseli，Florx，195；吴
12、玉鸣，205）和变系数的地理加权回归模型（GographiclWightdRgrin，GWR）来实现。空间计量经济模型及估计技术*空间计量经济学是计量经济学的一个分支，研究的是在横截面数据（Cros-ectionalDt）和面板数据（PanelDta）的回归模型中如何处理空间交互作用（空间自相关）和空间结构（空间非均匀性）（Anseli，198）。*最近发展起来的空间统计学和空间计量经济学不仅解决了标准统计方法在处理空间数据时的失误问题,更重要的是为测量这种空间联系及其性质、并在建模时明确地引入空间联系变量以估算与检验其贡献提供了全新的手段（应龙根，宁越敏，205）。*实际上，早在1970年代
13、欧洲就展开了空间计量经济学研究，并将它作为一个确定的领域。Paelinck&Klasen（）定义了这个领域，包括：空间相互依赖在空间模型中的任务；空间关系不对称性；位于其他空间的解释因素的重要性；过去的和将来的相互作用之间的区别；明确的空间模拟。*Ali（198）将空间计量经济学定义为：处理由区域科学模型统计分析中的空间所引起的特殊性的技术总称。换句话说，空间计量经济学研究的是明确考虑空间影响（空间自相关和空间不均匀性）的方法。*目前，空间计量经济学研究包括以下四个感兴趣的领域：计量经济模型中空间效应的确定；*合并了空间影响的模型的估计；空间效应存在的说明、检验和诊断；*空间预测。空间计量经济
14、学模型有多种类型（Anseli，tal.204）。*首先介绍纳入了空间效应（空间相关和空间差异）、适用于截面数据的空间常系数回归模型，包括空间滞后模型（SpatilLgModel，SLM）与空间误差模型（SpatilEroMdel，SEM）两种，以及空间变系数回归模型地理加权回归模型（GeographiclWeighdRegsin，GWR）。*适用于时间序列和截面数据合成的空间面板数据计量经济学模型将在以后予以介绍。空间滞后模型*空间滞后模型（SpatilLgModel，SLM）主要是探讨各变量在一地区是否有扩散现象（溢出效应）。其模型表达式为：yWyX=+*参数反映了自变量对因变量的影响，空
15、间滞后因变量是一内生变量，反映了空间距离对区域行为 Wy的作用。*区域行为受到文化环境及与空间距离有关的迁移成本的影响，具有很强的地域性（Anselital.，196）。由于SLM模型与时间序列中自回归模型相类似，因此SLM也被称作空间自回归模型（SpilAutoregsiveodl，SAR）。空间误差模型E*空间误差模型（patilroMdel，SEM）的数学表达式为：yX=+W*式中，为随机误差项向量，为n*1阶的截面因变量向量的空间误差系数，为正态分布的随机误差 向量。*SEM中参数反映了自变量X对因变量y的影响。参数衡量了样本观察值中的空间依赖作用，即相邻 地区的观察值y对本地区观察值
16、y的影响方向和程度。*存在于扰动误差项之中的空间依赖作用，度量了邻近地区关于因变量的误差冲击对本地区观察值的影响程度。由于SEM模型与时间序列中的序列相关问题类似，也被称为空间自相关模型（SpatilAutocrelationdel，SAC）。第4页共7页估计技术*鉴于空间回归模型由于自变量的内生性，对于上述两种模型的估计如果仍采用OLS，系数估计值会有偏或者无效，需要通过IV、ML或GS、M等其他方法来进行估计。*Anseli（198）建议采用极大似然法估计空间滞后模型（SLM）和空间误差模型（EM）的参数。空间自相关检验与S、E的选择*判断地区间创新产出行为的空间相关性是否存在，以及和E那
17、个模型更恰当，一般可通过包括MoransI检验、两个拉格朗日乘数（LagrneMultipler）形式LR、LMAG及其稳健（Robust）的R-LE、R-LMAG）等形式来实现。*由于事先无法根据先验经验推断在S和E模型中是否存在空间依赖性，有必要构建一种判别准则，以决定哪种空间模型更加符合客观实际。*Anseli和Florax（195）提出了如下判别准则：如果在空间依赖性的检验中发现LMAG较之LMER在统计上更加显著，且R-LMAG显著而R-LME不显著，则可以断定适合的模型是空间滞后模型；相反，如果LMER比L在统计上更加显著，且-LR显著而R-LMAG不显著，则可以断定空间误差模型是
18、恰当的模型。*除了拟合优度2检验以外，常用的检验准则还有：自然对数似然函数值（Loglikelihod，Log）、似然比率（LikelihodRatio，L）、赤池信息准则（Akaieinformatincritin，AIC）、施瓦茨准则（Schwartzcritin，SC）。*对数似然值越大，AI和值越小，模型拟合效果越好。这几个指标也用来比较OLS估计的经典线性回归模型和LM、E，似然值的自然对数最大的模型最好。空间变系数回归模型及估计*就目前国内外的研究来看，大多直接假定横截面单元是同质的，即地区或企业之间没有差异。传统的OLS只是对参数进行“平均”或“全域”估计，不能反映参数在不同空间
19、的空间非稳定性（吴玉鸣，李建霞，206；苏方林，207）。*当用横截面数据建立计量经济学模型时，由于这种数据在空间上表现出的复杂性、自相关性和变异性，使得解释变量对被解释变量的影响在不同区域之间可能是不同的，假定区域之间的经济行为在空间上具有异质性的差异可能更加符合现实。*空间变系数回归模型（SpatilVrying-CoeficntRegrsionMdel）中的地理加权回归模型（GeographiclWeightdRegsi，GW）是一种解决这种问题的有效方法。地理加权回归模型G*地理加权回归模型是一种相对简单的回归估计技术，它扩展了普通线性回归模型。在扩展的GR模型中，特定区位的回归系数不
20、再是利用全部信息获得的假定常数，而是利用邻近观测值的子样本数据信息进行局域（Local）回归估计而得的、随着空间上局域地理位置变化而变化的变数，W模型可以表示为：()()k01, ,i ii jiiijijyuvuvx=+*式中，系数的下标j表示与观测值联系的阶待估计参数向量，是关于地理位置的j 1m (),iiuvk+1元函数。GWR可以对每个观测值估计出k个参数向量的估计值，是第i个区域的随机误差，满足零均值、同方差、相互独立等球形扰动假定。*实际上，以上模型可以表示为在每个区域都有一个对应的估计函数，其对数似然函数可以表示为：()()()()0 2021 1log,|, ,kn ki i
21、i iiii jLuvuvMyvuvx= = + *式中，为常数，。由于极大似然法（ML）的解不是(),2,1,2,iijiiMyxuvnjk = 唯一的，Hastie和Tibshirani（193）认为用该方法求解是不恰当的。Tibshirani和Hstie（1987）提出了局域求解法，原理与方法如下：第5页共7页*对于第s个空间位置，任取一空间位置与其位置邻近，构造一个简单的(),1,2,suvsn= ()0,uv回归模型：01ki jiijyx=+*式中，每个为常数且为GWR模型中的近似值，通过考虑与点相邻近的点来校正经j (),jsuv ()0,uv典回归模型中的解。一个基本的方法就是
22、采用加权最小二乘法（WLS），寻找合适的使得下式最小：() 2001 1n kii jii jWdyx= = *式中，为位置和之间的空间距离，0i ()0,uv(),iiuv为空间权值。令为的估计值，可得GWR模型在空间位置上的估计值()0id j(),jsv (),suv。对上式求的一阶偏导数，并令其等于0，可得：()()() 01, ,s s ksuvuvu j()()2 20 0jXWXWY =*式中，为的对角线矩阵。可以看出，0 ()()01020,nddWd j(1,2,)k=的GR估计值是随着空间权值矩阵的变化而变化的，因此的选择至关重要，一般由观测值的空ij ijW间（经纬度）坐
23、标决定。*实际研究中常用的空间距离权值计算公式有三种（LeSag，204）。（1）高斯距离权值（GausinDistance）()/ij ijWd=*（2）指数距离权值（Exponetialistnce）()exp/ij ijdq=*（3）三次方距离权值（TricubeDistance）()31/ij ijWd=*其中为第i个区域与第j个区域间的地理距离,为标准正态分布密度函数，q为观测值i到第q个最近邻ijd 居之间的距离，为距离向量的标准差，为衰减参数（窗宽）。*在空间权值矩阵中，d和非常关键。如果d较大，则局域模型的解越趋向于全域模型的解；如果d等于所研究空间任意两点间的最大距离，则全域
24、和局域两个模型将相等，反之则相反。*若趋于无穷大，任意两点的权重将趋于1，则被估计的参数变成一致时，GWR就等于以OLS估计的经典线性回归；反之，当带宽变得很小时，参数估计将更加依赖于邻近的观测值。计算适当的窗宽或衰减函数的原理方法很多，最小二乘法仍然是一般常用的方法，其原理是：()210niiiDy=*式中，是用窗宽计算所得的的拟合值。()i iy空间权值矩阵的确定*在区域经济管理研究中，将空间效应因素引入经济管理过程的研究，建立空间计量经济模型进行空间统计分析时，一般要用空间权值矩阵来表达空间相互作用。*空间权值矩阵是一种与被解释变量的空间自回归过程相联系的矩阵。在实际的区域分析中，该矩阵
25、的选择设定是外生的，原因是n维的W包含了关于区域i和区域j之间相关的空间连接的外生信息，不需要通过模型来估计得到它，只需通过权值计算出来就行了。第6页共7页*W中对角线上的元素被设为0，而表示区域i和区域j在空间上相连接的原因。为了减少或消除区域ijWijW间的外在影响，权值矩阵被标准化（）成行元素之和为1。对于变量x，这种转换意味1nijij ijjww=着定义成空间滞后变量的*仅仅表示邻近观测值的加权平均数。xW空间权值矩阵W确定的方法有多种，根据空间统计和空间计量经济学原理，一般可将现实的地理空间关联或者经济联系考虑到模型中来，以达到正确设定权值矩阵的目的。其中，衡量地理联系的方法通常有
26、两种主要方式：邻近指标和距离指标。*按照这两种方法确定的为二进制的邻近空间权值矩阵，表示其中的任一元素，采用邻近标准或距离ij标准，其目的是定义空间对象的相互邻近关系，便于把地理信息系统数据库中的有关属性放到所研究的地理空间上来对比。*可用矩阵表示如下：ijW111212122212nmmmnwwww = 基于邻近概念的空间权值矩阵*根据相邻标准，为：ijW10ij jWi=当区域和区域相邻；当区域和区域不相邻；*式中，；或。1,2,1,2,injm= ；n=mn基于邻近概念的空间权值矩阵（ContiguityBasedSpatilWeights）有一阶邻近矩阵和高阶邻近矩阵两种。*一阶邻近矩
27、阵（theFirstOrdetiitMtrix）是假定两个地区有共同边界时空间关联才会发生，即当相邻地区i和j有共同的边界用1表示，否则以0表示。一般有Rok邻近和Quen邻近两种计算方法（Anseli，203）。*Rok邻近定义为仅有共同边界来定义邻居，而Quen邻近则除了共有边界邻区外还包括共同顶点的邻居。由此可见，基于Quen邻近的空间矩阵常常与周围地区具有更加紧密的关联结构（拥有更过的邻区）。*当然，如果假定区域间公共边界的长度不同（如10km和10k）其空间作用的强度也不一样，则还可以通过将共有边界的长度纳入权值计算过程中，使这种邻近指标更加准确一些。*空间矩阵不仅仅局限于第一阶邻近
28、矩阵，也可以计算和使用更高阶的邻近矩阵。Anseli&Smirnov（196）提出了高阶邻近矩阵的算法，其目的是为了消除在创建矩阵时出现的冗余及循环。*二阶邻近矩阵（theScondOreContiguityMatrix）表示了一种空间滞后的邻近矩阵。也就是说，该矩阵表达了邻近地区的相邻地区的空间信息。*当使用时空数据并假设随着时间推移产生空间溢出效应时，这种类型的空间权值矩阵将非常有用。在这种情况下，特定地区的初始效应或随机冲击将不仅会影响其邻近地区，而且随着时间的推移还会影响其邻近地区的相邻地区。当然，这种影响是几何递减的。*可以看出，邻近空间权值矩阵因其对称与计算简单而最为常用，适合于测
29、算地理空间效应的影响。K值最邻近空间矩阵（K-NearsteighborSpatilWeights）*Anseli（203）介绍了一种值最邻近空间矩阵（K-NesighborSpailWeights）。之所以提出这种距离矩阵，主要是因为一般使用的基于门槛距离（TldDitnc）的简单空间矩阵常常会导致一种非常不平衡的邻近矩阵结构。*譬如，在空间单元的面积相差甚大的情况下，就会出现小一些的地理单元具有很多邻近单元，而较大的第7页共7页地理单元则可能很少有邻近单元，甚至没有邻近单元而成为“飞地”。K值最邻近空间矩阵（K-NearsteighborSpatilWeights）*在这种情况下，考虑最近
30、邻居是一种可供选择的常用方法，这也是创建空间距离权值矩阵的第二种选择。*一般在给定空间单元周围选择最邻近的4个单元（亦可选4个以上，根据实际的空间关联情况由研究者确定），来计算K值最近邻居权值的大小。基于距离的空间权值矩阵*根据距离标准，为：ijW1()0ij jdijWdi=当区域和区域在距离之内（即区域和区域相邻）；当区域和区域在距离之外（即区域和区域不相邻）；*基于距离的空间权值矩阵（DistanceBasdSpatilWeights）方法是假定空间相互作用的强度是决定于地区间的质心距离或者区域行政中心所在地之间的距离，是一种在实践应用中常用的空间权值矩阵。*在这种情况下，不同的权值指标
31、随距离ij的定义而变化（如大城市圈、公路之间的距离1ijijwd=），其取值取决于选定的函数形式（如距离的倒数或倒数的平方，以及欧氏距离等）。21ij ijwd=*当然，还需要定义一个门槛距离，超过了某给定的门槛距离则区域间的相互作用可以忽略不计。如果输入的时空数据库中有x、y经纬度坐标数据，可以通过x、y坐标计算两点（两个地区的质心）之间的距离而获得空间权值矩阵。坐标的度量有欧氏距离（EuclideanDistance）和弧度距离（ArcDistance）两种，度量坐标系上任意两点间的距离可以通过具有地理坐标（坐标、y坐标）的变量的点来计算。*值得注意的是，对于经过投影的地理坐标只能计算欧氏
32、距离，而未经投影的经纬坐标适合于计算弧度距离。经济社会空间权值矩阵*除了使用真实的地理坐标计算地理距离外，还有包括经济和社会因素的更加复杂的权值矩阵设定方法。比如，根据区域间交通运输流、通讯量、GDP总额、贸易流动、资本流动、人口迁移、劳动力流等确定空间权值，计算各个地区任何两个变量之间的距离。空间权值矩阵的选择*从理论上来看，较之邻近矩阵，距离矩阵在空间效应测算中应该是比较科学和理想的一个指标。但是，在实际应用中，这种方法实行起来比较困难，原因有二：*一是社会经济距离的实际统计数据难以获得；二是模型中权值的计算是外生的。*当然，基于经济、社会因素的权值计算方法更加接近区域经济的现实，因而在数
33、据可得和模型结构清晰的情况下，可以考虑选择这种类型的权值。*尽管二进制的空间邻近权值矩阵并非适用于所有的空间计量经济模型，但是，处于某些情况下的实用性，空间统计学家在构建空间计量模型时的首选就是从二进制的邻近矩阵开始的。*一般是先从空间邻近的最基本二进制矩阵开始，逐步选择确定空间权值矩阵。关于各种权值矩阵的选择，没有现成的理论根据，一般可考虑空间计量模型对各种空间权值矩阵的适用程度，检验估计结果对权值矩阵的敏感性，最终的依据实际上就是结果的客观性和科学性。*本研究所有具体应用中涉及的空间邻近矩阵的计算，主要采用Anseli（19，203）研制开发的空间统计分析软件GeoDa095i，直接生成邻近矩阵来测算并确定地区之间的空间效应。
