在for循环中调用x[i,i]异常：代码返回0寻求解决方案

嘿，我来帮你拆解这个问题——你的代码里有两个核心问题，直接导致结果一直是0，咱们一步步理清楚：

先看你写的函数代码：

example <- function (x,p) { 
  sum <- 0 
  for (i in 1:ncol(x)) { 
    for (j in i:nrow(x)){ 
      sum<-(abs(x[i,i] - x[i,j])^p) ^ (1/p) 
    } 
  } 
  return (sum) 
} #x is a matrix

问题1：累加变量被反复覆盖 #

你初始化了sum <- 0，但每次进入内层循环时，都直接把sum赋值为当前i,j对应的计算值，而不是累加。比如第一次内层循环算出一个数，第二次循环就把这个数彻底替换掉了，最后函数返回的只是最后一次内层循环的结果。如果最后一次循环刚好是i=j的情况（主对角线元素自己减自己），那结果自然是0。

问题2：计算逻辑冗余且不符合需求 #

另外，(abs(x[i,i] - x[i,j])^p) ^ (1/p)其实等价于abs(x[i,i] - x[i,j])（只要p>0），这部分嵌套的幂运算完全没必要。而且从你的描述来看，你应该是想计算主对角线元素和同一行其他元素的Lp距离之和，但当前代码既没有累加，也没有正确实现Lp距离的计算逻辑。

修正后的代码示例 #

假设你想要的是：对矩阵的每一行，计算主对角线元素到该行所有其他元素的Lp距离，然后把所有行的距离加起来。这里给你写个修正版：

example <- function(x, p) {
  total_sum <- 0  # 换个变量名，避免和R内置的sum()函数冲突
  n <- nrow(x)
  
  # 确保输入是方阵（行列数一致）
  if (ncol(x) != n) {
    stop("输入的矩阵必须是方阵！")
  }
  
  for (i in 1:n) {
    row_p_sum <- 0
    # 遍历当前行的所有列
    for (j in 1:n) {
      # 跳过对角线元素自己和自己的计算
      if (i != j) {
        row_p_sum <- row_p_sum + abs(x[i, i] - x[i, j])^p
      }
    }
    # 计算当前行的Lp距离，累加到总和里
    total_sum <- total_sum + row_p_sum^(1/p)
  }
  
  return(total_sum)
}

修正点说明 #

• 用total_sum代替sum：避免和R内置函数重名，减少代码混淆。
• 真正实现累加逻辑：每行先计算所有元素与对角线元素的p次幂和，再开1/p次方得到该行的Lp距离，最后把所有行的距离加起来。
• 增加了方阵校验：如果输入不是方阵（行列数不等），直接抛出错误，避免后续循环出错。
• 可选的i != j判断：如果你确实不想计算对角线元素自己和自己的距离，这个判断可以保留；如果需要包含的话，删掉这个判断即可。

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