[R-br] Otimizar tempo de processamento

[Tiago Vieira]

Estou tentando rodar o programa de simulação de percolação (rede de bethe), mas esta demandando um tempo de processamento muito alto, pois há preciso de muitas iterações. O computador que estou utilizando tem processador i3 e 4gb de memória. Alguém pode me dar uma dica de como posso diminuir o tempo de processamento? Abaixo segue o código:
 
programa <- function(f,n,m) {
# f: número máximo de filhos
# n: número de gerações
# m: número de iterações
S <- matrix(0,101,2)
for (p in seq(0,1,by=0.01)) {
pp <- 0
for (g in 1:m) {
y <- 1
for (k in 1:n) {
y <- sum(rbinom(y,f,p))
if (y == 0) break
}
if (y > 0) pp <- pp + 1
}
S[p/0.01+1,1] <- p
S[p/0.01+1,2] <- pp/m
}
return(S)
}
 
prog <- programa(2,100,1000)
plot(prog, xlab = "p", ylab = "Probabilidade de Percolação")
 
Desde já agradeço pela atenção!
 
Tiago Vieira.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110818/161b27bc/attachment.html>