0

谁能在下面的 r 代码(automap 包)中找到错误,或者给我一些关于自动交叉验证的建议?

library(automap)
mydata<-read.table(".../mydata.txt", header=T, sep=",")
x<-mydata[,1]
y<-mydata[,2]
data1<-mydata[,3]
mydata.grid<-read.table(".../mydata.grid.txt", header=T, sep=",")
coordinates(mydata)=~x+y
gridded(mydata.grid)=~x+y
mykr.cv<-autoKrige.cv(log(data1)~1, mydata, model=c("Ste"), nfold=10)

但它显示了这个错误:

0% 错误:尺寸不匹配:位置 120 和数据 64

我使用的空间数据有 3 列 (x, y, mydata) 64 个点,2 列 (x, y) 中的网格数据包含 97868 个点。

4

1 回答 1

1

加载 mydata.txt 后,尝试重命名数据框的列而不实例化新向量:

colnames(mydata) = c("x","y","data1")

然后:

coordinates(mydata)=~x+y
mykr.cv<-autoKrige.cv(log(data1)~1, mydata, model=c("Ste"), nfold=10)

我试图在没有可重复的数据集示例的情况下回答。下次请提供一份。

于 2014-08-16T11:19:45.967 回答