
11 Set
2011
11 Set
'11
22:12
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