Skip to content
Prev 30652 / 398506 Next

How to do the significant test on Local Moran's I

On Mon, 14 Apr 2003, Danlin Yu wrote:

            
The ideas are in the documentation and references of the functions you
refer to in the spdep package: localmoran(), localG(), and
localmoran.sad(). You need to recall that doing lots of local
"significance" tests on the same data means that you have to apply
corrections, as in p.adjust(), to any p-values you might compute. If you
are just testing a single relationship (values of x in Rhode Island are
correlated with values of x in its contiguous neighbours), you can do that
in the standard way, but you cannot extend this to gat a map of p-values -
they will be very misleading, as the references point out - Ord, J. K. and
Getis, A. 1995 Local spatial autocorrelation statistics: distributional
issues and an application. Geographical Analysis, 27, 286-306 - have a 
table of corrections. 

The functions in the package let you compute the pieces you need to do the
test, but do not provide any p-values, because the function cannot know
how many tests you are doing on the same data - you have to do that. That
is also why localmoran.sad() returns a list of "htest" objects, to point
up the fact that you should decide yourself what you are trying to test.

Please also be aware that by modifying the boundaries of the aggregations 
you may be analysing, you can often choose the test results you might like 
(the Modifiable Areal Unit Problem), so your "significance" tests may not 
actually be very informative.

Roger