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 <[email protected]>
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
[email protected]
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
[email protected]
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
[email protected]
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.
