[R-sig-Geo] Note on variogram fitting with polynomial degree
caspar.hallmann at gmail.com
Tue Sep 17 11:36:13 CEST 2013
I read in ASDAR: "Residual variograms are calculated by default when a
more complex model for the trend is used..."
However, I noticed that when using fit.variogram in gstat with an
additional option: degree=2, but keeping the formula as an intercept
only, does not account for the higher degree of complexity in the
model. I suspect this is not the intended behavior. To honour a
dependence structure (polyn degree) in my observations I need to state
the dependence explicitly in the formula.
coordinates(meuse) <- c("x", "y")
R version 3.0.1 (2013-05-16)
More information about the R-sig-Geo