similar to: matrix to gal object

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 everybody, I would like to use the SpDep-package (especially the Local Moran index analysis and the Getis-Ord statistics) in R for analysing my data. However, I don't have x-y coordinates, but my data is in a distance matrix format. Is it possible to use the SpDep package with predefined distances as well instead of letting the program determining the distance itself by the function

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

Greetings: I am trying to work with spdep (everything is "brand new" downloaded this morning). OS = Windows 2000 (also up to date). The code I am using follows: #example gal.county=read.geoda("lnpilnd.GAL", row.names=NULL, skip=0) summary.nb(gal.county) Error in summary.nb(gal.county) : Not a neighbours list #end The gal file works just fine in GeoDa (also up to date).

Dear All, I am a beginner creating R packages. I followed the Leisch (2009) tutorial and the document ?Writing R Extensions? to write an example. I installed R 2.12.2 (I also tried R2.13.2), the last version of Rtools and the recommended packages in a PC with Windows 7 Home Premium. I can run R CMD INSTALL linmod in the command prompt and the R CMD check linmod. The following outputs are

Hi all, The survivals functions can be tested by the Log-rank test and others, for example the Gehan-Breslow. The graham breslow work with the alpha values. But I don't know how is the Gehan-Breslow test with R. Somebody know a type function?.. or other suggestions? Any help will be really appreciated Jos? Bustos Marine Biologist Master Apllied Stat Program University of Concepci?n

Hi everyone, I'm studying the manual name: Analysis of Epidemiological Data Using R and Epicalc, maked by: Virasakdi Chongsuvivatwong and Edward McNeil. And I can't find the data base that they use in some examples, this are the names: Chapter7.Rdata,Chapter8.Rdata,Chapter9.Rdata Somebody can tell me how can I have this files? Thk! Jos? O__ ---- Jos? Bustos M. c/ /'_ ---

Hi everyone, We are looking for some data sets working with relative risk mortality. so, someone know where can I find the data.mgus dataset and the data.mgus? Using 1384 records from Minnesota. Thank! O__ ---- Jos? Bustos M. c/ /'_ --- Master Apllied Stat Program (*) \(*) -- University of Concepci?n

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,

Hi all, I am attempting to create a weights object and perform a Moran I test as well. I have a very large spatial weights matrix (roughly 22,000x22,000) that was created in Excel and read into R, and I'm now trying to implement: library(spdep) SW=mat2listw(matrix) I am getting the following error: Error in if (any(x<0)) stop ("values in x cannot be negative"): missing value

Hello R Community, I've been run the following codes. However, I've been getting an unusual segfault that I'm unable to trace its origin. Please give me a light to decipher the "caught segfault" Thanks for you attention. Bernardo. > options(STERM='iESS', editor='emacsclient') > rm(list = ls()) > > source("fgenIGLD.R") #RNG

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

Howdy R gurus ! I am newbie to R I use R 2.0.1 in Windows XP. When I run R I got the follwoing memory error. My physical memory size is 3 Gb. My R got the memory problem when it reached to about 2 Gb. Thanks in advance, > library(spdep) > sfr.lagsarlm <- lagsarlm(sfr.data$Bldgsqft ~ sfr.data$Ncounty + sfr.data$Nugb + sfr.data$Ngroup, data=sfr.data, listw=sfr.listw,

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

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

Dear all, I would like to implement a spatial quantile regression using instrumental variable estimation (according to Su and Yang (2007), Instrumental variable quantile estimation of spatial autoregressive models, SMU economics & statistis working paper series, 2007, 05-2007, p.35 ). I am applying the hedonic pricing method on land transactions in Luxembourg. My original data set contains

Hello, I'm new to R and to Spatial Analysis and got a problem trying to create a Spatial Weights Matrix. *I us the following code to create the Neighbourslist:* >library(maptools) >library(spdep) >library(rgdal) >location_County<- readShapePoly("....") >proj4string(location_County)<- CRS("+proj=longlat ellps=WGS84") >location_nbq<-

dear R users, I have troubles performing Moran.I test as suggested on http://www.ats.ucla.edu/stat/r/faq/morans_i.htm my spatial data are longitude and lattitide of communities. The calculation of the inverse distance matrix according to the homepage (using my data) datAL <- read.csv2("C:\\Konvergenz AL.csv", header=T) ALdist <- as.matrix(dist(cbind(datAL$L?nge,

Hello, Is there any plan to use the client identifier field as part of the PXE search order ? This will be helpful when using PXELINUX on Infiniband network which do not have 6 byte Ethernet MAC address. Thanks, _________________________________________________ Nir Gal | +972-9-9717685 (o) | +972-54-499-6633 (m) Manager, Customer Support Voltaire - The Grid Backbone www.voltaire.com

Hello, I'm trying to run a lagsarlm (maximum likelihood estimation of a spatial lag model) in the spdep library ; but R gives following error message: Error in solve.default(inf, tol = tol.solve) : system is computationally singular: reciprocal condition number = 4.20137e-12 I get the same message when I try to run de lagsarlm with a bigger data set (4333 regions). The command