I am trying to perform an ordinary kriging with the above data using the 'krige.conv' function in the geoR package. I am however, encountering problems with successfully completing this. I have provided the data and code I have been running and would welcome some suggestions.
# Defining a prediction grid
# predicting by ordinary kriging
vario <- variog(sb)
vf <- variofit(vario)
kc <- krige.conv(sb, loc=grid, krige=krige.control(obj.m=vf))
I never worked with geoR package for interpolation pourposes. Have you already checked the gstat package ? Further, do you have any idea of which model would be the best for fitting the data you're working on (exponential, spherical, gaussian, etc ?).
Give me more information so I can help you better.
PhD italian National Council of Research - Institute of Methodologies for Environmental Research CNR - IMAA
85050 Tito Scalo - Potenza - Italy