similar to: Geodata object border

Dear R users: I want to model mosquito count data based on landcover attributes and meteorological variables using a Poisson GLSM in the geoRglm package. I have monthly mosquito counts over more than 20 years with repeated observations from individual trap sites over time. I have used as.geodata() to successfully read my dataset into the geodata format utilized by geoR and geoRglm,

Hi, Im using geoR and I'm trying to do some predictions, based on an external trend. I'm having some problems specifying my model.control, specifically how do I define my model, and also the source of the covariate data at the prediction locations? I am assuming that the covariate data at the prediction locations should be imported to a geodata object along with the prediction location

Dear R users. I'm trying to fit a generalised linear spatial mode using the geoRglm package. To do so, I'm preparing my data (geodata) as follow: geoData9093 = as.geodata(data9093, coords.col= 17:18, data.col=15,* covar.col=16*) where covar.col is a factor variable (years in this case 90-91-92-93)). Then I run the model as follow: / model.5 = list(cov.pars=c(1,1),

Hello, I'm new to this list. Sorry if my question or parts of it already came up before. For my research in geostatistics, I am working with large sets of data in R (basically large matrices containing discrete x and y coordinates and a value for a certain parameter). These sets are obtained by kriging. The operation I'd like to perform is smoothen the output data set. I want to do it by

Hey, I have a data set of 1.20MB. I used read.table( ) to read in the data and then tried to convert it to geodata type. But it failed and I got "Error in vector("double", length) : cannot allocate vector of length 1351974000" message. All the commands I used were: rawdata<-read.table(file=" ") datag<-as.geodata(rawdata) BTW, what is the biggest origin

I've enconuntered this problem with the last cran version of geoR: > library(geoR) > day <- rep(1:2, each=5) > coords <- matrix(rep(runif(10),2), 10, 2) > data <- rnorm(10) > data[1] <- NA > as.geodata(cbind(coords, data, day), realisations=4) as.geodata: 1 points removed due to NA in the data Errore in as.geodata(cbind(coords, data, day), realisations = 4) :

Hi, I'm a new R user. My goal is to do a variogram using geoR. I started by trying to do the example in the geoR Illustrative Session using my own data. I am able to read in my Ascii data using: D <- matrix(scan("file.dat", n=530*3), 530,3, byrow=TRUE). Then I use: as.geodata(D, coords.col=1:2, data.col=3) to make the object D geodata. I check the descriptive statistics,

Hello, I install R under Mandrake Linux 8.2. R itself work fine, but I had an error to install packages, i.e. fields, geoR and geoRglm : $: R CMD INSTALL geoRglm_0.4-3.tar.gz Installing *source* package `geoRglm' ... libs gcc-3.0.1 -I/usr/lib/R/include -I/usr/local/include -mieee-fp -D__NO_MATH_INLINES -fPIC -O3 -fomit-frame-pointer -pipe -mcpu=pentiumpro -march=i586 -fno-fast-math

Hello, I install R under Mandrake Linux 8.2. R itself work fine, but I had an error to install packages, i.e. fields, geoR and geoRglm : $: R CMD INSTALL geoRglm_0.4-3.tar.gz Installing *source* package `geoRglm' ... libs gcc-3.0.1 -I/usr/lib/R/include -I/usr/local/include -mieee-fp -D__NO_MATH_INLINES -fPIC -O3 -fomit-frame-pointer -pipe -mcpu=pentiumpro -march=i586 -fno-fast-math

Dear, we have a problem with geoR. We try to read an ASCII table (x,y,z) with 40000 lines. With read.geodata we get the error: Error in vector("double", length) : cannot allocate vector of length 799980000 We can read the file without any problem with read.table, but trying to convert it to the geodata class gets the same error. Any help/tips would be appreciated. Kris, Dave --

Hola Javier, si, y es un factor de Y. Pero ya entendí lo que tenía mal, me falta convertir Nt a valores de geodata antes de hacer el plot, es decir: date<-as.geodata(Nt) plot.geodata(date) Con relación a la actualización de mi versión de R, no me he animado porque he tenido problemas con paquetes que no tienen las actualizaciones y no me corren bien. Gracias y saludos Beatriz 2013/8/6

I'm using R for Mac OSX version 1.14 (2129) and the geoR package version 1.6-5 (the current version in the R repository). I'm running R in OS 10.4.6 on a Mac G4 iBook (933MHz, 640 MB DDR SDRAM). I searched the R archive and did not find a posting on this issue. I want to use the variog and variog4 functions of geoR to characterize the pattern of spatial autocorrelation of tree

Hola a todos, estoy trabajando con dos archivos csv. En uno de ellos extraigo dos variables (que son coordenadas), del segundo archivo extraigo otra variable que transformo. Con estas tres variables construyo un data.frame (las tres variables tienen la misma longitud). Hasta ahí todo bien, solo que cuando quiero usar este data.frame para hacer un plot.geodata me marca el siguiente error Error en

ComRades: I am geeting the error message Error in ldots[[which(MET)]] : attempt to select less than one element when I try to fit the geostatistical model with the likfit() function of geoR. I have tried with old data for which likfit() successfully maximised the likelihood in previous versions of geoR, and yet the current version fails. I have tried in Windows Vista and Windows XP (I haven't

Hi, I'm writing some R code that requires a massive amount of looping and would ideally like to write it so that it avoid the use of "for" loop ... however I'm having some trouble. Very briefly, the basic idea is to implement a binary partitioning algorithm to determine the optimal cutpoint based on deviance measures obtained from likelihood estiamtes. This is in the

Hi, everybody! i created a imagem by kriging using geoR package. I imported points from GRASS("zn", after converted to geodata "zn_geo"), the border "zn_border" and a raster mask. Then i interpolated the points by kriging and created a raster image. Now, i need export this image back to GRASS to use it in the module r.mapcalc. I can't do it. I tried use

Dear GeoR-er, If I use the variog function in the latest release of geoR, the first lag is always ignored. For instance, if you read in geodata, calculates the variogram using the variog function and give in a uvec like uvec=seq(0,max,by=2.44), it only starts giving results from distance=4.88 and ignores 2.44! This wasn't the case in former versions of geoR. Is this done on purpose? Why?

I used points graphic (geoR) in this way: points(zn,pt.divide=c("quintiles"), main="Location map of Zn",lambda=1,col="gray",yl="coord y",xl="coord x",x.leg=0. 3, y.leg=5.5,dig.leg=2,cex.min=0.5,cex.max=1,) Can I reduce legend characters (a sort of cex)?

Dear colleagues I found duplicate coordinates trying to convert a dataframe object (see the attached file) into a geodata one ('geoR' package is needed). My dataframe object consists of 2 columns with UTM geographical coordinates and a 3rd column with fish densities . To overcome the handicap I must add a tiny value generated by rnorm() to each geographical coordinate. And it must be in

Hi R-people: Working in geoRglm, it shows me, according to AIC criterion, that the non-spacial model describes the process in a better way. It's the first time that I'm facing up to. These are my results: OP2003Seppos.AICnonsp-OP2003Seppos.AICsp #[1] -4 (OP2003Seppos.lf0.p<-exp(OP2003Seppos.lf0$beta)/(1+exp(OP2003Seppos.lf0$beta))) #P non spatial #[1] 0.9717596