The simplest solution might be to just (simple) krige the residuals, and add them to your known varying mean.
Please be aware that this "model" suggests residuals have mean-dependent variance (p * (1-p)), and thus using a single variogram model for the full set of residuals is not in accordance with theory. I guess that this problem is larger than the one you mention -- covariances are simply derived from variograms that reach a sill by C(h) = gamma(inf) - gamma(h). -- Edzer facri...@libero.it wrote: > Dear guys, > > I want to use the indicator simple kriging with varying local mean on gstat. > > Do you have a script example about it? > > I do not understand which is the local mean that I have to put in the command > of simple indicator kriging..... > > I do not also well understand whether if it is posssible to use the > semivariogram model in SK, because theory tells that it need the covariance > function....... > > Any suggest? > > Thank you. > > Cristiano > > _______________________________________________ > R-sig-Geo mailing list > R-sig-Geo@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
