地理加权回归

gwr.sel <- function(formula, data = list(), coords, adapt=FALSE,  
                    gweight=gwr.Gauss, method="cv", verbose=TRUE, longlat=NULL, 
                    RMSE=FALSE, weights, tol=.Machine$double.eps^0.25, 
                    show.error.messages=FALSE) { 
  if (!is.logical(adapt)) stop("adapt must be logical") 
  if (is(data, "Spatial")) { 
    if (!missing(coords)) 
      warning("data is Spatial* object, ignoring coords argument") 
    coords <- coordinates(data) 
    if (is.null(longlat) || !is.logical(longlat)) { 
      if (!is.na(is.projected(data)) && !is.projected(data)) { 
        longlat <- TRUE 
      } else { 
        longlat <- FALSE 
      } 
    } 
    data <- as(data, "data.frame") 
  } 
  if (is.null(longlat) || !is.logical(longlat)) longlat <- FALSE 
  if (missing(coords)) 
    stop("Observation coordinates have to be given") 
  mf <- match.call(expand.dots = FALSE) 
  m <- match(c("formula", "data", "weights"), names(mf), 0) 
  mf <- mf[c(1, m)] 
  mf$drop.unused.levels <- TRUE 
  mf[[1]] <- as.name("model.frame") 
  mf <- eval(mf, parent.frame()) 
  mt <- attr(mf, "terms") 
  dp.n <- length(model.extract(mf, "response")) 
  #	mt <- terms(formula, data = data) 
  #	mf <- lm(formula, data, method="model.frame", na.action=na.fail) 
  #	dist2 <- (as.matrix(dist(coords)))^2 
  weights <- as.vector(model.extract(mf, "weights")) 
  # set up default weights 
  if (!is.null(weights) && !is.numeric(weights))  
    stop("'weights' must be a numeric vector") 
  if (is.null(weights)) weights <- rep(as.numeric(1), dp.n) 
  if (any(is.na(weights))) stop("NAs in weights") 
  if (any(weights < 0)) stop("negative weights") 
  y <- model.extract(mf, "response") 
  x <- model.matrix(mt, mf) 
  #	if (