[R-br] Ajuda 3 loop for

[Wagner Bonat]

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á
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140604/b2ade178/attachment.html>