A maneira mais rápida de fazer esse duplo somatório?
Qual é a maneira mais rápida de fazer esse somatório em R?
Isto é o que eu tenho até agora
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])
}
}
Ondeomega_2
esigma_2
são ômega ^ 2 e sigma ^ 2 respectivamente.