[R-sig-Geo] indicator simple kriging with varying local mean on gstat

[Edzer Pebesma]

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

facrimas at 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 at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>