r - 使用 R 中的交叉验证 krigecv 在 gstat 中出现长度错误

Question

今天我遇到了一个我找不到解决方案的问题。我想在 R 中使用 gstat 使用 krigcv 交叉验证，但我收到一条错误消息：

model.frame.default(terms(formula), as(data, "data.frame"), na.action = na.fail) 中的 0% 错误：可变长度不同（为 'x' 找到）

这是我的代码行：QCini.uk.cv <- krige.cv(QCini$var1 ~ x + y + I(y^2) + I(y^3) + I(x * y^2) + x: y, QCini, Vresfit.fi.vgm)

Answer 1

避免在公式参数中使用 data.frame，然后尝试

QCini.uk.cv <- krige.cv(var1 ~ x + y + I(y^2) + I(y^3) + I(x * y^2) + 
  x:y, QCini, Vresfit.fi.vgm)