Caros, Preciso de uma ajuda com a funcao abaixo. O que quero fazer e exatamente o que a minha funcao faz, porem sem usar esses for, pq e muito lento e eu preciso avaliar isso para matrized W realmente grande talvez 20000 x 20000, eu fiz algumas tentativas usando expand.grid e apply, mas tive problemas com alocacao de memoria ... Tenho certeza q deve ter uma forma de fazer isso rapido .... Alguma idea ?? Obrigado A = c(1,2,3) W = matrix(rnorm(9),3,3) res = c(1,2,3) cov.prod <- function(A, W, mu3){ soma = 0 n <- dim(W)[1] for(i in 1:n){ for(j in 1:n){ for(k in 1:n){ temp <- A[i]*W[j,k]*r[i]*r[j]*r[k] soma <- soma+temp } } } return(soma) } -- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná