[R-sig-Geo] universal kriging in gstat

Edzer J. Pebesma e.pebesma at geo.uu.nl
Mon Jul 3 17:40:27 CEST 2006

Alí Santacruz wrote:
> Dear list members,
> I am trying to perform universal kriging in gstat, and according to the help 
> of 'krige' function, the names of the coordinates of 'newdata' argument must 
> be the same as defined in 'locations'.
I suspect it is time to update this documentation; If you're using sp
objects, the names of coordinates does not play a role anymore, as
we can retrieve them using the coordinates() method.

Does this help?
> I am creating an SpatialPixels object for the prediction locations, using 
> the next code:
>> pts<-spsample(myShapePoly, n=3000, "regular")
>> gridded(pts) <- TRUE class(pts)
> [1] "SpatialPixels"
> attr(,"package")
> [1] "sp"
> And x1 and x2 are the names of the coordinates of the prediction locations 
> grid.
> So, to perform universal kriging, I have to change the names of the 
> coordinates of my observed data to x1 and x2, but I would like to know if I 
> can change the name of the coordinates of the prediction grid ('pts').
> I appreciate any comment. Sincerely,
> Alí M. Santacruz
> M.Sc. Geomatics, Student
> National University of Colombia
> Bogotá D.C.
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

More information about the R-sig-Geo mailing list