Самый быстрый способ сделать это двойное суммирование?
Какой самый быстрый способ сделать это суммирование в 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
омега ^ 2 и сигма ^ 2 соответственно.