r - how to do cross-validation for block kriging? -


i have written code in automap package cross-validate different kriging techniques. have cross-validated of them, cannot write code block kriging. shows error:

unused argument (block=c(400,400))

library(automap) mydata<-read.table(".../mydata.txt",header=t,sep=",") colnames(mydata)=c("x","y","data1") library(gstat) coordinates(mydata)=~x+y mygrids<-read.table(".../grids.txt",header=t,sep=",") gridded(mygrids)=~x+y block_kriging_cv<-autokrige.cv(log(data1)~x+y, mydata, mygrids, block=c(400,400), model=c("ste"), nfold=10) 

now, please let me know how can fix it. final step after fixing compare.cv of cv s, right?

this means autokrige.cv command not take in account blocks of data. performs cross-validation point-by-point , not blocks.

cross validation takes in account accuracy of estimates of interpolation (or prediction) points while block kriging smoothing method divides whole area several blocks , calculate local average of estimations each of area. in other words, area 'block' don't have 'value' compare estimation made kriging


Comments

Popular posts from this blog

javascript - Jquery show_hide, what to add in order to make the page scroll to the bottom of the hidden field once button is clicked -

javascript - Highcharts multi-color line -

javascript - Enter key does not work in search box -