[R] still have problem with krige and border option
l.bargeot at educagri.fr
Tue Apr 16 18:08:04 CEST 2002
I would like to estimate rainfall with the krige function. First, I produce a
polygon of my region of interest (where poly_test.txt is a x,y suite of
points defining a polygon, obtained with the grass v.out.ascii command) :
So now I have a list of x,y coordinates to limit my work area.
If I call the function krige with extrap=FALSE like this
[ where grid2.point is a point object to store prediction, logtab.point a
point object containing points data, 'lz' is the variable name in
logtab.point, logtab.Exp is a variogram object ]
all works quiet well, but I am limited to a convex hull inside the limite
designed by points of measure.
If I set extrap=TRUE, all works, but values are estimated all over my area.
Now if use the polygon I defined, and specifie it in the function call
all points in the area are estimated.
I've tried to put extrap=FALSE with border, but then there is no result, and
no points are found.
Here is my polygon definition:
 2272048 2275274 2296044 2293826 2280517 2272048
 781270.4 801031.7 795385.6 773003.0 765945.4 781270.4
So my question is, do you see any stupid mistake I can make ?
last point (optional!) about krige.G
I also tried to use krige.G from the grass package, but for the moment I have
a strange error message : (where G is obtained with gmeta() command)
Error in surf.gls(0, expcov, x = point.obj$x, y = point.obj$y, z = at.val, :
Rank failure in Choleski decomposition
I've just reported my point and variogram definition in the options. Perhap's
is there a small difference between parameters when they are passed to krige
function but I don't see wich.
Thanks in advance for all kind of help
4 rue champs-prevois
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