[R] krige and polygon limit problem

[Lionel]

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
-- 
CNERTA-ENESAD
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._