No subject

moran.test(my_shape$variable_of_interest, listw=nb2listw(my_nb, style="W", 
zero.policy=TRUE), zero.policy=TRUE)

roughly, depending on whether you have cells that have zero neighbours. 
But we can get back to this once the shapefile is on board.

Roger

> 
> > How had you considered handling the unconnected islands? It looks as
> > though most of the map has low values, and only a few higher values - this
> > is effectively saying that the variable is predicted by its mean, and
> > since most of the map is below the mean, and the rest above, the result is
> > a foregone conclusion, that is the data do not look mean-stationary. So I
> > don't think Moran's I will tell you much you don't know, but I could be
> > wrong.
> 
> I don't know about the shape files, I guess that they could either be 
> neighbours to the nearest cells or then as separate.
> You are right in this case, but not all cases are so easily interpreted. 
> I'd like to learn how to do this on this material (I have done the 
> analysis in ArcInfo, but would like to know how to do it with free tools.
> 
> I think that I should also reveal that I am a programmer, and I have very 
> little geoinformatics or geography education, so please bear with me. At 
> the moment I need this for a course in geoinformatics, and I'm trying to 
> learn more.
> 
> --Wolf
> 
> 

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no