[R-br] Eficiênica de algoritmo

Mauro Sznelwar sznelwar em uol.com.br
Domingo Setembro 11 21:58:27 BRT 2011


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>

    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
    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/37ded643/attachment.html>


Mais detalhes sobre a lista de discussão R-br