automap包autofitVariogram函数cutoff选项不生效问题咨询

问题原因 #

• automap::autofitVariogram 未设计直接接收cutoff的入参，你直接传入的cutoff参数会被函数的省略号（...）捕获但不会传递给内部调用的gstat::variogram，因此不会生效也不会报错
• miscFitOptions列表仅用于控制变异函数模型的拟合逻辑（如 sill 取值范围、是否强制包含块金值等），和半方差样本点对的计算逻辑无关，放入该列表的cutoff也不会被识别

正确配置方式 #

你需要将cutoff参数放入variogram.args列表中传入，该参数的作用就是将用户自定义设置传递给内部执行半方差计算的variogram函数，示例代码如下：

library(automap)
library(gstat)
data(meuse)
coordinates(meuse) = ~x+y

# 正确传入cutoff的写法
autoZn_cut = autofitVariogram(
  formula = log(zinc)~1, 
  input_data = meuse,
  variogram.args = list(cutoff = 1000)
)
# 查看结果验证cutoff生效
summary(autoZn_cut)
plot(autoZn_cut, pch=19, col="black")

运行后查看summary输出的实验变异函数最大距离，可确认截断设置已经生效。

内容的提问来源于stack exchange，提问作者Leprechault