使用R语言中的GWmodel进行GWR模型的运算

一、安装R语言

1、R安装包下载地址：https://cran.r-project.org/
2、RStudioa安装包下载地址：https://rstudio.com/products/rstudio/

二、使用步骤

1.安装GWmodel包

打开RStudioa界面

仅安装GWmodel包载入时报错，并不能使用bw.gwr和gwr.basic函数

library(“GWmodel”)
载入需要的程辑包：spatialreg
Error: package or namespace load failed for ‘spatialreg’:
loadNamespace()里算’units’时.onLoad失败了，详细内容：
调用: udunits_init(file.path(.get_ud_xml_dir(), “udunits2.xml”))
错误: Error in function udunits_init: Can’t open installed, default, unit database
错误: 无法载入程辑包‘spatialreg’

安装GWmodel包和udunits2包

> install.packages("GWmodel")
> install.packages("udunits2")

或者Tools->install packages…

2.载入包和数据

载入GWmodel包和udunits2包

> library(udunits2)
> library(GWmodel)

加载数据

> temp<-read.table(file = "E:/temp.csv", header = TRUE, sep = ',')

GWR计算

没有构建dataframe会报错

Error in xxxxxxxxxxxxxxxxxxxx
Given regression data must be Spatial*DataFrame

构建dataframe

> map = SpatialPointsDataFrame(data = temp, coords = cbind(x = temp$X, y = temp$Y))

计算带宽

> bw1 = bw.gwr(GKDM~POPU+GDP+ONE+TWO+THREE+GDZCTZ, data=map, approach = "CV", kernel = "gaussian")

地理加权回归计算

> res<-gwr.basic(GKDM~POPU+GDP+ONE+TWO+THREE+GDZCTZ, data = map,  bw = bw1, kernel = 'gaussian')
> res

导出数据

results <-as.data.frame(res$SDF)
write.csv(x = results, file = "E:/res.csv")