[R-sig-Geo] Getting regression coefficients when universal kriging

Steinbuch, Luc luc.steinbuch at wur.nl
Tue May 27 10:50:50 CEST 2014


Suppose I have this  (working) code:

library(gstat) 
library(rgdal) 
library(raster) 

data(meuse) 
coordinates(meuse) = ~x+y 
r <- raster(system.file("external/test.grd", package="raster")) 
projection(r)<-projection(meuse) 
names(r) = 'elev' 

vu<-variogram(log(zinc)~elev,meuse) 
mu<-fit.variogram(vu,vgm(1,"Sph",300,1)) 
gUK<-gstat(NULL,"log.zinc",log(zinc)~elev,meuse,model=mu) 
UK<-interpolate(r, gUK, xyOnly=FALSE)

How do I find the finally used regression parameters, calculated by the universal kriging function of 'interpolate'?



More information about the R-sig-Geo mailing list