Hello library(fields) x<-cbind(dt1$V1,dt1$V2) # V1:east , V2=north coordinates of points Y<-cbind(dt1$V4,dt1$V5) # V4,V5 cos and sines components fit2 <- Krig(x,Y,cov.function=exp.cov, scale.type="range") : this gives "Error in if (abs(f2 - f1) < tol) { : missing value where TRUE/FALSE needed" what should I have done ? is interpolation of x,y components as a vector ok ? or what would you suggest for interpolation of the vector Y ? I will appreciate , if you tell me what I have to do Ahmet Temiz General Directory of Disaster Affairs Ankara TURKEY ______________________________________ Inflex - installed on mailserver for domain @deprem.gov.tr Queries to: postmaster at deprem.gov.tr ______________________________________ The views and opinions expressed in this e-mail message are ...{{dropped}}