search for: krigged

hello I try to build DEM using Krig function of fields package. And I get this error message. here is the procedure I followed: > dt<-read.table("/usr/local/bartin/stat/topostat1",header=F,sep="|") > names(dt) <-c("x","y","z") > coord<-cbind(dt$x,dt$y) > elevation<-cbind(dt$z) >

Hi, I have a 4 column file (long/lat/elev/variable) and I tried to fit the values of my variable to the XYZ space using Tps and I keep getting the following message: Warning messages: 1: GCV search gives a minumum at the endpoints of the grid search in: Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, 2: GCV search gives a minumum at the endpoints of the grid search in:

I'm new at this, so please bear with me. I am trying to get Krig (actually tps) results into a data frame at a set of points. Frankly, all I really need is to query the Krig on a specific XY pair and return a result. I can carry the rest. I am using tps from the fields package. Here's the code I use to create the Krig object: =========== Snip! =============================== dCurr

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Mike J Smith > Sent: 09 June 2005 09:58 > To: r-help at stat.math.ethz.ch > Subject: [R] krig.image help > > > Hi > > I have recently been experimenting with the use of kriging, primarily > through Goldensoftware's Surfer.

Dear R-helpers, I'm having trouble compiling R packages from source on Win2K. I installed R 1.7.1beta [1] on my win2k machine [2], downloaded the fields package as source [3] and tried but failed to install the package [4]. I get the same problem with other packages, so it doesn't appear to be limited to fields. Installation of precompiled packages seems to work fine. I've been

In the code to follow, I'm trying to label points with their corresponding values but have been unsuccessful after many variations to the following code. The code below creates the plot I want, I simply cannot get the black points ("+") to display the actual value. I'm guessing the problem is somewhere in the second to last line of code (starts with "pts<-..."). I

I am returning to using R after a little time away. Saw some great new functions in the fields package and decided to upgrade. Everything seems to go fine until I tried to actually load the package, when I get an error: >library(fields) Error in parse(file, n, text, prompt) : syntax error on line 227 And (of course) I am unable to access any of the functions in the fields package. Below is

Dear list, I get the following error message when installing the package GRASS. missing link(s): interp.new missing link(s): vgram Krig predict.Krig predict.se.Krig See also below for the full installation message. What does this mean exactly? Regards, Ulrich ------------------------------------------------------------------ > R CMD INSTALL

I have a question to everybody. After log10 transfprmation, I have done ordinary kriging in gstat in R? I need to back trnasform the prediction values to orgiginal scale. How do I do this in gstat in R? Thanks Zia -- Zia Uddin Ahmed 915 Brad Field Hall Department of Crop and Soil Cornell University Ithaca NY 14850 USA

FreeBSD 5.2.1 Upgraded to dovecot 0.99.12.1 After authentication failed, I also upgraded OpenLDAP to 2.2.19 Dovecot and OpenLDAP are now the most recent versions in the ports. I can access the LDAP info through the OpenLDAP client without a problem, so slapd seems to be working fine. Both dovecot and OpenLDAP worked fine before the upgrade. Any help will be appreciated, since several people

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints

A majority of my data makes a kriged map perfectly using an exponential model for the semivariogram to fit my data and then going through the commands variofit() to define the model and then krige.conv() to use the model to predict values in a grid. But?one set of my data appears to be linearly correlated for the first 5000 meters and not correlated beyond that. I have been having problems

Hi ! Sorry for my ignorance. I have two questions: 1) which is in R the function to make a covariogram of spatial data ? 2) does anyone know if there exist the possibility in R of performing kriging with an arbitrary covariance function ? If I have well understood the Krig procedure offers only three possibilities (expcov, gauscov and sphercov) that are the common models for stationary and

Thanks for all the helpful responses. I include the data file and the syntax file for reference. Again, if I use the fields function, as is, I get the message: Error in svd(tempM) : error 159 in dsvdc using traceback, I get: > traceback() 4: stop(paste("error ", z$info, " in dsvdc")) 3: svd(tempM) 2: Krig(x, Y, cov.function = rad.cov, m = m, decomp = decomp,

Hi, I am trying to do a bayesian prediction for soil pollution data above a certain threshold, using geoR. Everything is working fine until i am doing the krig.bayes. I tried to do the prediction on a grid 67 by 113 cells and my computer is freezing to death. At larger numbers of cells it tells me after a while that it reaches the max. memory of 511 Mb. My computer has only 512 Mb of RAM.

Hello to everybody ! I don't have a question directly connected with R, but it takes starting point from the use of the R function Krig for kriging interpolation and from the idea of making a different implemetation of the function for my purposes. Doas anyone know where can I find free C or C++ codes of kriging ? As always, thanks a lot in advance, Paola.

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

Hello R users, I have numerical data sampled on two grids, each one shifted by 0.5 from the other. For example: grid1 = expand.grid(x=0:3, y=0.5:2.5) grid2 = expand.grid(x=0.5:2.5, y=0:3) gridFinal = expand.grid(x=0.5:2.5, y=0.5:2.5) plot(gridFinal, xlim=c(0,3), ylim=c(0,3), col="black", pch=19) points(grid1, xlim=c(0,3), ylim=c(0,3), col="red", pch=19) points(grid2,

Hi R community I need to interpolate precipitation data for a natural park. I have precipitation data from some climate stationts. (I know the table is not complete but I only need to show you X,Y, Altitude and PrepJul) X Y Altitude PrepJan PrepFeb PrepMar PrepAp PrepMay PrepJun PrepJul 597706 4093438 41 0 9 77,8 63,1 17,5 0 2,6 597535

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