如何在 R 中使用应用而不是循环

Question

我正在尝试编写一个计算量非常大的算法，因此我想用一个 apply 函数替换每个循环。但是，我在下面的循环中很累，我想知道是否有人有任何想法。 bc 是一个由 0 和 1 组成的矩阵，ac 是一个由 0 组成的矩阵，与 bc 具有相同的暗度。

for (i in 1:nr){                               
  for (j in 1:nr){
    ac[i,j] <- (bc[i,j]+bc[j,i])/2          
  }
}

Answer 1

我们可以得到bc的总和和bc的转置并除以2。

(bc+t(bc))/2

或者

Reduce(`+`, list(bc, t(bc)))/2

数据

bc <- matrix(1:25, 5,5)
nr <- 5
ac <- matrix(,5,5)