Apologies for cross-posting, its unclear if this is a spatial or a mixed models question. Maybe its both?
I am interested in extracting the residuals from level 2 of a mixed model (created using lme()), building a spatial weights matrix, and then testing for spatial autocorrelation using Moran's I. While I have found methods of incorporating spatial effects into a mixed model using corStruct, I am interested in first evaluating *if* that is an appropriate model for a given dataset by examining the level 2 residuals' spatial patterning - or lack thereof. Which slot contains the residuals for level 2 and how are they ordered? Best, Dexter [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo