[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com

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@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@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@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@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.