similar to: calculating variograms (gstat) with large data sets

Hi, I found the same question in the mailing list already a few months ago - but there was no answer to it - so I'll try it again Could somebody help me to solve this following problem? I just begin to learn how to connect my Oracle database with R. > library(DBI) > library(ROracle) Warning message: DLL attempted to change FPU control word from 8001f to 9001f >

Hello All! I calculate a variogram using the function variog (package geoR) afterwards I use variofit to fit a spherical model (see code below). Now I just changed the units of the variable (in this case MPa to kPa just a factor of 1000). If I do so, I get a different fit and therefore different ranges etc. Why? The semi-variance is of course 6 orders of magnitude higher but the values

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

Hello, We are looking to adjust the font size of the axis annotation on the graph that results from use of the metaplot() function. Metaplot seems to respond to cex and cex.lab to change those graphical parameters, but it doesn't respond to cex.axis. Is there a way to work around this by creating a customized x-axis, and if so, how? Thanks for all your help. Syntax is below. Best, Dawn

Hello group, For my PhD I'm working on a spatial sampling grid. I do have two data sets which I'd like to compare using cross-correlograms or cross-variograms. Is this an option in one of the R-packages? I've been searching the R-help archive and the available package-documentations, but I can't find how to do this. Thanks in advance, Ren?.

Hello Trying to develop variograms and kriged surfaces from a point file. Here is what I've done so far. library(gstat) # also loads library(sp) library(lattice) soilpts$x <- soilpts$UTM_X soilpts$y <- soilpts$UTM_Y soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN, TC, Total_Mg)) dim(soil.dat) [1] 1292 7 coordinates(soil.dat) <- ~ x+y

Jacques, provided that X and Y are colocated (i.e., have exactly the same observation locations), you get the cross variogram right; the definition of this cross variogram is however: gamma(h)= E[(X(s)-X(s+h))*(Y(s)-Y(s+h))] also, where you select: cv <- v$gamma[1:14] you may be better off using the more general v$gamma[v$id == "X.Y"] Best regards, -- Edzer

Hello R list, I need to compute empirical variograms using data from a large geographic area (~10^6 km2). Although I could not find a specific reference, I assume that both geoR and gstat calculate distances among data points assuming points are on a flat surface (using the Pythagorean Theorem). Because the location of my data is large and located near the pole, assuming that latitude and

Hi, I would know how to plot two variograms on a same graph. I can plot one by one but I would draw both on the same one. Is it possible? Do i need any special package? Thanks! Cordialement Damien Landais

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 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, I am new in R and trying to implement an algorithm which makes ordinary kriging by using gstat library and the predict method. I use the predict method as following: First, I create an object g: g <- gstat(id="tec", formula=TEC ~ 1, data=data) ## Create gstat object called g And then I use this object in the predict. p <- predict.gstat(g, model=mod, newdata=predGrid,

Hi, I am a student at the University of Twente ( ITC). I am using the R packages for my data analysis, but for the last few weeks now , I have been getting the error message when trying to work on variograms or krigging. Error : .onLoad failed in loadNamespace() for 'gstat', details: call: fun(...) error: .Random.seed is not an integer vector but of type 'list' Error:

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

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

Hello everybody, I am a new R & Grass user. I am interested in doing a kriging of the Temperature in a basin, based on pixel elevation (DTM). I use Grass 6.2.1 and R 2.5.1 and have installed the spgrass6 package. The DTM is called "DTM_trentino_100m" and is like this: north:5158000.0 south:5058500.0 east:1729700.0 west:1611600.0 rows:995 cols:1181 -9999.0 -9999.0 -9999.0

Hi there, I have a problem using gstat when attempting to fit a regression-kriging model to my data. I have successfully carried out a stepwise regression on my data and I have the model residuals. However when I try to apply a variogram to the residuals I get the following error message; *>>>Error in grid.newpage() : * * VECTOR_ELT() can only be applied to a 'list', not a

dear all i have this dataset: x,y, datavalue > dati[,c(1,2,5)] [,1] [, 2] [,3] [1,] 2.386 3.077 1.740 [2,] 2.544 1.972 1.335 [3,] 2.807 3.347 1.610 [4,] 4.308 1.933 2.150 [5,] 4.383 1.081 1.565 [6,] 3.244 4.519 1.145 [7,] 3.925 3.785 0.894 [8,] 2.116 3.498 0.525 [9,] 1.842 0.989 0.240 [10,] 1.709 1.843 0.625 [11,] 3.800 4.578 3.873 [12,] 2.699 1.199 1.425

Hello, # I am able to install.packages("gstat") #but when I try to upload I get an error message library(gstat) #Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : # there is no package called ‘intervals’ #In addition: Warning message: #package ‘gstat’ was built under R version 3.0.2 #Error: package or namespace load failed for ‘gstat’ what do

Schlatter Christian wrote: >Dear list members > > > >I'm very new to R but a little informed about geostatistics. > >As I was looking for possibilities of geostatistical analysis in R I encountered at least two very interesting packages: > > > >Rgeo and gstat > > > >And of course I'm wondering now about the one which fits better my needs