R语言地理加权回归、主成份分析、判别分析等空间异质性数据分析

在自然和社会科学领域有大量与地理或空间有关的数据,这一类数据一般具有严重的空间异质性,而通常的统计学方法并不能处理空间异质性,因而对此类型的数据无能为力。以地理加权回归为基础的一系列方法:经典地理加权回归,半参数地理加权回归、多尺度地理加权回归、地理加权主成份分析、地理加权判别分析是处理这类数据的有效模型。本次培训从局部加权回归开始,详细讲述了基于R语言的空间异质性数据分析方法。

图片

专题一:地理加权回归下的描述性统计学
1.R语言操作简单回顾
2.局部加权的基本原理
3.带宽与核函数选择
4.局部加权的均值,标准差和相关系数
5.分位数及基于分位数的稳健估计

图片

图片

专题二:地理加权主成分分析
1.普通的主成分分析,因子载荷与因子得分分析
2.主成分个数的选择,碎石图
3.地理加权的主成分分析
4.主成分的空间载荷
5.空间主导因子分析

图片

专题三:地理加权回归
1.线性回归:高斯-马尔科夫假设
2.地理加权回归:基本方法与稳健方法,异常值的检验
3.带宽选择:修正的赤池信息法
4.系数检验:F1,F2,F3检验
5.空间稳定性检验:蒙特卡洛方法
6.共线性与变量选择:地理加权回归中的岭回归与Lasso回归
7.时空地理加权回归:GTWR
8.QGIS中的地理加权回归

图片

专题四:高级回归与回归之外
1.多尺度地理加权回归:可变带宽的选择
2.异方差模型 
3.广义地理加权回归:链接函数,泊松回归与二项式回归
4.空间权重矩阵与半参数地理加权回归
5.分位数回归与地理加权分位数回归
6.判别分析与地理加权判别分析

图片

 

 

  • 11
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
R语言中可以使用gwmodel包进行地理加权回归分析。 首先安装gwmodel包: ```r install.packages("gwmodel") ``` 然后加载包: ```r library(gwmodel) ``` 使用数据集meuse作为例子,该数据集包含了荷兰一个小地区的铜、铅等重金属元素浓度以及地理位置信息。首先加载数据集: ```r data(meuse) ``` 进行地理加权回归分析需要用到的变量包括自变量和因变量,以及空间变量(例如坐标)。在这个例子中,因变量为铅(lead)的浓度,自变量为距离河流的距离(dist.m)和土壤pH值(pH),空间变量为坐标(x、y)。我们先将数据集按照空间变量进行排序: ```r meuse <- meuse[order(meuse$x, meuse$y), ] ``` 然后将空间变量和自变量提取出来: ```r coords <- data.frame(x = meuse$x, y = meuse$y) xdata <- data.frame(dist.m = meuse$dist.m, pH = meuse$pH) ydata <- meuse$lead ``` 接下来使用gwr函数进行地理加权回归分析: ```r fit <- gwr(coords = coords, x = xdata, y = ydata) ``` 默认情况下,gwr函数使用的是高斯核函数,可以通过bandwidth参数来调整核函数的带宽。带宽越小,核函数的权重越集中,模型的拟合效果会更加精细,但是模型的复杂度也会增加。带宽越大,核函数的权重分布越平均,模型的拟合效果会更加平滑,但是模型的泛化能力会降低。我们可以通过交叉验证来选择最优的带宽值: ```r cvfit <- gwr.cv(coords = coords, x = xdata, y = ydata) cvfit$bw ``` 最优的带宽值会保存在cvfit$bw中。 除了高斯核函数之外,gwr函数还支持其他的核函数,例如指数核函数和双三次核函数。可以通过kernel参数来指定核函数的类型。例如,使用指数核函数: ```r fit <- gwr(coords = coords, x = xdata, y = ydata, kernel = "exponential") ``` 使用双三次核函数: ```r fit <- gwr(coords = coords, x = xdata, y = ydata, kernel = "bisquare") ``` gwr函数的输出结果包括每个样本的系数估计值、t值、p值和R2值等等。我们可以将结果可视化,例如绘制系数估计值的热力图: ```r library(ggplot2) library(RColorBrewer) ggplot(fit$beta, aes(x = coords.x1, y = coords.x2, fill = est)) + geom_raster() + scale_fill_distiller(palette = "RdBu", type = "div") + coord_fixed() + theme_void() + theme(legend.position = "bottom") ``` 这里使用ggplot2和RColorBrewer库来制作热力图。结果如下图所示: ![image](https://user-images.githubusercontent.com/26807215/128842726-6f09b6d9-6de1-4011-8356-2e6d47d16e84.png) 可以看到,距离河流的距离对铅的浓度有显著的正向影响,而土壤pH值对铅的浓度则有显著的负向影响。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值