[R-sig-Geo] globalG.test - x vector with different lenght as the neighbours list

Francesco Perugini francescoperugini07 at gmail.com
Tue Feb 2 16:15:17 CET 2016

Dear all,
I've a routine to create spatial weight matrix for different k and then to
implement, for each of the created spatial weighted matrix, the Global G
test for autocorrelation.
The point is that I'll try to do this test on multiple variables but then
lenghts are different

> length(dlwknn.B$neighbours)
[1] 49
> length(var)
[1] 1

Any help on how to code this?
Thanks a lot.

This is the routine:

coord <- coordinates(columbus)

z <- c(1,2,3,4,5,6,7,8,9)
neighbors.knn <- list()

for (val in z)
  neighbors.knn <- c(neighbors.knn, list(knn2nb(knearneigh(coord, val,
longlat=F), sym=F)))


plot(neighbors.knn[[1]], coord)

plot(neighbors.knn[[2]], coord)

ndists <- length(z)
ZG <- vector(mode="list", length=ndists)
names(ZG) <- as.character(z)

f <- c(CRIME, INC)
for (val in z) {
  dlwknn.B <- nb2listw(neighbors.knn[[val]], style="B", zero.policy=TRUE)
  temp <- list()

  for (var in f)
    temp <- c(temp, list(globalG.test(var, dlwknn.B, zero.policy=F)))
  ZG[[val]] <- temp

t(sapply(ZG, function(var) c(var$estimate[1], var$statistic,

	[[alternative HTML version deleted]]

More information about the R-sig-Geo mailing list