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

R 如何对块克里金进行交叉验证?,r,cross-validation,automap,R,Cross Validation,Automap,我已经在automap包中编写了一个代码来交叉验证不同的克里格技术。我已经对它们进行了交叉验证,但我无法编写块克里金法的代码。它显示了这个错误: 未使用的参数(block=c(400400)) 库(自动映射) mydata这意味着autoKrige.cv命令不考虑数据块。它逐点而不是按块执行交叉验证 交叉验证考虑点插值(或预测)估计的准确性,而块克里格法是一种平滑方法,它将整个区域划分为几个块,并计算每个区域估计的局部平均值。换句话说,对于面积“块”,你没有一个“值”来比较你通过克里格法做出的估

我已经在automap包中编写了一个代码来交叉验证不同的克里格技术。我已经对它们进行了交叉验证,但我无法编写块克里金法的代码。它显示了这个错误:

未使用的参数(block=c(400400))

库(自动映射)

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


交叉验证考虑点插值(或预测)估计的准确性,而块克里格法是一种平滑方法,它将整个区域划分为几个块,并计算每个区域估计的局部平均值。换句话说,对于面积“块”,你没有一个“值”来比较你通过克里格法做出的估计

谢谢@Fabio,但毕竟,我怎样才能确保这种方法给出的估计不比简单、普通或通用克里格法好。因为块克里格方法的SSE图的值比其他方法的少。因此,我不知道如何决定和判断哪一种方法是最合适的插值方法,我只会使用简单、普通和通用的克里格法。并通过这三种方法的交叉验证,得到信息统计。块克立格法的交叉验证并不简单,我也这么做了。非常感谢@Fabio。亲爱的@Fabio,你能不能也写一个代码:1.计算空间平均数,2。计算克里格结果的标准递减或方差?谢谢
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)