similar to: non-linear trends in kriging model

Hi I have two geostatistical models from geoR. An ordinary kriging model with AIC=-148.6 and a universal kriging model with AIC=-156.7, there are 345 data points. The improvement shown by the AIC by adding a trend component to the model seems quite small given the number of data points, is there a test to see if the improvement to the model fit is significant? Thanks David

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

All, I am new to using R and know some basics. I wish to use kriging in R to do the following: given data Y =f(X1,X2,X3,.....,Xn) --1000+ irregular measured data set. I would like to be able to get a single value y given sinle input set (x1,x2,x3,...xn) A google search on this takes me lierally to the same example on involving analysis with soil sampling and I cannot figure out how to

R Experts, I'm looking for some help with the geoR package. I'm trying to krig some data without using a global neighborhood. I would like to set my moving neighborhood to a distance, say 100 meters, where I know my data is spatially correlated. I have tried the ksline function, but that only allows my moving neighborhood to be set to a number of data points. But, since my data

Dear all, When a conventional kriging and then a contour plot is limited with a polygon (as possible with krige.conv and contour.kriging), the polygon border is displayed in black by default. > kc<-krige.conv(CZdata,loc=pred.grid,borders=czpoly,krige=krige.control(obj.m=ls)) > plot(CZcoord,xlab="x",ylab="y",type="n",asp=1) >

Hi, I'm using the function krige.conv from package geoR but I'm getting a warning that I do not understand: Warning messages: 1: NaNs produced in: sqrt(variance) 2: NaNs produced in: sqrt(variance) I'm performimg an ordinary kriging with a log transform (lambda=0). Does anyone knows where this comes from ? Thanks and regards EJ

Dear Aaron, I am really a tool user and not a tool maker (actually an ecologist doing some biostatistics)... so, I take the liberty of sending a copy of this e-mail to the r-help list where capable computer persons and true statisticians may provide more relevant information and also to Paulo Ribeiro and Peter Diggle, the authors of geoR.. I really feel that your huge matrix cannot be

Hi I'm using geoR package to perform linear spatial interpolation(OK). The function likfit() fails to compute REML. The error meassage is : Error in solve.default(v$varcov, xmat); How I can find out that likfit() is failed to process and retrieving the error message ? Thank you so much for your help.

Hi list, Well, this time I’ve a doubt with mapping generation. I was already able to read and plot shapefiles, plot point on this map. All this without any problems, but now I want to do something that I think, only Golden Software Surfer is capable of. I would like to plot a kriging result on the existing map (map script below). Well I looked for on the net, but I was not able to be sure of

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

We want to compute a pseudo R-squared for a model whose parameter estimation was based on maximum likelihood (function likfit, package geoR). I tried to compute the R2 proposed by Maddala (1983) which compare the maximized likelihood for the model without any predictor and the maximized likelihood for the model with all predictors. I got a really low value (0,01%). Did I miss something? Are

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,

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060628/2ff7714d/attachment.pl

Hi, I want to automatically fit variograms to a large number of different sample data sets, and call the funtion "likfit" (in package geoR) from within a for-loop. "likfit" does again call "optim". After ssuccessfully fitting variograms to some of the data sets, the procedure crashes and I get the error message: Error in optim(par = ini, fn = negloglik.GRF,

Hi, all, Have anyone used kringing included in R? How is it? Does it handle anisotropy data well? How does it compare with Kriging in Arc/Info? or other geostatistics software customized to do kriging or other geostatistics functions? I tried Easykriging, a geostatistics tool developed for Matlab. It has very nice GUI, but it does not provide library which i can call in my programs. so it is

Hi all, I have to draw a map of microorganisms repartition in a sample. I use the "image" function of geoR package to build the map but there are three problems: 1) the legend is Horizontal (with or without vert.leg=TRUE) 2) I want to plot the position of organisms and level curves of density on the map but I don't find the script to insert these data on the map. 3)the grey range

Hello All! I like to do some indicator kriging with R. So far I used geoR for simple and ordinary kriging. Does anybody know which package I should use? Thanks for your help! Sascha! -- Jetzt dabei sein: http://www.shortview.de/wasistshortview.php?mc=sv_ext_mf at gmx

Hello, I'm interested in doing Bayesian kriging using R. I see that the package geoR has a function that will allow one to do this (krige.bayes). However, my data are not in the form of points, but rather they are blocks that represent spatial averages (i.e., the number of fishing hooks per month in a given lat x long square). I am therefore interested in treating the data as

Dear R list, I've just installed R and then followed the instructions to install the package "geoR" from within R by using the following line of code: install.packages("geoR", contriburl = "http://www.est.ufpr.br/geoR/windows") This installed okay under the folder C:\documents and settings\ ... But when I follow the instructions in the illustrative session

Hi R-Users Is there any package in R that work Indicator Kriging? Or somebody is working in a link between GSLIB library in FORTRAN for R? Thank you very much for your help. --