[R-br] Eficiênica de algoritmo

[Mauro Sznelwar]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110911/ff1aa520/attachment.html>