空间地理加权回归stata_如何用stata做地理加权回归（GWR）

地理加权回归是空间计量经济学中处理空间异质性的重要模型，现有的stata目前可以支持GWR分析，该命令包为gwr，其语法格式为：

gwr depvar [varlist] [if exp] [in range] , east(varname) north(varname) [options]

其中，east存放个体的经度坐标或X轴位置坐标的变量，north存放个体的纬度坐标或Y轴位置坐标的变量.

options里选项有：

saving(filename) dots reps(#) double eform family(familyname)

link(linkname) [ln]offset(varname) test replace noconstant

nolog scale(x2|dev|#) disp(#) iterate(#) init(varname)

outfile(filename) comma wide bandwidth(#) mcsave(filename)

sample(#)

其中，常用到如下：

检验专用：

保存文件专用：

模拟次数：

其他option可以help gwr

以下具体以co13_d90等数据为例。

shp2dta using co13_d90, data(georgia_data) coor(georgia_coor)

replace

use georgia, clear

经典回归命令:

regress PctBach TotPop90 PctRural PctEld PctFB PctPov PctBlack

结果如下:

而地理加权回归命令则为：

gwr PctBach TotPop90 PctRural PctEld PctFB PctPov PctBlack, north(Latitud) east(Longitud) reps(1000) saving(georgia_gwr.dta) replace

需要注意的是，在reg回归下，每个变量只有一个参数估计值，而在GWR模型中则每一地区对应一个系数，与每个地区有关的系数beta_i，已经放入georgia_gwr的文件中，因此，为了对其进行描述，可以使用sum命令，也可以结合spmap命令对不同系数在不同地区的差异性情况进行描述。reps(1000)导致运行时间会非常长，测试时可以减少次数。

结果如下：