[R-sig-Geo] Performing a kriging

Komine momadou at yahoo.fr
Thu Jun 13 21:17:41 CEST 2013

Thank you Saman and Paulo, 

Certainly , I use Awkwardly the example advised by Paulo.  I don't
understand the code :var(Humidity) and 0.05*max(dist(cbind(X,Y))).
However, after the change of my data to geodata, I have not an error message
but the result is a black plot (useless).  
 The followed code was used: 

> library(geoR) 
> tab<-as.geodata(tab,coords.col=1:2,data.col=3,data.names=NULL) 
> ml <- likfit(tab, 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(tab, loc = pred.grid, krige = krige.control(obj.m = ml)) 
> image(kc, loc = pred.grid, col=gray(seq(1,0.1,l=30))) 

Thank you in advance for your advice. 

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

More information about the R-sig-Geo mailing list