similar to: How to handle NA-values in raster-based Geary´s C test?

Dear Users: I have been trying to use the geary.test() function in *R*, but am having slight difficulty understanding how I am to apply it in my context. I have 2 matrices: 1) *n x p* matrix of *n* observations with *p* measurements each. It may be noted that this matrix has a spatial dimension to it, as the *n*observations are at different geographical locations on a map. 2) *n x n* spatial

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 All, We aim to remove the spatial structure of our data using Moran Eigen Vectors and spdep package . Our data has 3694 samples and 13 variables. The computer stop working after almost 4 days of processing (we found it emitting a sharp sound and with all colors on the screen. No wories, it was restared without problem!). And we are left with nothing: no result file was produced since the

-----Original Message----- From: r-help-request at stat.math.ethz.ch [mailto:r-help-request at stat.math.ethz.ch] Sent: Tuesday, July 01, 2003 3:09 AM To: r-help at stat.math.ethz.ch Subject: R-help Digest, Vol 5, Issue 1 At 18:48 30/06/2003, Martin Wegmann wrote: >hello, > >I want to do a test for spatial correlation. >I tried it with geary.test() but I don't understand the

Hello, I have been using the package "spdep" to run spatial regressions on a data set with about 2500 observations. It has performed well up until now, but the following code resulted in an error: > load("Panel.90s.ok.R") > attach(Panel.90s.ok) > neighs<-dnearneigh(cbind(x,y),0,50000) > help(nbdists) > dists<-nbdists(neighs,cbind(x,y)) >

David Geary will be speaking at the Pikes Peak Java Developers group meeting, Tuesday, November 8. The meeting starts at 6:00pm. So if you live in the Colorado Springs area, come out and help support rails! For the detailed agenda, visit, http://www.ppjdg.org Also, if you are interested in starting a Ruby Group in Colorado Springs, please send me an email and let me know. -- Joe Clarke

Hello, i' m looking for a function to compute geary statistic for spatial correlation. Thanks. -- St?phane DRAY --------------------------------------------------------------- Biom?trie et Biologie ?volutive - Equipe "?cologie Statistique" Universite Lyon 1 - Bat 711 - 69622 Villeurbanne CEDEX - France Tel : 04 72 43 27 56 Fax : 04 78 89 27 19 04 72 43 27 57 E-mail

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 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

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 all, I would like to calculate a Moran's I statistic using the moran function in the spdep package. The problem I'm having deals with how to create the listw object. My data stems from the area of social network analysis. I have list of poeple and for each pair of them I have a measure of their relationship strength. So my dataset looks like: Jim; Bob; 0.5 This measure 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 all I am trying to run samba-3.0.28 on a solaris 10 sparc machine, as a member server for an AD 2003R2. I use nss_ldap-254. I compile samba with: GCC_EXEC_PREFIX=/opt/csw/gnu/ LDFLAGS="-L/opt/csw/lib" CPPFLAGS="-I/opt/csw/include" ./configure --prefix=/opt/samba --with-acl-support --with-quotas --with-included-popt --with-privatedir=/opt/samba/etc/samba/private

Hola Comunidad se que este grupo es para asuntos de R. Pero me asalta la duda... *Estoy comenzando un trabajo de estadística espacial sobre el delito robo, es mi trabajo de grado ;* Estoy recolectando documentación, que me sirva para el marco teórico de la tesis ,trabajaré sobre los software gvsig y crimestat (ambos libres): Hasta ahorita tengo pensado obtener: ü Estadísticos descriptivos

Hola, buen día a todos. He estado mirando cómo interpolar unos datos anuales con el método IDW, pero no consigo aún realizarlo. idw_models <- list() for (i in 1:ncol(data_y_mt)) { #Extract the current variable current_var <- names(data_y_mt)[i] #Create an IDW model for the current variable idw_model <- idw(current_var, sta_mt, grid, p = 1) idw_models[[i]] <- idw_model }

Can someone tell me what an object of class lm returned by lm means? I assumed it mean the regression model - but I'm not sure how to enter this in. I have tried y~a+b but this is not working. I have also tried saving the regression results and entering these, but again this is incorrect. This language is from the following: lm.LMtests(model, listw, zero.policy=FALSE,

Emilio, muchas gracias. Tuve que rehacer el código para entender mejor todo, pero me sale el siguiente error: *Error in as.data.frame.default(x[[i]], optional = TRUE) : cannot coerce class ?structure("RasterLayer", package = "raster")? to a data.frame* El código es el siguiente: *for (i in 1:nrows_mt) { data_it_mt <-

Hi, I need help! I am trying to iterate an iterative process to do cross vadation and store the results each time. I have a Spatial data.frame, called Tmese > str(Tmese) Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots ..@ data :'data.frame': 14 obs. of 17 variables: .. ..$ ID : int [1:14] 73 68 49 62 51 79 69 77 57 53 ... .. ..$

Dear list, I am trying to set up a propensity-weighted regression using the survey package. Most of my population is sampled with a sampling probability of one (that is, I have the full population). However, for a subset of the data I have only a 50% sample of the full population. In previous work on the data, I analyzed these data using SAS and STATA. In those packages I used a propensity weight