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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Le Lundi 15 Avril 2002 21:04, o.christensen at lancaster.ac.uk a ?crit :> It is difficult to say exactly what your problem is, because you haven't > provided the call you used.Yes, I'm confused, it was the end of the day.> 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'' .You're right ! My call was grid2.krige<-krige(grid2.point,logtab.point,"lz",logtab.Exp,extrap=T,border=li) without maxdist and point.obj Thanks for all 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._