[R-sig-Geo] Indicator Kriging R/gstat

[Pete Larson]

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