0

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

未使用的参数 (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 的 compare.cv ，对吗？

4

1 回答 1

2

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

交叉验证考虑了 POINTS 的插值（或预测）估计的准确性，而块克里金法是将整个区域划分为几个块并计算每个区域的估计的局部平均值的平滑方法。换句话说，对于“块”区域，您没有“值”来比较克里金法所做的估计

于 2014-08-18T12:58:35.060 回答