Moran's I for gwr objects

Description

The function returns Leung et al. (2000) three moment approximation for Moran's I, for a gwr object calculated with argument hatmatrix=TRUE. This implementation should not be regarded as authoritative, as it involves assumptions about implied methods and about estimated degrees of freedom.

Usage

1
gwr.morantest(x, lw, zero.policy = FALSE)

Arguments

x

a gwr object returned by gwr() with argument hatmatrix=TRUE

lw

a listw object created for example by nb2listw in the spdep package

zero.policy

if TRUE assign zero to the lagged value of zones without neighbours, if FALSE (default) assign NA

Value

a “htest” object with the results of testing the GWR residuals

Author(s)

Roger Bivand

References

Leung Y, Mei C-L, Zhang W-X 2000 Testing for spatial autocorrelation among the residuals of the geographically weighted regression, Environment and Planning A, 32, 871-890.

Examples

1
2
3
4
5
6
7
if (suppressWarnings(require(spdep)) ){
  data(columbus)
  bw <- gwr.sel(CRIME ~ INC + HOVAL, data=columbus, coords=coords)
  col0 <- gwr(CRIME ~ INC + HOVAL, data=columbus, coords=coords,
    bandwidth=bw, hatmatrix=TRUE)
  gwr.morantest(col0, nb2listw(col.gal.nb))
}