similar to: how to store estimates results as scalars of a matrix?

Hello, I am interested in using Morans I for different time intervals to detect disease clusters. Ultimately I would like to use CUSUM - or similar monitoring statistic to monitor the results of Morans I - similar to the work by Rogerson (2005) Spatial Surveillance and Cummulative Sum Methods in Spatial and Syndromic Surveillance for Public Health. Thus far - thanks to the list I have

I want to do "moran's I test" in R language. I try to use "gearymoran" in Package "ade4","moran" in Package "spdep", and Moran.I in Package "ape". But I do not know how to do it because data format is different. My data: x y dbh 111.03 10.7 7 118.11 0.28 1.2 165.36 0.36 8.4

Hello List, I'm a very new user to the R system. I'm only beginning to learn the basics, but so far I've been able to do little more than try a few examples, and of course begin reading the documentation. My primary motivation for exploring R is the availability of tools like the 'spdep' package for calculating spatial statistics such as Geary's C and Moran's

Hi R helpers, Is there a function (package?) in R available which tests "spatial autocorrelation" between points (e.g. vector layer of weather stations)? (e.g. Moran's I...) Via the archives we found out that there is a package 'spdep' which uses grid data for testing spatial autocorrelation. Thanks a lot, Jan

hi,friends, we all know that moran's I index and Geary'C index can be used to test spatial autocorrelation in both the area data and point data, but i only can find something on how to calculate on the data of area, and can't find the methods to perform it on the point data, could anybody give me some information, thanks in advance! -- Kind Regards, Zhi Jie,Zhang ,PHD Department of

Hi, I'm have some data on a grid (specifically high throughput assay plates) and am interested in evaluating measures of spatial autocorrelation to flag plates for corrections. I have been using moran.test and geary.test from the spdep package. My approach is as follows: ## plate is a matrix of data coords <- expand.grid(1:32, 1:48) x <- as.numeric(plate) moran.test(x,

Hi group, I'd like to do spatial analysis of my data using the spdep package. It appears that a file of class nb is necessary, but I do not find what that should look like and if there is a function that creates such file for me. How can I create a nb-object of my data points? Thanks, René Eschen. [[alternative HTML version deleted]]

hello, I want to do a test for spatial correlation. I tried it with geary.test() but I don't understand the required input. x= a numeric vector the same length as the neighbours list in listw (my sampled data, I assume) listw= a listw object created for example by nb2listw (well when I check nb2listw() I get to "neighbours - an object of class nb" - but I couldn't figure

Dear list, (also sent to Roger Bivand, but perhaps somebody of you can help me also) I am trying to use package spdep for fitting an SAR model with errorsarlm. However, I am not sure how to make a valid nb object out of my neighborhood. As far as I have seen, there is no documentation for nb.object. I have done the following: class(pschmid$nb) <- "nb" # pschmid is a prab object as

Hi, I have a question on testing spatial autocorrelation on raster data including NA-values. In particular, I like to calculate Moran?s I and Geary?s C indices by using inverse distance weighting matrices. Calculating Moran?s I with moran.test works fine, because the function contains the option "na.action=na.pass". Unfortunately, the function geary.test does not contain this

Dear R-users, I was wondering if someone could give me some advices on the following problem. I tried to apply moran’s test to a small dataset and couldn’t succeed, here is the error message:   mor <- moran.test(x, res2)   Avis dans moran.test(x, res2) :   Out-of-range p-value: reconsider test arguments   mor     Moran's I test under randomization data:  x  weights: res2 

Dear , This is Elaine. I am computing moran's I using moran.test for a generalized linear model (multiregression). The following contents are the results, and I cannot find the observed Moran's I mentioned as estimate in the manual. Please kindly help indicate if there is observed Moran's I did not notice or other method for calculation. Thanks Elaine data: residuals(modg1)

Dear R users, I have two data frames that consist of statistical information for most countries around the world. One dataframe consists of the latitude and longitude ("coord.csv") of each country, while the other consists of 100's of different attributes ("countryattri.csv") for each country (like, GDP, Population, etc.). The data is organized with a header and then

Hi Craig, I have cc''d your request to the mailing list. Thank you for the idea. :) Regards, Dan --- Craig Moran <craig.m.moran.ruby@gmail.com> wrote: > Hey, folks- > I saw the call for some development ideas. I would > most like to see > an admin utility that interfaces with Active > Directory. As a minor > stepping stone, it would be able to iterate (just

Hi, I am trying to do a Moran's I test on an ordinal logistic regression model. I have a simple spatial weights matrix listed below I would like to use. Y= 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 1 0 1 1 0 0 0 0 1 1 0 0 0 1 1 1 1 0 0 0 0 0 1 1 1 1 0 0 0 0 0 1 1 1 1 0 0 0 0 0 1 1 1 1 0 0 0 1 1 0 0 0 0 1 1 0 0 1 0 0 0 0 1 1 I try to run the test as follows- moran.test(order$resid,

Hello all, I'm a new user of R and I have used the spdep package to calculate Moran's I estimates. I used the dnearneigh() function and moran.test() function and obtained what I wanted, but now I want to plot my results with the p-values. I regrouped my results for the Moran's I like >m=c(m1$e,m2$e...) #m1 contain the results of the moran.test for lag1

Estimation of spatial probit models is as far as I know not implemented functions in the spdep package of R. I am looking for a scripting that handles both endogeneity and heteroscedasticity in spatial probit models. This estimator was first suggested by Case (1992). If someone could provide me with such an R-script I would be most grateful. Moreover I would be interested in obtaining script for

Apologies for my previous effort in HTML which apparently was scrubbed Dear R-users I wonder if I could get advice on the above problem I have just installed V 2.12.0 (I chose only the 32-bit version) into a new directory (C:/R) on a 64bit Windows 7 machine > sessionInfo() R version 2.12.0 (2010-10-15) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_Australia.1252

Hi Are there any more efficient ways of calculating the neighbourhood object for large datasets? I am trying to compute Moran I statistics for a very large data set (over 14,000 points). I have been using moran.test from the spdep package and everything works fine for a small data set (200 points). However, applying the same script to the whole dataset is taking days to compute (it so far has

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