[R] krige and polygon limit problem

Ole Christensen o.christensen at lancaster.ac.uk
Mon Apr 15 21:04:28 CEST 2002

Dear Lionel

It is difficult to say exactly what your problem is, because you haven't
provided the call you used.
But see my guess below.

> I'm new on R and this mailing list. We work on spatial rainfall estimation
> with R and Grass.
> We have a problem with the krige function from the sgeostat package. We would
> like to limit the estimated area with a polygon limit.
> I use a 50 points polygon to describe my work area. The krige function work
> quiet well without limit. But if I use this option I have the following error
> messages :
> -------------------------------
> Using points within c(2210907, 2212817, 2214516, 2216568, 2216985, 2215795,
> 2217396,  units of prediction points.

This should be a number (``maxdist'' in the call) not a vector. 
My guess is that you have given the polygon-coordinates as input to
``maxdist'' instead of to ``border'' .

Cheers Ole


Ole F. Christensen
Department of Mathematics and Statistics
Fylde College, Lancaster University 
Lancaster, LA1 4YF, England
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch

More information about the R-help mailing list