任何人都可以在下面的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点的网格数据。
答案 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)
我试图在没有可重现的数据集示例的情况下回答。下次请提供一个。