R语言for循环中Matern模型拟合kappa=50时MSE结果异常差异问题

问题背景 #

使用geoR包的variofit函数拟合Matern模型时，fix.kappa=FALSE自动寻优kappa失败，转而手动遍历kappa取值寻找最小MSE。发现同一kappa=50在不同起始值的遍历序列中，拟合得到的MSE存在差异；但单独运行kappa=50的拟合代码时，结果始终一致。

核心原因 #

variofit函数在未显式指定初始参数（如ini.cov.pars、nugget）时，会默认基于当前环境中的上下文（或前一次拟合的结果）生成初始值。当循环遍历kappa时，前一次拟合的参数会被作为下一次拟合的初始值，导致同一kappa对应的初始参数不同，最终拟合出的MSE存在差异。而单独运行时，初始参数是基于原始变异函数重新计算的，因此结果一致。

解决方案 #

遍历kappa时，显式固定所有拟合参数的初始值，确保每次拟合的起始条件完全一致，避免前一次拟合的影响。

修正后的代码示例 #

假设已生成变异对象vgm：

library(geoR)

# 1. 预先定义固定的初始参数（可基于原始变异函数特征设置）
# 示例：取变异函数的块金值、基台值、半变异函数达到基台的距离作为初始值
init_nugget <- vgm$nugget
init_sill <- vgm$gamma
init_range <- vgm$dist[which.max(vgm$gamma)]

# 2. 定义要遍历的kappa序列
kappa_list <- c(10, 20, 30, 40, 50)

# 3. 遍历kappa，每次用固定初始参数拟合
mse_results <- sapply(kappa_list, function(k) {
  fit <- variofit(
    obj = vgm,
    cov.model = "matern",
    kappa = k,
    fix.kappa = TRUE,
    ini.cov.pars = c(init_nugget, init_sill),
    nugget = init_nugget,
    max.dist = init_range
  )
  fit$MSE
})

# 输出结果
names(mse_results) <- kappa_list
print(mse_results)

额外注意事项 #

• 如果变异函数无明显块金值，可将init_nugget设为0或基于经验值设置；
• 若使用lapply替代sapply，需注意返回结果的格式处理；
• 每次拟合前可通过rm(list = ls(pattern = "fit"))清理环境中前一次的拟合对象，避免潜在的上下文干扰。

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