Skip to content

KNN spatial weights matrix: asymmetric vs symmetric for spatial autoregressive models

3 messages · Roger Bivand, Dustin Duncan

#
On Thu, 16 Dec 2010, Dustin Duncan wrote:

            
The approximation to the Jacobian used in GeoDa requires that W be 
symmetric (or similar to symmetric, Ord, 1975). The MC approximation used 
in the Matlab toolbox does not. In spdep, you can use method="eigen" for 
smaller data sets, "LU" for exact Jacobian larger data sets, or "MC", with 
intrinsically asymmetric W. However, note that the differences between a 
symmetricised KNN W and the original one may not be large, as in SAR 
(simultaneous rather than CAR conditional) SSE terms, you do get W'W 
appearing, though its influence will vary. Try checking to see whether 
results vary enough to change your inferences:

library(spdep)
example(columbus)
knn4 <- knn2nb(knearneigh(coordinates(columbus), k=6))
knn4
knn4s <- make.sym.nb(knn4)
lw <- nb2listw(knn4)
lws <- nb2listw(knn4s)
summary(lagsarlm(CRIME ~ HOVAL + INC, data=columbus, listw=lw))
summary(lagsarlm(CRIME ~ HOVAL + INC, data=columbus, listw=lws))

So if you have to make the KNN weights symmetric to use - say - GeoDa, it 
may not matter that much in practice.

Hope this helps,

Roger