<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.3429" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2><FONT face=CMTT10 size=2>
<P align=justify><FONT color=#0000ff size=3>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?</FONT></P>
<P align=left> U=runif(3*10^4)</P>
<P align=left>U=matrix(data=U,nrow=3) #matrix for sums</P>
<P align=left> X=-log(U) #uniform to exponential</P>
<P align=left> X=2* apply(X,2,sum) #sum up to get chi
squares</P></FONT><FONT face=CMTT10 size=2><FONT face=CMTT10 size=2>
<P>test1=X</P>
<P>test2=rchisq(10^4,df=6)</FONT><FONT face=CMSS10 size=2></P>
<P align=left></FONT>system.time(test1);system.time(test2)<BR>
usuário sistema decorrido
<BR>
0
0 0 <BR>
usuário sistema decorrido
<BR>
0
0 0 <BR></P>
<P align=left><FONT color=#0000ff size=3>No livro tem este resultado:</FONT></P>
<P align=left>system.time(test1());system.time(test2())</P>
<P align=left>user system elapsed</P>
<P align=left>0.104 0.000 0.107</P>
<P align=left>user system elapsed</P>
<P>0.004 0.000 0.004</P></FONT></FONT></DIV></BODY></HTML>