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á
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
