r - 做这个双重求和的最快方法？

Question

进行此求和的最快方法是什么在R？

这是我到目前为止

ans = 0
for (i in 1:dimx[1]){
  for (j in 1:dimx[2]){
    ans = ans + ((x[i,j] - parameters$mu)^2)/(parameters$omega_2[i]*parameters$sigma_2[j])
  }
}

哪里 omega_2 , 和 sigma_2分别是 omega^2 和 sigma^2。

Answer 1

没有什么花哨:

# sample data
m <- matrix(1:20, 4)
sigma <- 1:ncol(m)
omega <- 1:nrow(m)
mu <- 2

sum(((m - mu) / outer(omega, sigma))^2)