在R三个变量的剩余总和

Question

我想找到三个变量的剩余总和，我有以下代码：

ans = 0
for (i in 1:dim(x)[1] ){
  for (j in 1:dim(x)[2]){
    ans = ans + (x[i,j] - alpha[i] - beta[j])^2
  }
}

其中x是维数为n的矩阵，m和alpha和beta分别是长度为n，m的向量。

最快的方式是什么？

Answer 1

sum((t(x - alpha) - beta) ^ 2)