[R] krige and polygon limit problem

Lionel l.bargeot at educagri.fr
Mon Apr 15 17:03:48 CEST 2002

```Dear all,

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.

Using points within c(785291, 789980, 794030, 798282, 805920, 813151, 817000,
820902,  units of prediction points.
Predicting.Error: (list) object cannot be coerced to vector type 14
In addition: Warning messages:
1: X matrix was collinear in: lsfit(xmat, y, wt = w, intercept = FALSE)
2: X matrix was collinear in: lsfit(xmat, y, wt = w, intercept = FALSE)
3: X matrix was collinear in: lsfit(xmat, y, wt = w, intercept = FALSE)
4: X matrix was collinear in: lsfit(xmat, y, wt = w, intercept = FALSE)
-------------------------------

Do you have an idea of the origine of the problem ?

Thanks,

Lionel
--
4 rue champs-prevois
batiment grand-champs
21000 Dijon
tel:03.80.77.28.49
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._

```