ordinary kriging using 'krige.conv' function

[djrwicks]

Moran_XY.csv

List readers,

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.

data<-read.csv("Moran_XY.csv")
sb<-as.geodata(data)
# Defining a prediction grid
grid<-expand.grid(seq(300000,360000,l=60),seq(3613000,3650000,l=60))
# predicting by ordinary kriging
vario <- variog(sb)
vf <- variofit(vario)
kc <- krige.conv(sb, loc=grid, krige=krige.control(obj.m=vf))
#generating image
image(kc)

Thanks in advance,

Dan.

[Giuseppe Calamita]

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.
Cheers

Giuseppe Calamita