package = "spData")) 模型形式 form <- Z ~ PEXPOSURE + PCTAGE65P + PCTOWNHOME 带宽 带宽(bandwidth)确定了局部的范围,该包的gwr.sel()函数提供了两种确定带宽的方法:交叉验证法和AIC信息准则法。语法结构如下: gwr.sel(formula, data = list(), coords, adapt = FALSE, gweight = gwr.Gauss, me...
formula:regression model formula as in lm,表示R语言里面进行回归的方程式 data:model data frame as in lm, or may be a SpatialPointsDataFrame or SpatialPolygonsDataFrame object as defined in package sp 指定数据 coords: matrix of coordinates of points representing the spatial positions of the observati...
复制 library(rgdal)NY8<-readOGR(system.file("shapes/NY8_utm18.shp",package="spData"))form<-Z~PEXPOSURE+PCTAGE65P+PCTOWNHOMEbw<-gwr.sel(formula=form,data=NY8,gweight=gwr.Gauss,method="cv")model<-gwr(formula=form,data=NY8,bandwidth=bw,gweight=gwr.Gauss)model ## Call:##gwr(formul...
为了使用spgwr包将GWR的模型参数应用于更精细的空间比例:1.以粗略比例计算GWR 1.使用参数fit.points、predictions和fittedGWRobject再次应用步骤1。代码:
2、导入数据data(columbus,package="spData") 3、查看数据View(columbus) 进行回归,代码为: 4、进行传统OLS回归col.lm<-lm(CRIME~INC+HOVAL,data=columbus)summary(col.lm) 结果为: 2、带宽选择 ---# 计量经济学服务中心#5、带宽选择col.bw <- gwr.sel(CRIME ~ INC + HOVAL, data=columbus,coords=cbin...
2、导入数据data(columbus,package="spData") 3、查看数据View(columbus) 进行回归,代码为: 4、进行传统OLS回归col.lm<-lm(CRIME~INC+HOVAL,data=columbus)summary(col.lm) 结果为: 2、带宽选择 ---# 计量经济学服务中心#5、带宽选择col.bw <- gwr.sel(CRIME ~ INC + HOVAL, data=columbus,coords=cbin...