追踪我国制造集聚的空间来源:基于马歇尔外部性与新经济地理的综合视角 中国区蜮经济论坛 表1我国地级及以上城市各空间变量 变量法或广义矩(GMM)估计。但由于工具变量法 其他变量的样本统计值 小值工最大值在很大程度上依赖于工具变量的选取,而对于不同 造业集聚水平 变量和模型,很难找到合适的工具变量,这势必影 FDI存量万元) 15011504216895161.227252019716 响到模型的稳健性。因此,本文采用面板GMM进 厂商数(家 50276112450338.00001847400 32691328274709202374382463 行估计。面板GMM估计方法更适合于“大N小T” 6.332253.1313 6478139.955 特征的微观数据,其估计偏误在T(时间)固定情况 均路平方*人817551.7040m下,随着N(截面)的增加而减小,本文采用200 万人电话数(台/万人 18340114739473737745260726 2009年284个地级及以上城市的面板数据,样本结 专业化劳动力可得性(d=1)2200136136530000730481 投人可得性(=1) 构符合以上特征。另外,面板广义矩估计法有差分 技术外溢I(=1 147281286109349义矩估计(DFCM)和系统广义矩估计(sYs 空间技术外溢Ⅱ(6=1) 0005763900660000786103930 国内市场潜力(8=1 181391890162+246M)之分,由于差分CMM估计量的有限样本特性 1173008497050195364518243 较差,尤其是当滞后项和随后的一阶差分项存在非 中投人可200141210101常弱的相关性时,工具变量较弱(Rodm00 空间技术外溢1(=241073719460403004343713971 此时采用系统广义矩估计更有效、偏差也更小。以 3体3Fmm下估计结果均基于系统GM (二)空间外部性作用范围的检验 100公里、100~300公里、300-500公里、500-1000公 确定空间外部性作用范围有助于把握城市间 里、10002000公里、20003000公里以及3000公里不同经济联系随距离变化的空间特征。由于各城 以上各范围内的空间外部性变量,以检验和确定各 市国际市场潜力不受国内距离影响,因而需要检验 空间变量的作用范围。表1报告了我国地级及以上 除国际市场潜力外其他空间变量的作用范围。为 城市要素供给范围在100公里以内、商品市场遍布 与现有文献有可比性,我们使用0-50公里、50~100 全国的各空间变量及其他变量的样本统计值 公里、100~300公里、300-500公里、500-1000公里、 五、计量检验与结果分析 1000-2000公里、2000~3000公里以及3000公里以 上各范围内的数据根据(16)~(20)式构建空间变量 (一)计量策略 并对方程(14)进行估计。本文采用面板系统GMM 进行计量检验之前首先要选择适宜的面板数 模型估计了分别包含距离衰减参数为1或2的空间 据模型。由于本文基本模型的解释变量中含有被变量(包括空间供给外部性和市场潜力)的计量方 解释变量的滞后项,即使假定随机扰动项不存在序程。估计结果中,衰减参数为2的空间变量参数估 列相关,被解释变量滞后项与非观测固定效应也可计的显著性及模型拟合程度均不及衰减参数为1时 能存在相关性。另外,根据区域和城市经济学原的情况,表2中报告了衰减参数为1时不同距离范 理,不仅地区存在的专业化劳动力、中间投入的规 围内各空间变量的估计结果。 模经济、技术外溢和市场需求规模有利于制造业集 表2的各个模型中,系数联合显著性的Wad检 聚,而且制造业在同一地区的集中分布对该地区与验均在1%水平上拒绝了解释变量系数为0的原假 其他地区的空间联系也具有加强作用。因而空间设。 Sargan和 Hansen统计量对应的p值均大于 供给和需求的外部性与制造业集聚之间可能存在10%,不能拒绝工具变量有效的原假设。Arla 联立内生性。出于数据可得性,我们在设置计量模 no- Bond AR(1)统计量显示各计量方程残差的一阶 型时也遗漏了诸如自然条件、资源禀赋等变量,而差分项存在一阶负相关,但AR(2)统计量的伴随概 这些变量可能致使相关解释变量与随机扰动项之率均在10%以上,因此不能拒绝原模型中残差无自 间存在相关性。无论使用固定效应、随机效应模型相关的假设。可见,工具变量的选择是合适的,模 还是FGIS模型,可能都会导致估计系数有偏和不型估计结果是可取的 致。为了得到无偏、一致的估计量,可采用工具 就各控制变量而言,各方程中上一期的制造业 62i l 琮我国制 i t 此集聚的全间来源：基于马跃余外部性与淅经涔地理的综合祝 角 中国区城经涔论 坛 表 1 我国地级及以上城市各空间变 量 及其他变量的样本统计 值 变量 均值 标 准 差 最小值 最 大 值 制造业集聚水 平 0.32958 0.26315 0.00135 1.98945 劳均工资(元） 17611.37 6606.099 1895.148 134432.1 FDI存量(万元） 1501150 4216895 161.2272 52019716 厂商数(家） 502.761 1245.033 8.00000 18474.00 土地价格(万元/公顷） 326.913 282.747 0.920237 4382.463 财政收人比重(％) 6.33225 3.13135 0.464781 39.95523 人均道路面积(平方米/人） 8.13779 5.584187 0.78000 64.00000 中学及以上在校人数比重(％) 21.38746 11.06082 3.061811 88.42524 万人电话数(台/万人） 18340.1 14739.47 3737.745 260726.2 专业化劳动力可得性 (5=1 ) 2.280012 3.613653 0.0000109 77.30481 中间投入可得性(5=1 ) 0.048195 0.041116 0.0017866 0.396793 空间技术外溢 1^=1 ) 731.4728 2429.826 0.2079904 36449.09 空间技术外溢 n 0.0057639 0.0036036 0.0000786 0.039308 国内市场潜力(<y=i) 731881.5 337990.2 189270.3 6247246 国际市场潜力(<y=i) 117300.8 497701.5 0.1195364 5182437 专业化劳动力可得性(<y=2 ) 0.2351961 0.7510085 0.00000119 27.28197 中间投人可得性垆 = 2 ) 0.0041811 0.0021023 0.0004111 0.0198886 空间技术外溢 1^=2 ) 41.07737 94.60403 0.0443437 1397.176 空间技术外溢n(<y=2 ) 0.0008132 0.0008568 0.0000149 0.0078361 国内市场潜力(<y=2) 14918.41 13696.44 1392.958 398470.5 国际市场潜力( 5 = 2) 103396.59 366273.37 0.04180367 4175988.53 100 公里、100-300 公里、300~500 公里、500~ 1000 公 里、1000-2000公里、2000-3000公里以及3000公里 以上各范围内的空间外部性变量，以检验和确定各 空间变量的作用范围。表1报告了我国地级及以上 城市要素供给范围在100公里以内、商品市场遍布 全国的各空间变量及其他变量的样本统计值。 五、计量检验与结果分析 (一）计量策略 进行计量检验之前首先要选择适宜的面板数 据模型。由于本文基本模型的解释变量中含有被 解释变量的滞后项，即使假定随机扰动项不存在序 列相关，被解释变量滞后项与非观测固定效应也可 能存在相关性。另外，根据区域和城市经济学原 理，不仅地区存在的专业化劳动力、中间投入的规 模经济、技术外溢和市场需求规模有利于制造业集 聚，而且制造业在同一地区的集中分布对该地区与 其他地区的空间联系也具有加强作用。因而空间 供给和需求的外部性与制造业集聚之间可能存在 联立内生性。出于数据可得性，我们在设置计量模 型时也遗漏了诸如自然条件、资源禀赋等变量，而 这些变量可能致使相关解释变量与随机扰动项之 间存在相关性。无论使用固定效应、随机效应模型 还是FGLS模型，可能都会导致估计系数有偏和不 一致。为了得到无偏、一致的估计量，可采用工具 变量法或广义矩（GMM)估计。但由于工具变量法 在很大程度上依赖于工具变量的选取，而对于不同 变量和模型，很难找到合适的工具变量，这势必影 响到模型的稳健性。因此，本文采用面板GMM进 行估计。面板 GMM 估计方法更适合于 “ 大 N 小 T ” 特征的微观数据，其估计偏误在T(时间）固定情况 下，随着N(截面）的增加而减小，本文采用2003~ 2009年284个地级及以上城市的面板数据，样本结 构符合以上特征。另外，面板广义矩估计法有差分 广义矩估计（DIF GMM) 和系统广义矩估计（ S Y S GMM)之分，由于差分GMM估计量的有限样本特性 较差，尤其是当滞后项和随后的一阶差分项存在非 常弱的相关性时，工具变量较弱（Roodman，2006)， 此时采用系统广义矩估计更有效、偏差也更小。以 下估计结果均基于系统GMM。 (二）空间外部性作用范围的检验 确定空间外部性作用范围有助于把握城市间 不同经济联系随距离变化的空间特征。由于各城 市国际市场潜力不受国内距离影响，因而需要检验 除国际市场潜力外其他空间变量的作用范围。为 与现有文献有可比性，我们使用0~50公里、50~100 公里、100-300 公里、300-500 公里、500-1000 公里、 1000-2000公里、2000~3000公里以及3000公里以 上各范围内的数据根据（16)~(20)式构建空间变量 并对方程（14)进行估计。本文采用面板系统G M M 模型估计了分别包含距离衰减参数为1或2的空间 变量（包括空间供给外部性和市场潜力）的计量方 程。估计结果中，衰减参数为2的空间变量参数估 计的显著性及模型拟合程度均不及衰减参数为1时 的情况，表2中报告了衰减参数为1时不同距离范 围内各空间变量的估计结果®。 表2的各个模型中，系数联合显著性的Wald检 验均在 1 % 水平上拒绝了解释变量系数为 0 的 原 假 设 。 Sargan 和 Hansen 统计量对应的 p 值均大 于 m %，不能拒绝工具变量有效的原假设。Arella - no-Bond AR(1)统计量显示各计量方程残差的一阶 差分项存在一阶负相关，但AR(2)统计量的伴随概 率均在10%以上，因此不能拒绝原模型中残差无自 相关的假设。可见，工具变量的选择是合适的，模 型估计结果是可取的。 就各控制变量而言，各方程中上一期的制造业