[R-br] Eficiênica de algoritmo
Leonard de Assis
assis.leonard em gmail.com
Domingo Setembro 11 23:39:42 BRT 2011
test1=function(){
U=runif(3*10^4)
U=matrix(data=U,nrow=3) #matrix for sums
X=-log(U) #uniform to exponential
X=2* apply(X,2,sum) #sum up to get chi squares
}
system.time(test1())
Código certo.
Repare a última linha
[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
Em 11/09/2011 21:58, Mauro Sznelwar escreveu:
> Eu fiz, mas continua não dando certo, qual o erro?
> > test1=function(){
> + U=runif(3*10^4)
> + U=matrix(data=U,nrow=3) #matrix for sums
> + X=-log(U) #uniform to exponential
> + X=2* apply(X,2,sum) #sum up to get chi squares
> + }
> > system.time(test1)
> usuário sistema decorrido
> 0 0 0
> >
>
>
>
> O test1 e test2 dele são funções, e não variáveis... tu tem que
> criar funções com esses passos e dai sim usar o system.time()
>
>
> []'s
>
> 2011/9/11 Mauro Sznelwar <sznelwar em uol.com.br
> <mailto:sznelwar em uol.com.br>>
>
> Pessoal, eu tenho esta simulação que quero testar a eficiência
> em relação ao rchisq, mas não estou chegando no resultado do
> livro. Alguém poderia dizer o que estou fazendo de errado?
>
> U=runif(3*10^4)
>
> U=matrix(data=U,nrow=3) #matrix for sums
>
> X=-log(U) #uniform to exponential
>
> X=2* apply(X,2,sum) #sum up to get chi squares
>
> test1=X
>
> test2=rchisq(10^4,df=6)
>
> system.time(test1);system.time(test2)
> usuário sistema decorrido
> 0 0 0
> usuário sistema decorrido
> 0 0 0
>
> No livro tem este resultado:
>
> system.time(test1());system.time(test2())
>
> user system elapsed
>
> 0.104 0.000 0.107
>
> user system elapsed
>
> 0.004 0.000 0.004
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto:R-br em listas.c3sl.ufpr.br>
> 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.
>
>
> ------------------------------------------------------------------------
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110911/aebd44bb/attachment.html>
Mais detalhes sobre a lista de discussão R-br