similar to: gstat 0.9-12: cokriging cross validation and class name incompatibilities

Dear All, I run following code to estimate the blocks using cokriging ( my data set has more than 50,000 data points). All the things run finely but Once I run the predict.gstat function it gave the error message - "memory.c", line 58: can't allocate memory in function m_get(). I run this code on LINUX sever but result is same. Would any one please be able to give a solution for

Hi, please help!! I've been trying to fit variogram models using the REML method in the gstat package. Every time the Windows GUI crashes. For example library(gstat) data(meuse) x <- variogram(zinc ~ 1, ~x + y, meuse) v <- vgm(140000, "Sph", 800, nug = 10000) plot(x, model = fit.variogram(x, model = v, fit.method=5)) Other fit methods are non problematic (eg. fit.method=7

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

The following came up when Roger Bivand and I discussed R's name spaces and overlap in packages, in the bus to the field trip during StatGIS, last Tuesday: If two packages create the same class, say "variogram", and both are loaded, then using a method for an object of class "variogram" cannot discriminate between them and will call the method in the package loaded last.

Much like time, spatial locations are not likely to change for observed data. Still, R has no infrastructure or knowledge about spatial locations. Many packages that deal with spatial data in R use their own classes to deal with spatial data. On the workshop on spatial statistics software held during DSC2003 (organized by Roger Bivand), we decided that a base class that defines classes for

Dear list, I am having troubles increasing the fontize when plotting a variogram{gstat} and its model (vgm) with plot and using jpeg(). Also the resolution in the jpeg call does not work. I am using R2.4 on Mandriva 10.2 linux. I can change fontsize with cex.axis in a normal plot, so I presume it has to do with plotting the variogram model. Any help on how to increase the font size and

Hi all-- I'm new to R, statistics and programming, so sorry if this is a really basic question! I have plotted a directional variogram, and I want to a. overlay the omni-directional line over each directional panel b. display the directional variograms in a single panel with a legend that associated each line to each degree measurement. The line I'm using is

I wonder what is the default value for the argument 'cutoff' when not specified in the variogram.formula function of gstat. Computing variogram envelops within gstat, I am comparing the results obtained with variog in geoR and variogram in gstat, and it took me a while before understanding that the cutoff default value is not the maximum distance. Can Edzer tell us about it? All the

Dear all can you help me? i have this problem: i have a dataset in a text file in a matrix of 3 columns: x, y, z where x and y are the coordinates and Z are the mesurements. How can i do a variogram with R?

Hi ALL, With the pakage gstat How can one create a grid of points as in the meuse.grid example, from measuring sites as in the Meuse example Thanks in advance! Alain PONSERO Conservateur Réserve Naturelle de la Baie de Saint-Brieuc site de l'Etoile 22120 HILLION tel/fax : 02.96.32.31.40 <mailto:aponsero@cabri22.com> aponsero@cabri22.com

Hello, I'm using the gstat package within R for an automated procedure that uses ordinary kriging. I can see that there is a logical ("singular") atrtibute of some adjusted model semivariograms: .- attr(*, "singular")= logi TRUE I cannot find documentation about the exact meaning and the implications of this attribute, and I dont know anything about the inner calculations

Hello. I'm pretty much new to R and I'm trying to produce some figures. It seems to me, that R has some asynchronous way of plotting figures. When I run this code: #constructs the semivariogram of SC1929 vgm1 <- variogram(SC1929~1,~U+V,puerto.map$att.data) # trying to make new plot dev.set(which=dev.next()) plot(vgm1) title(main="Semivariogram",font.main=4)

The majority of the functionality present in the gstat stand-alone program (http://www.gstat.org/) is now available as a package/library for the S language (R, S-Plus), again called gstat. The package provides multivariable geostatistical modelling, prediction and simulation, as well as several visualisation functions. Gstat was started 10 years ago and was released under the GPL in 1996; the

The majority of the functionality present in the gstat stand-alone program (http://www.gstat.org/) is now available as a package/library for the S language (R, S-Plus), again called gstat. The package provides multivariable geostatistical modelling, prediction and simulation, as well as several visualisation functions. Gstat was started 10 years ago and was released under the GPL in 1996; the

Dear Members, I am trying to fit a variogram model using fit.variogram function from the gstat package. The figure showing my experimental variogram can be seen here: https://i.stack.imgur.com/UZXw4.png My code line for this operation is: > c2.vgm.fit<-fit.variogram(c2.vgm.exp,vgm(nugget=0, psill=400,model="Exp",range =40000),fit.method = 7) The system throws following error

Dear all, I am aiming to calculate variograms using variogram() from gstat. The problem is that some of my data-sets are very large (> 400000 points). Running the command takes some hours of time and does not give any error-message. Nevertheless the result seem not to be appropriate - the first few bins are ok (up to a distance of about 300) but then it gains lags which are much larger than

Hola, # Hacemos el KED. Ver función "krige()": KED.rad <- krige( formula=pluvPcp~layer, # covariable -> radar locations=lluvia.rad.pluv.spdf, newdata=radarGrid, # podría ser cualquier objeto Spatial model=v.fit, # modelo de semivariograma. maxdist=Inf

Sale plano sí. Ya se que sin tener los datos y el código es un poco difícil, pero es que mis datos ocupan mucho, es imposible. Seguiré mirando por internet. Muchas gracias Rubén. Un saludo, > To: r-help-es en r-project.org > From: rubenfcasal en gmail.com > Date: Thu, 6 Aug 2015 14:21:47 +0200 > Subject: Re: [R-es] Duda interpolación (package ' gstat ') > > Hola

HI, I wonder does anyone have experience with doing sequential gaussian simulation with krige() function in gstat? I find it VERY slow compared to use krige() to achieve kriging function itself.. I wonder why, is that because it has to model the variogram, and do the kriging separately for each point to be simulated? so it would be N times slower to achieve the simulation than the kriging if