地理加权回归-空间分析.docx

地理加权回归-空间分析

地理加权回归

空间回归模型中的回归系数不随空间位置而变化，因此空间回归模型是全局模型。但是由于空间异质性和空间非平稳性，不同空间子区域中自变量和因变量的关系很可能不同，因此就需要处理空间异质性的局部空间回归方法，因此就有了地理加权回归模型（GeographicallyWeightedRegression,GWR）的提出。地理加权回归同时考虑了空间的相关性与异质性。

地理加权回归模型一般形式如下：

；

其中，与是因变量与自变量在处的观测值，为第个采样点的坐标，是第个采样点上的第个回归参数，点被称为回归点。，当时，。简便记为：

地理加权回归会得出个回归方程，对应每个回归点，都会有一个回归方程。若，则地理加权回归模型退化为普通线性回归模型。

模型回归参数需要通过局部加权最小二乘进行估计。假设在位置的权重为，（），那么位置的参数估计为使达到最小的值。

上式中的为权函数，反映其他观测点的样本对回归点的影响。权函数值越大，影响越大。该值通常由对应两点之间的距离决定。常用的权函数有：

（1）Gauss函数：

（2）bi-square函数：

以上两式中的为点到点的距离。被称为带宽，是需要人工选择的参数。这两个函数中，距离越大，函数值都越小。这说明选择这两个函数时，都假设观测点越远，影响越小。

地理加权回归模型将数据的空间位置嵌入到回归参数中，因此考虑了空间的异质性。同时，由于不同位置的观测点对回归参数的影响大小不同（通常离回归点越近，影响越大），因此该模型也考虑了空间相关性。地理加权回归的使用，也当同时以空间相关性与异质性为前提。如果没有空间相关性，那么该模型就缺乏合理性。如果没有空间异质性，那么该模型就缺乏必要性。

【案例STYLEREF1\s7SEQ案例\*ARABIC\s12】

以2011年北京、天津、唐山各县（区）疾病发病率和同期各县（区）的人口密度、人均GDP、年平均风速、光照强度、相对湿度、年降水量等数据为例，研究社会经济和气象因素对该疾病发病率在各地区影响的不稳定性。

本节所用的数据为2011年某地区某疾病的发病率（C:\Example\Data\7.3WGR\地区汇总.csv），命名为“地区汇总.csv”。此处只展示前10条

表STYLEREF1\s7SEQ表\*ARABIC\s132011年某地区的某疾病发病情况

code

precp

relHum

sunShn

wndspd

PopuDen

PerGdp

rate

110101

58.06842

49.82301

206.0755

2.195716

15318.95

0.002438

347.1859

110102

62.59835

52.05048

202.8901

2.169381

14956.04

0.003675

447.7185

110105

57.65938

50.06015

205.5302

2.189571

8528.913

0.00096

1139.069

110106

55.73981

50.20182

205.5477

2.215026

7698.981

0.000419

1759.943

110107

52.39305

50.77886

205.6759

2.255836

8249.254

0.000541

1561.691

110108

52.13075

51.01494

207.5023

2.255693

8524.685

0.001021

1010.973

110109

44.66608

53.70504

218.3655

2.473628

224.6719

0.00038

758.6155

110111

46.62295

52.79277

208.7774

2.333494

532.4592

0.000457

2082.544

110112

55.23023

51.21378

201.2761

2.186171

1454.807

0.00036

2406.252

110113

50.22722

55.34219

195.5819

2.014273

967.6505

0.001211

1643.1

表STYLEREF1\s7SEQ表\*ARABIC\s14各变量的含义

变量名称

变量含义

单位

code

地理编码

无

precp

降水量

毫米

relHum

相对湿度

%

sunShn

日照强度

瓦/平方米

wndspd

风速

千米/时

popuDen

人口密度

千人/平方公里

PerGdp

人均国内生产总值

千元

rate

发病率

此外，还需要的数据是包含该地区所有区县的地图文件（.dbf文件和.shp文件从光盘中获取，C:\Exam