如何对块克里金进行交叉验证？

Question

我在automap包中编写了一个代码来交叉验证不同的克里金技术。我已经交叉验证了所有这些，但我无法编写Block kriging的代码。它显示了这个错误：

未使用的参数（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)

现在，请让我知道如何解决它。修复此问题后的最后一步将是所有cv s的compare.cv，我是对的吗？

Answer 1

这意味着autoKrige.cv命令不会考虑数据块。它逐点执行交叉验证，而不是按块执行。

交叉验证考虑了POINTS的插值（或预测）估计的准确性，而块克里金法是一种平滑方法，将整个区域划分为多个BLOCKS并计算每个区域的估计的局部平均值。换句话说，对于区域＆＃39;块＆＃39;你不具备“价值”＆＃39;比较你克里金的估计