难以在 for 循环中维护 x[i，i]

Question

我一直在

sum <- (abs(x[i, i] - x[i, j])^p) ^ (1/p) 

我希望它是例如x[1,1]，x[2,2]等。所以我认为 for 循环中的x[i,i]可以完成这项工作，但它只给了我 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

Answer 1

每次循环时都会替换总和，sum <- sum + (abs(x[i,i] - x[i,j])^p) ^ (1/p)应该做的。