如何修复autoKrige.cv错误?

时间:2014-08-16 10:07:57

标签: r automap

任何人都可以在下面的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点的网格数据。

1 个答案:

答案 0 :(得分: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)

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