Dear List,

I would like to make predictions by a gwr-kriging model and estimate it in a simultanious way. Does anyone know how to estimate the variogram and then run regression-kriging with gwr trend component using package 'gstat'? I could not find any source about how to add the gwr model as a trend to kriging in R. 


I attach the two model components that I've already run separatly and I would like to combine and estimate in a simultanious way:


# GWR model:
bw<-gwr.sel(r_reg_mn~r_vallalk+mobilit+szaporul+r_fiatal+r_ffi, mn.adat, adapt=TRUE, gweight=gwr.Gauss, method="cv", longlat=FALSE)
mn.cv.gwr<-gwr(r_reg_mn~r_vallalk+mobilit+szaporul+r_fiatal+r_ffi, mn.adat, adapt=bw, gweight=gwr.Gauss)

# ordinary kriging:
v<-variogram(r_reg_mn~1, mn.adat)
v.fit.ok<-fit.variogram(v, vgm(0.0475, model="Lin", nugget=0.003))

ok<-krige(r_reg_mn~1, locations=mn.adat.m, newdata=loc, model=v.fit.ok)

Thanks,
Laszlo 

	[[alternative HTML version deleted]]

