Hello,
I would like to do indicator kriging in R/gstat. I have dichotomous 0/1
data and have performed ordinary kriging and universal kriging, but get
predctions that are far from 0 and 1.
Am I doing something wrong? Here is the code I have been using:
#### Ordinary Kriging with krige function #####
# ordinary kriging:
x <- krige(z~1, ~x+y, model = v.fit, data = estand, newd = grd)
### universal block kriging:
uk <- krige(z~x+y+DistHF+RiverDist+RiverDist2, ~x+y, model = v.fit, data
= estand, newdata =
grd)
Any help would be appreciated.
Pete
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
