¿La forma más rápida de hacer esta doble suma?
¿Cuál es la forma más rápida de hacer este resumen?
Esto es lo que tengo hasta ahora
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])
}
}
dóndeomega_2
ysigma_2
son omega ^ 2 y sigma ^ 2 respectivamente.