Indicator kriging is nothing but ordinary or simple kriging on binary (0/1) data; there's quite a list of packages that can do kriging: gstat, geoR, sgeostat, fields, spatial, ... Be aware that indicator kriging outcomes only estimate probabilities and may lie outside [0,1] To my knowledge there's no package that ports GSLIB to R. -- Edzer