我想跑来krige.cv()获取一组 LOOCV 数据。
它给了我一个错误,我看不到问题。这可能是什么问题?
这是代码:
## LOAD DATA
comp_150_sg <- read_excel("Geos_descanso/150/comp_150_sg.xlsx")
View(comp_150_sg)
str(comp_150_sg)
## SPATIAL OBJECT
comp_150_sg_sp2 <- comp_150_sg
coordinates(comp_150_sg_sp2) <- ~x + y
SG=comp_150_sg_sp2$SG_150
f=SG~1
variog = autofitVariogram(f,comp_150_sg_sp2)
krige.cv(f,
locations=comp_150_sg_sp2,
model=variog$var_model)
这是 Excel 文件的样子:
它有201个观察点
