Трудность в поддержании x [i, i] в цикле for

Question

У меня возникла проблема в

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

Я хочу, чтобы это было, например, x[1,1], x[2,2] и т. Д. Поэтому я подумал, что x[i,i] в цикле for сделает эту работу, но это только возвращает мне 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) должен делать.