[R-sig-Geo] Performing a kriging

[Komine]

Dear,
I use the data (.csv) shown below for performing a kriging.
       X             Y       Humidity
1 719134.1 1433671       90
2 736867.3 1431585       50
3 751471.1 1422197       32
4 749384.8 1414895       75
5 732694.8 1413851       12
6 719134.1 1412808       65

I use this code retreived from internet and adapted with my data.

> tab<-read.table("C:\\Users\\Kriging.csv",sep=";",dec=",",header=TRUE)
> tab
> library(geoR)
> ml <- likfit(Kriging, ini = c(1,0.5), fix.nugget =T)
> pred.grid <-  expand.grid(seq(0,1, l=51), seq(0,1, l=51))
> kc <- krige.conv(Kriging, loc = pred.grid, krige = krige.control(obj.m =
> ml))
> image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30))

When I run the code, I have an error message from the fourth command line. 

Thank you in advance for your advice, I am not familiar with in spatial
statistics.  

Komine  



--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Performing-a-kriging-tp7583773.html
Sent from the R-sig-geo mailing list archive at Nabble.com.