[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
--
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list