使用GWmodel进行GWR模型相关运算

最新推荐文章于 2024-07-21 11:05:11 发布

serenahxr

最新推荐文章于 2024-07-21 11:05:11 发布

阅读量3.4k

点赞数 3

分类专栏：空间分析

本文链接：https://blog.csdn.net/tjhxrrr/article/details/105105277

版权

空间分析专栏收录该内容

6 篇文章 7 订阅

订阅专栏

未完但不待续

安装GWmodel包

不要试图在anaconda的环境的RStudio中安装！

引入GWmodel包

读取数据

guangzhouershoufang<-read.table(file="/Users//广州/4-数据分析/广州二手房Sample2CSV.csv", header = TRUE, sep=',')

构建dataframe

map = SpatialPointsDataFrame(data=guangzhouershoufang, coords=cbind(x=guangzhouershoufang$lon, y=guangzhouershoufang$lat))

计算带宽

bw=bw.gwr(price~area+room+hall+toilet+direction+floor+age+volumeRatio+plantRatio+location+subway+pSchool+mSchool+university+hospital, data=map, approach="CV",kernel="gaussian")

个人建议用spgwr里的gwr.sel函数算带宽。用GWmodel里的这个bw.gwr函数算出来的带宽在下一步构建GWR模型时有可能失败（报错：matrix seems singular），我猜测可能的原因是这个函数在计算带宽时，在迭代寻找最优带宽的过程中，一旦有的带宽对应cv值或aic值算不出来，它就会使用一些近似值，然后为下一步埋下隐患。

构建GWR模型

 gzesf.res<-gwr.basic(price~area+room+hall+toilet+direction+floor+age+volumeRatio+plantRatio+location+subway+pSchool+mSchool+university+hospital, data=map,  bw=bw, kernel='gaussian')