加权一阶局域法matlab程序_地理加权回归模型介绍

写这篇文章的目的：

最近看了一篇删失数据中的地理加权COX模型，COX模型中考虑加入患者地理位置信息。该文章主要是理论推导类问题文章，文章中算法对应R包还未发布。发布后，在模型中添加患者地理信息也是个不错的点。

在地学空间分析中，n组观测数据通常是在n个不同地理位置上获取的样本数据，全局空间回归模型就是假定回归参数与样本数据的地理位置无关，或者说在整个空间研究区域内保持稳定一致，那么在n个不同地理位置上获取的样本数据，就等同于在同一地理位置上获取的n个样本数据，其回归模型与最小二乘法回归模型相同，采用最小二乘估计得到的回归参数户既是该点的最优无偏估计，也是研究区域内所有点上的最优无偏估计。而在实际问题研究中我们经常发现回归参数在不同地理位置上往往表现为不同，也就是说回归参数随地理位置变化，这时如果仍然采用全局空间回归模型，得到的回归参数估计将是回归参数在整个研究区域内的平均值，不能反映回归参数的真实空间特征。为了解决这一问题，国外有些学者提出了空间变参数回归模型(SpatiallyVarying-Coeffi Cient Regression Model)(Fosterand Gorr，1986；Gorrand Olligschlaeger，1994)，将数据的空间结构嵌入回归模型中，使回归参数变成观测点地理位置的函数。Fortheringham等(Brunsdonetal，1996；Fortheringhamet al，1997；Brunsdonet al，1998)在空间变系数回归模型基础上利用局部光滑思想，提出了地理加权回归模型(GeographieallyWeighted Regression Model-GWR)。

地理加权回归模型(GWR)是对普通线性回归模型(OLR)的扩展，将样点数据的地理位置嵌入到回归参数之中，即：

     式中：(u_i，v_i)为第i个样点的坐标(如经纬度)；β_k(u_i，v_i)是第i个样点的第k个回归参数 ；E_i是第i个样点的随机误差。为了表述方便，我们将上式简写为：

若，则地理加权回归模型(GWR)就退变为普通线性回归模型(OLR)。

Fotheringhamet al依据“接近位置i的观察数据比那些离i位置远一些的数据对的估计有更多的影响”(Fotheringhamet al，1996)的思想，利用加权最小二乘法来估计参数，得

其中：