0

我想跑来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个观察点

4

0 回答 0