similar to: How to use your own data in gstat and sp?

Hello I wonder if anybody can help, I am using the package adehabitatHR to estimate the potential distribution of a species using the command "domain" In the example given in the AdehabitatHS manual a map containing elevation information is loaded (class= spatial pixels data frame) as well as the GPS points of the animals being tracked, these are then plotted on each other and

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

Dear all, I'm working with satellite images in R and plotting them via the code below. I was wondering whether coordinates (spatial["x"], spatial["y"]) are used as centre coordinates of the pixels in the GRID? In this script; spatial["x"] & spatial["y"] are the centre coordinates of the satellite image pixels. I'm asking this because some

I have been using the sp package to deal with gridded data. I would like to know how to extract the vertices (x,y) of polygons that outline areas of like-valued cells in a grid. Here is a simple 3x3 grid: 2 2 1 1 2 1 1 1 3 x <- c(1,1,1,2,2,2,3,3,3) # define a 3 x 3 array of points y <- c(1,2,3,1,2,3,1,2,3) h <- SpatialPoints(cbind(x,y)) # make these an sp object of the points

R Users: ...been working with the sp and gstat packages for the past couple of days in an effort to analyze a set of ~ 200 soil samples collected from various eastings, northings, and depths and containing a wide range of measured hydrocarbon concentrations. Thus far, I've managed to import the data, log-transform the concentrations, assign coordinates, generate and fit a variogram model and

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

I uploaded gstat 0.9-12 to CRAN, which has a few important changes: 1. Cokriging cross validation Cokriging cross validation is now possible with the function gstat.cv: you simply pass a multivariable gstat object, and cross validation is done for the first variable in the object. Optionally, secondary variable records at locations coinciding with the validation locations are removed. 2. Class

I uploaded gstat 0.9-12 to CRAN, which has a few important changes: 1. Cokriging cross validation Cokriging cross validation is now possible with the function gstat.cv: you simply pass a multivariable gstat object, and cross validation is done for the first variable in the object. Optionally, secondary variable records at locations coinciding with the validation locations are removed. 2. Class

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

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,

Now I have the forest plot data with x, y locations, and I measured the DBH for every indivicuals. The data looks like that: x=runif(100,0,100) y=runif(100,0,100) dbh=runif(100,1,100) rdata=data.frame(x,y,dbh) > rdata[1:5,] x y dbh 1 99.5354145 1.412844 34.10112 2 0.8259361 87.737036 39.12914 3 6.5678613 65.699032 22.55990 4 67.2987881 72.053877 45.83978 5 2.2491372

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

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;

hello i'm trying the script of gstat course of munich "Spatial and spacetime classes in R " in http://geostat-course.org/Topic_Bivand_2012 the code is follow ... but a in timeIsInterval(dts) the script failed.... the message is "timeIsInterval function not found" i'm updated R and package spacetime but has not worked any suggestions? thanks setwd("C:/Documents

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

Hi, I'm using R 2.5.1 and gstat 0.9-39 library and I'm working with the jura data set provided by gstat library. I tried to plot a graph of metal concentrations (let's say Cd) with the command spplot, but I realized that the default lags are equally distributed between Min and Max. I did: library(gstat) data(jura) jura.pred.xy <- jura.pred names(jura.pred.xy)[1:2] <-

Ciao a tutti ho appena iniziato ad utilizzare R per ora per attuare un'analisi geostatistica di dati. Volevo sapere come poter creare un oggetto gstat partendo da un file testo(che ho gia importato con read.table)e che contiene 3 colonne: x,y,value. Mi servirebbe far questo per costruire un variogramma. So che la domanda molto probabilmente per voi sara' banale....vi ringrazio