Hi to all. I made some research about this issue and found some helpful information (maybe only for me, because Iâm not an expert in this field). Automap mainly depends on gstat. In a very helpful book (Bivand, Pebesma, and Gómez-Rubio 2008. Applied Spatial data Analysis with R), borrowed from a friend of mine, I found this example for performing Indicator Kriging with gstat library:
ind.f <- I(zinc < 500) ~ 1  ind.fit <- fit.variogram(variogram(ind.f, meuse), vgm(1, "Sph", 800, 1))  ind.kr <- krige(ind.f, meuse, meuse.grid, ind.fit) [using ordinary kriging] The continuous variable is transformed into a binary one (TRUE/FALSE) and used in the interpolation process. I used a similar procedure with the function âautoKrigeâ in automap and got the prediction map with no error message. I suppose everything went well Cheers. Dr. Francesco Carotenuto (Ph.D.) Department of Earth Sciences Federico II University Largo S. Marcellino 10, 80138 NAPLES (ITALY) http://francesco-carotenuto.blogspot.it/ http://3dpaleontology.blogspot.it/ ________________________________ A: "r-sig-geo@r-project.org" <r-sig-geo@r-project.org> Inviato: Venerdì 13 Settembre 2013 11:48 Oggetto: Indicator Kriging in automap Dear list, What parameter set allows one to perform automatically Indicator Kriging in automap?  Dr. Francesco Carotenuto (Ph.D.) Department of Earth Sciences Federico II University Largo S. Marcellino 10, 80138 NAPLES (ITALY) http://francesco-carotenuto.blogspot.it/ http://3dpaleontology.blogspot.it/ [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
