R语言嵌套循环迭代后变量值重置问题求助

你的代码瞬间执行但无任何有效输出，核心问题是循环范围的写法错误：

• for(i in 252) 并非让i遍历1到252的所有值，而是仅循环一次，i固定为252；同理for(j in 10000)和for(k in 10000)也只执行单次迭代，而非遍历1到10000的序列。
• 若252超出countries$avdist的索引范围，或计算时的坐标索引超出points的行数，就会导致变量完全无变化。

先把循环范围改为序列形式（如1:252），同时优化索引计算避免重复运算：

# 确保预声明向量的初始值正确
subtotal <- numeric(10000)
total <- numeric(10000)

for(i in 1:252){
  start_idx <- (i-1)*10000  # 提前计算当前组的起始索引，减少重复计算
  for(j in 1:10000){
    j_coord <- c(points$POINT_X[start_idx + j], points$POINT_Y[start_idx + j])
    for(k in 1:10000){
      k_coord <- c(points$POINT_X[start_idx + k], points$POINT_Y[start_idx + k])
      subtotal[k] <- distGeo(j_coord, k_coord)
    }
    total[j] <- mean(subtotal)
  }
  countries$avdist[i] <- mean(total)
}

你原计划的三重循环总次数达252×10000×10000=2.52e10次，这在R中几乎不可能跑完（哪怕每秒执行100万次，也需要7000+小时），必须用向量化运算大幅优化：
利用geosphere包的distm函数直接生成一组内所有点的距离矩阵，再批量计算均值：

library(geosphere)

for(i in 1:252){
  start_idx <- (i-1)*10000
  end_idx <- start_idx + 10000
  # 提取当前组的坐标子集
  current_points <- points[start_idx:end_idx, c("POINT_X", "POINT_Y")]
  # 一次性生成所有点对的距离矩阵
  dist_matrix <- distm(current_points, current_points, fun = distGeo)
  # 先算每行均值（每个点到组内其他点的平均距离），再算组整体均值
  row_means <- rowMeans(dist_matrix)
  countries$avdist[i] <- mean(row_means)
}

这种方式将循环次数从2.52e10降至252次，效率提升数个数量级，完全具备可执行性。

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