<div dir="ltr"><div><div><div><div><span class=""><span></span></span>Caros,<br><br></div>Preciso de uma ajuda com a funcao abaixo.<br></div><div>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 ??<br>
<br></div><div>Obrigado<br></div><div><br></div>A = c(1,2,3)<br></div>W = matrix(rnorm(9),3,3)<br></div>res = c(1,2,3)<br><div><div><div><div><div><br>cov.prod <- function(A, W, mu3){<br>    soma = 0<br>    n <- dim(W)[1]<br>
    for(i in 1:n){<br>        for(j in 1:n){<br>            for(k in 1:n){<br>                temp <- A[i]*W[j,k]*r[i]*r[j]*r[k]<br>                soma <- soma+temp<br>            }<br>        }<br>    }<br>    return(soma)<br>
}<br><br clear="all"><br>-- <br>Wagner Hugo Bonat<br>LEG - Laboratório de Estatística e Geoinformação<br>UFPR - Universidade Federal do Paraná
</div></div></div></div></div></div>