基于自相关的R语言缺失值插补问题求助（附错误代码）

原代码的问题 #

• 函数Fun_mv无返回值，无法生成结果矩阵
• 错误使用outer函数，outer逐元素传入参数的逻辑和函数内部嵌套循环冲突
• 直接修改全局的AirPassengers对象，易引发意外问题
• 未处理acf计算时的缺失值处理逻辑

正确实现代码 #

# 加载数据集并复制一份避免修改原数据
data(AirPassengers)
ap <- AirPassengers
ap[43] <- NA
ap[100] <- NA

# 定义取值范围
val_range <- 104:622
n <- length(val_range)

# 初始化自相关结果矩阵
acf_matrix <- matrix(nrow = n, ncol = n, dimnames = list(val_range, val_range))

# 遍历所有组合计算滞后1阶自相关系数
for (i in seq_along(val_range)) {
  ap[43] <- val_range[i]
  for (j in seq_along(val_range)) {
    ap[100] <- val_range[j]
    # 计算滞后1阶ACF，关闭绘图、设置缺失值处理规则
    acf_res <- acf(ap, lag.max = 1, plot = FALSE, na.action = na.pass)
    acf_matrix[i, j] <- acf_res$acf[2]
  }
}

# 筛选最大自相关对应的插补组合
max_acf_val <- max(acf_matrix, na.rm = TRUE)
best_pos <- which(acf_matrix == max_acf_val, arr.ind = TRUE)[1, ]
best_43 <- val_range[best_pos[1]]
best_100 <- val_range[best_pos[2]]

# 输出结果
cat("最优插补值：位置43 =", best_43, "，位置100 =", best_100, "\n")
cat("对应最大滞后1阶自相关系数：", max_acf_val, "\n")
# 若需要查看完整自相关结果矩阵，取消下方注释
# print(acf_matrix)

代码说明 #

• 先复制原数据集，避免修改全局变量导致的意外问题
• 使用索引遍历取值范围，避免直接用数值作为矩阵索引引发的维度不匹配问题
• 调用acf时关闭绘图、明确缺失值处理规则，确保计算过程稳定
• 最后定位到最大自相关对应的插补组合，同时生成要求的自相关结果矩阵

注意：104到622共包含519个取值，嵌套循环会执行约27万次计算，运行耗时较长，可结合时间序列的趋势特征缩小取值范围（比如只取对应时段的合理区间值）以提升效率。

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