Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-12 Thread Paul Hiemstra
till does not work. cheers, Paul > > # Error in predict(k, testpoint) : argument "k" is missing, with no default > > > regards, > Jeff > > -------------------- > >> Date: Fri, 11 Jun 2010 15:04:32 +0200 >> From: p.hiems...@geo.uu.nl >> To: vivac...@hotmail.sg >&g

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-12 Thread v v
dict(k, testpoint) : argument "k" is missing, with no default regards, Jeff > Date: Fri, 11 Jun 2010 15:04:32 +0200 > From: p.hiems...@geo.uu.nl > To: vivac...@hotmail.sg > CC: r-sig-geo@stat.math.ethz.ch > Subject: Re: [R-sig-Geo] H

Re: [R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-11 Thread Paul Hiemstra
Dear Jeff, I think you cannot escape discretizing your kriging surface as there is no mathematical expression that captures the kriging surface. To use optim you could make an objective function that was two inputs, x and y location. The functions calls the fields package to estimate the krigi

[R-sig-Geo] How to locate local minima points on a Krige-fitted surface?

2010-06-10 Thread v v
Hi, I am using the 'fields' package to conduct Kriging interpolation on my dataset. I have managed to generate the 3-dimensional surface plot of the Krige-fitted surface. There are several local minima points on the surface. My main interest is to pinpoint the exact x-y coordinates of the low