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())

Código certo.
Repare a última linha

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


Em 11/09/2011 21:58, Mauro Sznelwar escreveu:
> 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
>     <mailto: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 <mailto: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.