Hi,

I am trying to calculate Moran's I test for the residuals for a
regression equation, but I have trouble converting my  coordinates into
nb format.

I have used the dnearneigh() funtion now with an arbitrarily high upper
distance to make it include all plots. However, when I do the
lm.morantest() I get a Moran's I value which is the same as the expected
value and a P-value of 1.

I have used the following code:
>attach(lowland)
>coords <- as.matrix(cbind(long, $lat))
>coord.nb <- dnearneigh(coords, 0, 10000 longlat=TRUE)
>coord.list <- nb2listw(coord.nb, style="W")
>lianasp.lm <- lm(lianasprich ~ log(averdist) + dsl + lianadens +
wooddens)
>lm.morantest(lianasp.lm, coord.list, alternative="two.sided")

And get the following result

        Global Moran's I for regression residuals

data:  
model: lm(formula = lianasprich ~ log(averdist) + dsl + lianadens +
wooddens)
weights: coord.list
 
Moran I statistic standard deviate = 0, p-value = 1
alternative hypothesis: two.sided 
sample estimates:
Observed Moran's I        Expectation           Variance 
     -1.754386e-02      -1.754386e-02       2.125706e-16 

Can anyone tell me where I went wrong or if there is another way to
generate a nb list? 

Many thanks,
Geertje

~~~~
Geertje van der Heijden
PhD student
Tropical Ecology
School of Geography
University of Leeds
Leeds LS2 9JT

Tel: (+44)(0)113 3433345 
Email: g.m.f.vanderheijden04@leeds.ac.uk



	[[alternative HTML version deleted]]