Displaying 1 result from an estimated 1 matches for "hatmatrix".
Did you mean:
fatmatrix
2008 Oct 30
1
A question about pairs()
...h together.
For example:
data(columbus)
# compute gwr bandwidth
col.bw <- gwr.sel(crime ~ income + housing, data=columbus,
coords=cbind(columbus$x, columbus$y))
# estimate gwr
col.gauss <- gwr(crime ~ income + housing, data=columbus,
coords=cbind(columbus$x, columbus$y), bandwidth=col.bw, hatmatrix=TRUE)
# summarize gwr results
col.gauss
# compute local correlation matrix among gwr coefs
col.cor <- gw.cov(col.gauss$SDF,vars=2:4,bw=col.bw,cor=T,longlat=F)
# extract gwr coefficients and local correlations
col.gauss.mat <- as(col.gauss$SDF, "data.frame")[,2:4]
col.cor.mat <- a...