The function finds a bandwidth for a given geographically weighted regression by optimzing a selected function. For crossvalidation, this scores the root mean square prediction error for the geographically weighted regressions, choosing the bandwidth minimizing this quantity.
1 2 3 
formula 
regression model formula as in 
data 
model data frame as in 
coords 
matrix of coordinates of points representing the spatial positions of the observations 
adapt 
either TRUE: find the proportion between 0 and 1 of observations to include in weighting scheme (knearest neighbours), or FALSE — find global bandwidth 
gweight 
geographical weighting function, at present

method 
default "cv" for drop1 crossvalidation, or "aic" for AIC optimisation (depends on assumptions about AIC degrees of freedom) 
verbose 
if TRUE (default), reports the progress of search for bandwidth 
longlat 
TRUE if point coordinates are longitudelatitude decimal degrees, in which case distances are measured in kilometers; if x is a SpatialPoints object, the value is taken from the object itself 
RMSE 
default FALSE to correspond with CV scores in newer references (sum of squared CV errors), if TRUE the previous behaviour of scoring by LOO CV RMSE 
weights 
case weights used as in weighted least squares, beware of scaling issues — only used with the crossvalidation method, probably unsafe 
tol 
the desired accuracy to be passed to 
show.error.messages 
default FALSE; may be set to TRUE to see error messages if 
If the regression contains little pattern, the bandwidth will converge to the upper bound of the line search, which is the diagonal of the bounding box of the data point coordinates for “adapt=FALSE”, and 1 for “adapt=TRUE”; see the simulation block in the examples below.
returns the crossvalidation bandwidth.
Use of method="aic" results in the creation of an n by n matrix, and should not be chosen when n is large.
Roger Bivand Roger.Bivand@nhh.no
Fotheringham, A.S., Brunsdon, C., and Charlton, M.E., 2002, Geographically Weighted Regression, Chichester: Wiley; Paez A, Farber S, Wheeler D, 2011, "A simulationbased study of geographically weighted regression as a method for investigating spatially varying relationships", Environment and Planning A 43(12) 29923010; http://gwr.nuim.ie/
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26  data(columbus)
gwr.sel(crime ~ income + housing, data=columbus,
coords=cbind(columbus$x, columbus$y))
gwr.sel(crime ~ income + housing, data=columbus,
coords=cbind(columbus$x, columbus$y), gweight=gwr.bisquare)
## Not run:
data(georgia)
set.seed(1)
X0 < runif(nrow(gSRDF)*3)
X1 < matrix(sample(X0), ncol=3)
X1 < prcomp(X1, center=FALSE, scale.=FALSE)$x
gSRDF$X1 < X1[,1]
gSRDF$X2 < X1[,2]
gSRDF$X3 < X1[,3]
yrn < rnorm(nrow(gSRDF))
gSRDF$yrn < sample(yrn)
bw < gwr.sel(yrn ~ X1 + X2 + X3, data=gSRDF, method="cv", adapt=FALSE, verbose=FALSE)
bw
bw < gwr.sel(yrn ~ X1 + X2 + X3, data=gSRDF, method="aic", adapt=FALSE, verbose=FALSE)
bw
bw < gwr.sel(yrn ~ X1 + X2 + X3, data=gSRDF, method="cv", adapt=TRUE, verbose=FALSE)
bw
bw < gwr.sel(yrn ~ X1 + X2 + X3, data=gSRDF, method="aic", adapt=TRUE, verbose=FALSE)
bw
## End(Not run)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.