Как сделать перекрестную проверку для блочного кригинга?

Question

Я написал код в пакете automap для перекрестной проверки различных методов кригинга. Я перекрестно проверил их все, но не могу написать код для блочного кригинга. Это показывает эту ошибку:

неиспользованный аргумент (блок = с (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 всех cv s, я прав?

Answer 1

Это означает, что команда autoKrige.cv не учитывает блоки ваших данных. Он выполняет перекрестную проверку по точкам, а не по блокам.

Перекрестная проверка учитывает точность оценок интерполяции (или прогнозирования) для ТОЧЕК, в то время как блочный кригинг - это метод сглаживания, который делит всю область на несколько БЛОКОВ и вычисляет локальное среднее ваших оценок для каждой из этих областей. Другими словами, для области "блок" у вас нет "значения" для сравнения вашей оценки, сделанной кригингом