Dear list,
I did GWR residual plot using the same technique as spgwr package. When I do
for Georgia data, I get the nice plot (with the state map) with the residuals,
but same code applied to my data gives residuals plot (without the state map).
My code is given below, if someone can tell me what is the problem.
library(spdep)
library(spgwr)
data<-read.csv("GWRDATA.csv",header=TRUE)
attach(data)
map = SpatialPointsDataFrame(data=data, coords=cbind(Longitude,Latitude))
bw = gwr.sel(Y~x1+x2,data=map,adapt=T)
(gwr.model = gwr(Y~x1+x2,data=map, adapt=bw))
SD.adapt.gauss <- gwr.sel(Y~x1+x2,data=map,adapt=TRUE)
res.adpt <- gwr(Y~x1+x2,data=map,adapt=SD.adapt.gauss)
brks <- c(-0.25, 0, 0.01, 0.025, 0.075)
cols <- grey(5:2/6)
res.adpt$SDF$ols.e <- residuals(lm(yield~temp+precpt,data=map))
spplot(res.adpt$SDF, c("ols.e", "gwr.e"),main="Residuals from OLS and GWR")
Milan
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
